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Introduction 



Quantum field theory (QFT) is an essential tool in order to understand phenom- 
ena occuring from atomic (10~^° m) to the Planck scale (10~^^ m). It provides a 
theoretical framework to calculate scattering cross sections, particle lifetimes, and 
other observables of processes of the quantum world. Its most powerful strength 
is demonstrated by the fact that every elementary interaction known in particle 
physics can be described by a corresponding relativistic quantum field theory. The 
main reason of the introduction of the field description is that the processes of par- 
ticle physics changing the particle number can not be described by the simple ap- 
plication of quantum mechanics (e.g. see the Klein-paradox [1]). Furthermore, QFT 
combines quantum mechanics and special relativity naturally: it explains the relation 
between spin and statistics, and solves the causality problem of (relativistic) quan- 
tum mechanics. Quantum field theories also play an important role in statistical- and 
condensed-matter physics, e.g. in the description of critical phenomena and quantum 
phase transitions. It is extensively used in atomic and nuclear physics as well. 

In field theory we have infinitely many degrees of freedom: one assigns dynamical 
variables (typically scalars, vectors, spinors and/or matrices) to every point of the 
three (or more generally d) dimensional space and follow for the time evolution 
of them. The classical dynamics is determined through the Lagrange formalism as 
a generalization of classical mechanics. In the quantum version of the theory, the 
variables are considered to be linear operators acting on the Hilbert space of the 
quantum states in question (in QFT this is the so-called Fock space). The key of 
the transition to QFT from the classical version is the so-called canonical equal-time 
commutation relation. 

In order to obtain measurable quantities, one has to compute n-point functions, 
which arc vacuum expectation values of the time ordered product of the fields. The 
most successful method for obtaining them is perturbation theory, in which every 
quantity is expressed through power series (pcrturbative series) of small parame- 
ters of the theory, i.e. the coupling constants. Although general analyses proved 
that pcrturbative series are asymptotic [3], therefore the convergence radius is zero, 
very valuable results can be obtained within this framework, e.g. the perturbatively 
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solved quantum electrodynamics (QED) has extremely accurate predictive power. 
For example its prediction for the anomalous magnetic moment of the electron agrees 
with the experimentally measured value to more than 10 significant digits (this is 
the most accurately verified prediction in the history of physics) [1]. However, for 
strongly coupled theories perturbation theory completely breaks down; the series 
do not seem to show even practical convergence at all. At low enough energy, this 
is the case of the fundamental theory of the strong interaction, the quantum chro- 
modynamics (QCD) [2]. One of the most challenging task of contemporary particle 
physics is to solve QCD in different setups relevant for phenomena from accelerator 
experiments to internal structure of compact stellar objects. 

The most productive and mathematically well controlled method to attack the 
problem of large couplings is lattice field theory. It is a fully non-perturbative ap- 
proach, in which the degrees of freedom are reduced in a way that the quantum 
equations of motion are solved in a lattice, instead of the continuum formulation 
[5]. The method typically requires very huge computational resources, nevertheless, 
it seems to be the most competitive method in order to describe non-perturbative 
features of the spectra of strongly coupled theories. 

A more traditional alternative way to deal with the problem is the application of 
functional methods. The best known technical framework is offered by the infinite 
set of Dyson-Schwinger equations. In the first chapter, we will discuss the main 
properties of this approach, in particular a variant called two-particle-irreducible 
(2P1) formalism, which is the central topic of this thesis. 

Before we proceed to review the literature on 2P1 formalism, let us shortly discuss 
a central problem of every quantum field theory. QFT in its original formulation is 
ill-defined: divergences are obtained even at the lowest orders of perturbation theory. 
These divergences come from the ultra-violet regime of the momenta, or in other 
words from spacctimc points which are very close to each other. This implies that 
QFT at very short distances must break down, which gives no surprise; we should 
not expect that in an arbitrary small neighborhood of a given point infinitely many 
physically relevant degrees of freedom can appear. We should think of QFT as an 
effective model of some more fundamental theory (e.g. string theory), and expect 
that the details of the dynamics at this scale are not relevant for obtaining results 
to a much larger scale. This expectation leads us to the idea of renormalization. 

In order to obtain well-defined quantities, first we have to regularize QFT. The 
most natural regularization is to build up equations on a lattice, or the use of a 
cutoff in the momenta. On the regularized theory, we can perform the so-called 
renormalization programme. This means that our lack of knowledge of the short 
range dynamics is hidden into a finite number of physical quantities determined 
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by appropriate measurements or fixed by some conditions. If tliis can be done in a 
consistent way, wliere the cutoff can really be thought as a large quantity, therefore 
formally(!) can be sent to infinity with obtaining finite results, then the theory is 
said to be renormalizable. 

Renormalization and renormalizability have an extended literature and in per- 
turbation theory it is worked out in great details. Usually the way we treat the 
theory is to think of the model parameters as cutoff (defined as the highest al- 
lowed momentum) dependent objects, which can be separated into renormalized 
parts and counterterms. In perturbation theory we treat also the counterterms as 
perturbations, which can be determined from the requirement of the cancellation 
of the appearing divergences. Pioneers of perturbative renormalizaton with rigorous 
proofs built upon Feynman's diagram technique were Bogoliubov, Parasiuk, Hepp 
and Zimmermann [8, 9], together with Weinberg [10]. A very detailed description of 
their work and the general theory of perturbative renormalization can be found in 
[3]. 

2PI formalism generalizes the idea of the IPI effective action [4] . The so-called 2PI 
effective action has two types of variables: mean fields and propagators. Stationary 
conditions of this action lead to self-consistent equations for these quantities. The 
solution induces a partial resummation of the perturbative series (i.e. resums certain 
infinite classes of Feynman diagrams), therefore it is a rather nontrivial question, how 
the perturbative renormalizability is carried over into certain 2PI approximations. 
One of the goals of this thesis is to show various methods, which are capable to 
formulate 2PI renormalizability. 

The formalism was developed in the 1960s for non-relativistic field theories of 

condensed matter physics [11], and was generalized to relativistic theories in the 
1970s [15] by Cornwall, Jackiw and Tomboulis. Referring to the latter paper, the 
method is sometimes called CJT formalism. The ^-derivable approximation appel- 
lation is also used with reference to the paper [14], but these mean the very same 
formulation. Although the technique is known for a long time, the applications have 
become an especially active research topic in the last 10 years placing emphasis on 
non-equilibrium simulations. This is due to the fact, that 2P1 formalism is the only 
known solution to the so-called secular time-evolution problem. Considering scalar 
fields for instance, the simplest non-local contribution to the 2P1 effective action 
leads to non-secular time evolution showing thermalization at late times [21]. Same 
observations were made in theories coupling scalars to fermions [22]. Some of the 
non-perturbative effects of 2PI approximations can be captured by combining the 
approach with the large-N expansion. In out-of-equilibrium phenomena, the 2PI 1/N 
expansion at next-to-leading order proved to be very valuable in studying various 
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problems, including thermalization [23, 24, 25] and cosmological preheating [22, 26], 
non-thermal fixed points [27] and decoherence [28]. The lowest ordered 2P1 approx- 
imations in equilibrium also show good convergence of thcrmodynamical quantities, 
such as the pressure [29, 30]. Discussions of finite temperature gauge theories can 
be found in [31, 32]. There were attempts to determine transport coefficients in 
[33, 34, 35]. Phenomenological studies have also appeared using various 2P1 approx- 
imations, see [36, 37, 38, 39]. Very recently promising applications appeared for the 
evaluation of critical exponents for continuous phase transition [40]. 

The renormalization of approximations obtained from the 2PI technique was 
investigated first for scalar theories at finite temperature in real time formalism [41, 
42] and also in imaginary time formalism [43, 44]. The method was then extended to 
fermions [45] and to QED [46] . Renormalization with non- vanishing field expectation 
values and therefore the determination of the effective potential was discussed in 
[44, 47, 48, 49, 50]. In case of theories with more complicated global symmetries 
renormalization techniques were developed in [51]. Although many papers discussed 
2P1 renormalization and renormalizability, there are still some open questions, in 
particular the standardized qualification of its numerical realizations is missing. 

The present thesis deals with quantum field theories in equilibrium. Its goal is 
to map and solve some problems of renormalization both from an analytic and a 
numerical point of view, together with applications of some newly developed tech- 
niques to scalar quantum field theories. It is important to stress that publications 
discussing the accuracy of numerical realizations of 2P1 renormalization in a repro- 
ducible manner are virtually missing from the literature. One of the aims of this 
thesis is to fill the gap caused by this absence of interest. 

The greater part of the models examined here are used as effective theories of 
strong interactions. With this, the thesis would like to contribute to developing the 
treatments of models of phenomenological importance beyond perturbation theory, 
which - as already discussed - is crucial when large coupling constants appear. It is 
important to stress that this thesis has no task to do phenomenology. We wish to 
search for appropriate methods applicable for phenomenological investigations. 

The structure of the thesis is as follows. In Chapter 1, we provide a discussion 
of the functional methods of QFT. The task of this part is to present pedagogically 
the usual generator functionals of the Green's functions in order and end up at 
2PI formalism. As an intermediate step, we will also discuss the Dyson-Schwinger 
equations and establish connections to two-particle-irreducibility. In Chapter 2, we 
introduce the models investigated in the thesis. At this point we will discuss some 
fundamental aspects of the large-N technique, which will also be a central subject 
of our investigations during the next chapters. In Chapter 3, we deal with the 0(N) 
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model in 2PI formalism at next-to-leading order of the large-N expansion. We will 
use the so-called auxiliary field formulation, and demonstrate its renormalizability 
with constructing appropriate counterterms explicitly. The equivalence between the 
original and the auxiliary field formulation will also be discussed. The task of Chap- 
ter 4 is to present a careful numerical study of 2P1 renormalization, in which the 
one component 0^ theory was employed. The first part of this chapter nevertheless 
deals with a theoretical issue: we show that the minimal subtraction procedure used 
in Chapter 3 is equivalent to imposing appropriate renormalization conditions on 
some quantities. This will allow us in one hand to obtain explicitly finite 2P1 equa- 
tions, and on the other hand the opportunity to compare the convergence of the 
numerical solutions of these equations with the ones containing counterterms. We 
will discuss the features of various numerical algorithms and produce very accurate 
results. In Chapter 5 we turn again to a phenomenologically more important the- 
ory, the U{N) X U{N) symmetric meson model. We will present an approximate 
large-N solution in the broken phase, which turns out to be renormalizable and ful- 
fills Goldstone's theorem. Various symmetry breaking patterns will be investigated, 
and the renormahzed effective potential will be constructed. We will investigate the 
vacuum structure of the theory which reveals new extrema of the effective potential 
not observed in previous perturbative computations. 
The thesis is based on the following four pubhcations: 

• G. Fejos, A. Patkos, Zs. Szep: Renormalized effective actions for the 0(N) 
model at next-to-leading order of the 1/N expansion 

Phys. Rev. D80, 025015 (2009), arXiv:0902.0473 [hep-ph] 

• G. Fejos, A. Patkos: A renormalized large N solution of the U(N)x U(N) linear 
sigma model in the broken symmetry phase 

Phys. Rev. D82, 045011 (2010), arXiv:1005.1382 [hep-ph] 

• G. Fejos, A. Patkos: Spontaneously broken ground states of the U(N)x U(N) 
linear sigma model at large N 

Phys. Rev. D84, 036001 (2011), arXiv: 1103.4799 [hep-ph] 

• G. Fejos, Zs. Szep: Broken symmetry phase solution of the 0^ model at two-loop 
level of the ^-derivable approximation 

Phys. Rev. D84, 056001 (2011), arXiv:1105.4124 [hep-ph] 
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Chapter 1 



Functional techniques in quantum 
field theory 

In this chapter we review the standard functional techniques of quantum field 
theory. We will go through the ordinary generator functionals and show explicitly 
the properties of these quantities. The motivation for such a summary is that our 
goal at the end of this chapter is to arrive at the so-called two-particlc-irrcducible 
(2PI) formalism, which is a generalization of the standard one-particle-irreducible 
(IPl) formulation. The best way to achieve this is to build up the usual generator 
functionals from the beginning to see the correspondence between the different ap- 
proaches. In this introductory chapter we present the functional techniques in a less 
formal but nevertheless expressive way. The discussion is built in a way from which 
the 2P1 formalism can be obtained very naturally. Derivations will be presented 
mainly graphically, using the language of the Fey nman- diagram technique. For a 
greater transparency we illustrate the general line of thinking within the one com- 
ponent 0^ theory in the symmetric phase, nevertheless all steps can be performed 
easily in arbitrary theories. 

In this thesis we work in h = c = 1 units and this chapter will contain calculations 
only in Minkowski space. The metric is defined as x'^ = x^x^ = Xq — x^. In this 
thesis only scalar quantum fields will appear in 3 -|- 1 space-time dimensions, with 
Lagrangians parametrised as 

C[d,4>, 0] = -\(t>{d,d^ + - (1-1) 

where is a multicomponent c- number scalar field and describes the self- 

interaction. The n-point Green function is defined as a time-ordered vacuum expec- 
tation value of Heisenberg field operators 

Gn{xi,X2,...,Xn) = {0\T {^{xi)^{x2) . ..4>{Xn))\0) . (1.2) 

7 
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In functional integral formalism these functions can be obtained as [4] 

where refers to the action of the configuration (f){x): S[4>] = J (i'^x£[(9^0, 0]. 

The Green functions are not measurable quantities, but with the help of the so- 
called reduction formulae (Lehmann, Symanzik and Zimmermann, 1955), they can 
be turned into physical transition amplitudes, see details in [1]. Due to this result, 
one can prove that the Green functions incorporate all information which can be 
extracted from a quantum field theory. Therefore, it is worthwhile to concentrate 
on building up equations for these quantities, which we eventually do in the next 
sections: the corresponding Dyson-Schwinger and 2PI equations will be derived. 
Before turning to this task, let us discuss the generator functionals of the Green 
functions. The properties of these functionals summarized here are quite general 
and can be obtained in arbitrary theories. 



1.1 Generator functionals 

The generator functional of the scalar Green functions is a functional of a source 
J, which is an arbitrary c-number function: 

^[J] = ^ / d'^Xi...d^XnJ{xi)...J{Xn)Gn{xi...Xn), (1.4) 
n=0 ^' ^ 

therefore the Green functions can be obtained by functional differentiation as 



Gn[Xi,X2, ...,Xn) = (-^) " 



j=o "-^ 



(1.5) 

j=o 



SJ{Xi)...SJ{Xn) 

We will often not write down explicitly the spacetime arguments and use the short- 
hand notation as shown by the second equality. Z[J] can be represented with func- 
tional integrals: 

as it can be easily seen by expanding exp ii J Jcj)) in Taylor scries. To calculate Z[J] 
and eventually the Green functions we split the action into two parts: 

5[0] = 5o[0]+5/[0], (1.7) 
where 5*0 refers to the gaussian part (quadratic in 0), while Sj to the interaction: 

So[4>] d^^(t)ix){d^d^ + m')0(x), Si[(f)] = - y" d^xU[(l)ix)]. (1.8) 
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We write Z[J] as 

/ r'0e^(^M+/-^'^)/ / p^e^^ol'^] _ <exp{iSi + i J J(l))> 
^^^^ ~ J V(l)e'Sl^y J V(l)e'3o[<t>] = <exp {iSj)> ' ^^'^^ 

where < . > refers to the operation of the gaussian averaging. The evaluation of the 
numerator and the denominator can be made by expanding the exponentials in the 
brackets, since Green functions of the gaussian theory can be easily calculated [1]. 
Here we implicitly assumed that > in order to have a positive definite quadratic 
form for the gaussian part, therefore the expansion makes sense. For each term in 
the expansion one can associate a diagram with defining model specific Feynman 
rules. For ?7[0] = A0^/4! of a one component field these are the following [1]: 




Af{x - y) 



—iX I d z 



i I d^zJ{z), 



where Ap{x — y) is the Feynman (or perturbative) propagator: 

/j4l ■ ik{x-y) 

Every diagram must be divided by a symmetry factor, which is defined as the di- 
mension of the graph's symmetry group. Note that the 1-point vertex appears only 
in the expansion of the numerator, due to the term proportional to the source. In 
the averages every possible diagram without external legs appear, these are the so- 
called vacuum graphs. Using the Feynman rules (keeping also the symmetry factors 
in mind), in 0^ theory the first few terms of the denominator are 




We can observe that some terms contain disconnected diagrams, which offers a 
possibility for a simplification of the expression. Let V be the set of the internally 
connected diagrams: 
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We can check that every term in (1.11) can be written in the following form: 

i 

where Vi is the ith element of V and equals to the number of occurrences of this 
connected diagram in the term in question (the factorials are symmetry factors of 
the interchanging of disconnected diagrams). The denominator of (1.9) is therefore 



(1.12) 



where the summation goes through all ordered sets {ni,n2, ...}. It is easy to show 
that this expression can be factored: 



Til n2 na i rii 

We recognize the definition of the exponential function and obtain 



<e'^i>= JJexp^i =exp^Vi = exp (W + + O ^ ^ ^ ...).(1.14) 



This is the linked cluster theorem, which states that the sum of all possible vacuum 
diagrams is equal to the exponential of the sum of all the connected ones. The 
numerator of (1.9) can be evaluated the very same way. We have to keep in mind 
that there we also have a 1-point vertex, therefore in addition to the diagrams of V, 
we have new vacuum graphs. Let V be the extended set of the vacuum diagrams: 



V 




Note that the crosses do not correspond to external legs, these are 1-point vertices, 
therefore every diagram of V is indeed a vacuum graph. Following the method de- 
scribed above we arrive at 



<e 



iSi+i f J4> 



>= exp 



+ 




+ 



+ 




+ ... 



+ 



XI 



+ 



+ 



+ 



on 



+ 



+ ... 




+ 




+ 




...).(1.15) 
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Since Z[J] is the quotient of (1.15) and (1.14), only diagrams involving the source 
remain [1]: 



Z[J] = exp 8 ® + 8 — Q — » + 8 — Q — » + ( »00 (g , + ^ Q 



+ X + X A=A ^'-''^ 

If now we are interested in the value of the various n-point functions, then we 
should re-expand the exponential function and read off the coefficients of the power 
series in the source J. We immediately recognize that in our concrete example only 
even number of sources could appear in every term, therefore when we take the 
source to zero only Green functions of even variables are non-zero. (In the broken 
phase of 0^ theory also Green functions of odd variables are present.) For example 
the 2- and 4-point functions read as 



G2(xi,X2) = . .+ O + O + OO + . Q ....,(1.17) 



Gi{xi, X2, ^3, 0:4) — 





+ ... (1.18) 



The crosses from the endpoints disappeared, they indicated the presence of the 
sources (therefore neither multiplicative factors nor integration is needed correspond- 
ing to these endpoints). This means that these diagrams are not vacuum graphs 
anymore but have 2 and 4 external legs, respectively. 

We recognize that there are lots of disconnected diagrams appearing in the series 
of the Green functions. (We note that in the symmetric phase of the 0'^ theory G2 
contains only connected diagrams - see (1.17). This would not be the case in the 
broken phase, since then there would be a 1-point vertex in the original Lagrangian, 
which could also produce nonzero contribution to Gi.) Recalling (1.16) it is possible 
to build up a generator functional which generates only the connected pieces, which 
provides great simplifications in further calculations. We introduce W[J] as 



iW[J] = log Z[J]. 



(1.19) 
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Factoring an i is just a convention, which will prove to be useful later. Using (1.16) 
we have 



= ® ® + ^ — Q — + 8 — Q — g) + ^ Q O (g ) + ^ Q (g )... 

+ V + X +■••+ +- (1-20) 






If we think of W[J] as a power series of the source: 

iW[J] = — / d'^Xi...d'^XnJ{Xi)...J{Xn)Wn{Xi, Xa, (1.21) 

it is obvious that the Wn quantities refer to sums of connected diagrams. These are 
to be called as connected Green functions or cumulants and can be calculated as 



(1.22) 

j=o 



In our 0^ example due to the absence of nonzero Gi, we have Wi = and W2 = G2, 
however W4 and further Green functions are greatly simplified: 



W2{xi,X2) = . . + O + O + 00 + Q ...,(1.23) 



w,{x,,x2,x,,x,)= + y-y +... \x + X/ +...(1.24) 







It turns out that even more simplification can be made, which leads us to the 
idea of the quantum effective action. Let us introduce the so-called semi- classical or 
mean-field (f) as 

which is the vacuum expectation value of in the presence of the source J. We 
assume that this relation can be inverted uniquely: for a given J only one exists 
and vice-versa. The quantum effective action is the functional Legendre-transform 
of the generator of the connected Green functions W[J] [4]: 

r[0] = W[J] - [ J(f). (1.26) 
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Using (1.25) we have the usual relation 

5r[4>] r 5W5J f 5 J 



f5W5J fSJj ^ ^ 



50 

Let us try to calculate the diagrams contributing to this quantity. We have 



^ ^ = r^. -.ui . (1-28) 



therefore 



where in the second equality wc performed the following change of the integration 
variable in the numerator: 0^0 + 0. The shifted action in the numerator can be 
expanded around 

+ 0] = sw\ + + \s"m^ + (1-30) 

which gives new vertices. 

We are now equipped to calculate such a quantity like (1.29). As we have already 
seen, the denominator is the exponential of all possible connected vacuum diagrams 
with couplings of the original Lagrangian. The numerator contains the same type 
of graphs but with vertices of the shifted Lagrangian extended with an additional 
term for the 1-point vertex coming from the source. Diagrams of the denominator 
cancel the diagrams of the numerator which do not contain at least one of the new 
vertices or the source. This is a consequence of the normalization Z[0] — 1, which 
implies r[0] = 0. 

Before starting to calculate the remaining diagrams by brute force, let us point 
out an important simplification taking place. The key observation is that the 1-point 
function of the shifted action is zero by construction [7] : 



This gives a great simplicity of the structure of V[(f)\. Consider a diagram of the 
numerator in which two subdiagrams 7*^^-' and 7*^^-' are connected with one single 
fine. This is to be called 1-particle-reducible (IPR). Then find every possible diagram 



14 



Chapter 1. Functional techniques in quantum field theory 




Figure 1.1: 1-particle-reducible diagrams do not count in the effective action. 

where 7(2) is replaced with some other subdiagram. In the sum of these diagrams, 
the exact 1-point function will appear hanging on 7^^^. Since we saw that the 1- 
point function is identically zero, we may forget about diagrams with the previously 
indicated property (see Figure 1.1). In other words, we may omit all IPR diagrams 
and consider only one-particle- irreducible (IPI) ones. These do not fall into two 
parts if one internal line is cut. The expressions corresponding to these diagrams do 
not refer to the 1-point vertices at all, therefore the source does not appear in r[0], 
as it should. We may write 



(1.32) 

ipi 



When we are interested in the physical case (i.e. J — )■ 0), from (1.27) we have 



5(t) 



0, (1.33) 



which means that the semi-classical field is an extremal configuration of F. This is the 
reason for calling it quantum effective action. In agreement with the expectations, 
after expanding the shifted action in (1.32) around S'[0] and ignoring fluctuations, 
F[0] becomes the classical action: F[0] ^ 

In 0^ theory the shifted action introduces the following new 2 and 3-point vertices 
(the 1-point vertex is irrelevant due to the previous analysis): 




Note, that some sticks appeared at the end of the dashed lines. These refer to 
the multiplication(s) of and to the integration over spacetime. The first few terms 
in the series of F[0] turns out to be: 
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+ 




(1.34) 




We may think of r[0] as a power series of (p: 

r[0]=^^ J d^Xi...d'^XnTn{Xl,...,Xn)(f){Xl)...(}){Xr^ 



;i.35) 



where the r„ coefficients are the so-called proper vertices. From (1.34) we can read 
for example that 



^2{xi,X2) = -{d +m )6{xi - X2) - i 



.0. 



+ 




;i.36) 



r4(xi, X2, X3, X4) = \5{Xi — X2)S{Xi — Xz)5{Xi — X4) — i 



:<!>: 



+ 



+ 



+ 




+ 



1.37) 



Proper vertices with odd number of variables are zero in the symmetric phase of the 
(j)^ theory. Note that diagrams in (1.36) and (1.37) do not have sticks on their "legs", 
as the multiplicative factors (their integrals more precisely) were not included in 
the definition of the proper vertices. One also has to note that the factorials in (1.35) 
were "eaten up" by the symmetry factors of diagrams still containing these sticks 
(i.e. diagrams in (1.34)). 

r2 is strongly related to the connected 2-point function, since 



z) 6 



5J{y) 6J{y) <50(z) J, 6J{y)6J{z) S^{z) ' 
which applied to the function J(x) gives 

S^W[J] 5J{x) f 5^W[J] S^T[ 



:i.38) 



S{x - y) 



6J{y)6J{z) 5(l){z) 



6J{y)6J{z) 6^{z)6(P{x) 



;i.39) 
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This shows that in functional sense 



:i.40) 



(5J2 

which imphes that 

iiy-i = T2. (1.41) 

The squared bracket of (1.36) including the factor i is called the self-energy, and is 
denoted by S(xi, X2). With the help of (1.40) we arrive at the Dyson- equation of the 
propagator 

iW2 \xi,X2) = -{d^ + m2)5(xi - X2) - S(xi, X2), (1.42) 

to which we will refer several times in the forthcoming sections. 

It is possible to organize the series of F in terms of the number of loops. When 
one works in ordinary units, it turns out to be a series in h. Introducing the notations 

r[0] = ^[0] + r,[0] + r,2[0] + ..., (1.43) 

the leading order piece (beyond the classical value) F/j can be resummed in a rather 
compact form, in which one has to take into account every 1-loop diagram: 

<\ <\ <\ 

Taking into account the symmetry factors, the sum of these diagrams is 



2 V 2 
1/ .A\2 



X 



+ 4 ( ~ ' 2 ) j j ^^^""^ y)^p(y^ x)^\x)^^{y) 

+ ^(-4)V / [^F{x,y)AF{y,z)AF{z,x)^\x)4>\yW{z) + ..., (1.45) 



1/ .A\3 

'x J y J z 

where the integrals are over the whole space-time with respect to the corresponding 
letters. One can identify the series of the logarithm function: 



iVn = l^Tf2l[^^'A^Y = -^^Tlog{l + ^^'^A^). (1.46) 



n=l 



1„ , A^o 



The trace has to be taken in functional sense. Recalling that U" {(j)) = l^^, for F we 
have 

F[0] = S[4)] + 2 '^^ [^ + U {(j))iAp^ — i{ diagrams with at least 2-loops).(1.47) 



1.2. Dyson- Schwinger equations 
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This equation serves as a general formula applicable in arbitrary scalar theories, not 
just in the 0^ model. In the case when is a multicomponent field, the trace has to 
be taken over the internal indices as well. 

From the proper vertices every Green function can be built up, therefore it is 
sufficient to determine the r„ functions to solve fully a quantum field theory. We 
wish to mention that it is often more convenient to work with the Fourier transform 
of the proper vertices, denoted by F^. Due to the translation invariance, one can 
always factor out a delta function, therefore F„(pi, ...,p„) can be always written as 

fn(pi, ■■■,Pn) = (27r)^5(pi + ... +p„)r,(pi, (1.48) 



1.2 Dyson-Schwinger equations 

As already announced at the beginning of the chapter, in order to obtain all 
the information contained in a quantum field theory, it is sufficient to determine all 
the Green functions. In practice we need to calculate the r„ proper vertex functions 
generated by the effective action. In this section we derive an infinite tower of coupled 
equations from which any r„ can be obtained. We begin with recalling the translation 
invariance of the measure of the functional integral [4] : 

= j -D^j^i...). (1.49) 
Let us apply this identity to the function e**^'^['^]+-^'''^^: 

= / p0^e^(5M+/J0) ^ [ Ve^m\+!J<t>)(i^_^+u]. (1.50) 
J Sep J \ dcp J 

Since under the integral a multiplication with is equivalent to differentiating with 
respect to the source, we may write the formal expression 







5S 



+ jU^sm+!J^). (1.51) 



Now neither of the terms in the bracket contain therefore it can be taken outside 
the integral. Dividing (1.51) by J T>(f)e'^^^'^^ we obtain 



= 



6S 



+ JjZ[J]. (1.52) 

Since Z[J] — expil^[J], this equation is equivalent to 

e^^[-^l + J. (1.53) 
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Let T[J] be a test functional, then 

S_ 

This means that (1.53) can be also written as 

\54> 

where 1 is the unit functional. RecaUing (1.25) and (1.27), we have 

5r[0] 5S 



+ J 



, SWIJ] . s 
6 J ^SJ 



)l = 0. 



:i.55) 



:i.56) 



5(t) 54) 

Using (1.38) and (1.40) we arrive at the generator equation of the Dyson-Schwinger 
equations for r„: 

5r[0] 6S 



5(f) 5(1) 



[1.57) 



which shows how to compute the derivatives of the effective action and builds up a 
tower of coupled equations for the proper vertices. 

In (p'^ theory with the classical potential U{(p) = A0^/4!, the first two derivatives 

are: 



5m 

5(f){x) 



5^r 



5(f)(x)5(l)(y) 



= -(d^ + m^ + -{^\x) + 3W2(x,x))')^(x) 

- / W2{x,u)W2{x,v)W2{x,w)r3{u,V,w) 
J uvw 

= -(d'^ + m'^ + ^{(j)'^{x) + W2{x, x))^5{x - y) 

- i^4>{x) [ W2{x,u)W2{x,v)T3{u,v,y) 



:i.58) 



+ 



A 



W2{x, s)W2{x, u)W2ix, v)rs{u, V, w)W2{w, t)rs{s, t, y) 



stuvw 



— I 



A 
3! 



W2{x, u)W2{x, v)W2{x, w)Ti{u, V, w, y). 



:i.59) 



For a given n-point proper vertex, higher V^n functions always appear in their Dyson- 
Schwinger equations, therefore practically at some point we have to make an approx- 
imation and truncate the tower of the coupled equations. When the configuration 
is set to the physical one (i.e. J — > 0), we know that in the symmetric phase of 
0^ theory only proper vertices with even variables are nonzero. In this case (1.58) 
becomes trivial, and for (1.59) we have (note that r2 = iW2^) 

T2{x,y) = -{d^ + + ^W2ix,x))5{x-y) 



— I 



3! 



W2{x, u)W2{x, v)W2{x, w)Vi{u, V, w, y). (1.60) 



1.3. 2PI formalism 
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In agreement with the previously introduced Feynman rules, diagrammatically it 
can be written as 



;i.61) 



where the thick lines do not refer to Feynman propagators, but to full connected 
2-point functions {W2) and in the second term of the squared bracket the circle 
represents —iV/^ (one should keep also the symmetry factors in mind). (1.61) is very 
similar to (1.36). In fact, this is the same equation, but in a different form. In (1.61), 
diagrams of (1.36) are resummed into W2 and If we insert the perturbative 
expansion of W2 and r4 into the right hand side of (1.61), we obtain (1.36). In the 
next section, we will derive this equation in a third form using 2PI formalism. 

What makes the tower of Dyson- Schwinger equations (and therefore the form 
(1.60) of the Dyson-equation) in particular important is that they are valid and 
make sense also non-perturbatively (i.e. without referring to any Feynman diagrams). 
However, closing the hierarchy using appropriate approximation(s) typically lead to 
partially resummed perturbative series of the vertices yet to determine. This will also 
be a key feature of the 2PI formalism, which is introduced in the next section. Due 
to this, Dyson- Schwinger and 2PI formalisms are related in certain approximations 
(this will be highlighted later). 



1.3 2PI formalism 

Two-particle-irreducible (2PI) formalism generalizes the idea of the effective ac- 
tion: it treats also the propagator as a variable. The so-called 2PI effective action 
can be obtained easily in the framework showed before. 

We introduce into the usual Z generator functional of the Green functions a 
bifocal source K{x,y) with respect to the product of the fields [15]: 

rmei(s[0]+/J0+|//i^</'</') 
Z[J,K] = ^ — ^ — „ ^ . (1.62) 

The quadratic term including K can be viewed formally as an additional "mass 
term" , therefore our previous analysis can be directly applied. This means that 

iW[J,K]=\ogZ[J,K] (1.63) 

contains the connected diagrams and in 

TkW\ = W[J,K]- J J0 (1.64) 
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we find the IPI diagrams, where the (f) semi-classical field in the presence of sources 
J and K is defined as before: 

SW[J,K] 

-JJ^=m- (1-65) 

Since we have two sources, it is possible to perform one more Legendre-transformation. 
Therefore in addition we introduce G{x, y) (which will refer to the connected 2-point 
function) with defining the derivative with respect to X as 

f|f = i(»)+G(.,.)). (1.66) 

The J and K dependence of and G are not indicated explicitly, as usual. The 2PI 
effective action is the double Legendre-transform oiW[J,K]: 

mG]^W[J,K]- j J4>-\j j K{G + U)- (1-67) 
Using (1.65) and (1.66) we have 

--J{x)- K{x,y)(j){y), — r ^ -K{x,y). (1.68) 



Scpix) Jy ' ' SG{x,y) 

In the physical case (i.e. J — > 0, X — > 0) and G are determined from the stationary 
conditions 

G] 5r[4>, G] 

rl( X = 0, — r = 0. (1.69) 

50(x) oG{x,y) 

Similarly to the procedure of obtaining the IPI effective action, we have 
Performing the change of variable — >■ + 0, 

Wc may now try to calculate both the numerator and the denominator with the 
method introduced in the beginning of the chapter (divide both with J T)(f)e''^°^'^^ 
and use perturbation theory), however there are simplifications which should be 
taken into account. First we see again, that the 1-point function corresponding to 
the shifted action is zero: 



-0e^r[^.G] ^ ^rm (^^ffi^ - 0) = 0. (1.72) 



1.3. 2PI formalism 
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Figure 1.2: Self-energy insertions are resummed and a full propagator appears (thick 
line). This produces the absence of 2-particle-reducible diagrams. 

In the first equality we performed the change of variable —> — 0. As it was 
already argued, the 1-particle-reducible diagrams drop out. Furthermore, where the 
connected 2-point function (W2) appears, it can be substituted by the new variable 
G [7]. It can be seen by the following argument. Let us calculate the 2-point function 
of the shifted theory in two different ways. The first one is: 

j U^(p(pe V ^ ^ JJK(00+G)) ^ ^ JW\J,K] 

-00e^r[<^.G] ^ ^^mG]l ^''W[J,K] ^ gir[0,G]p^^_ ^.73) 
i oJ^ 

The other way is: 

-00e*^[^'^] = e^^t^'^] f 00 + G - = e^ri^'^'G. (1.74) 



The key part was to observe that the product of two fields could appear with 
differentiating with respect to K, or two times with respect to J. Comparing (1.73) 
with (1.74) we have W2 = G. Consider then a typical diagram of the numerator of 
(1.70) which contains two subdiagrams 7^^-' and 7*^2) connected by two perturbative 
propagators. Let us treat the subdiagram 7*^2) ^s a self-energy input. (We could have 
chosen 7*^^^ it does not matter.) Then find every possible diagram in which 7(1) is 
connected to an other self-energy insertion also with two lines. In other words, let 
us find all self-energy inputs, and connect them to 7*-^'' via two lines. In the sum 
of these diagrams, the connected 2-point function {W2) appears as it connects a 
part of 7^^'' with another one (see Figure 1.2). Since W2 can be replaced by G, we 
conclude that in r[0, G] diagrams with properties described above do not appear, 
when propagator lines are set to G. We say, that in r[0, G] only those diagrams are 
included, which do not fall into two parts if two internal lines are cut. These are 
to be called two-partide-irreducible (2PI) diagrams and they have to be calculated 
with vertices of the shifted action and with propagator lines equal to G. This is the 
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reason for calling r[0, G] the 2PI effective action, which can be calculated as: 



1.75) 



numerator: 2PI, W^—^^G 



In (1.75) there is no reference to the sources (since neither the 1- nor the 2-point 
functions are needed for the calculation), as it should. 

Careful readers may complain about the denominator of (1.75), which was used 
in IPI formalism to cancel diagrams of the numerator including only the original 
vertices of the theory (i.e. before the shift). Since now in the numerator these di- 
agrams were used to sum up the self-energy insertions, no cancellations occur at 
all. Therefore the denominator is just playing the role of normalization and is often 
neglected since it gives an additive constant for r[4>, G] and in practice we need only 
the derivatives of it. 

Prom the 2PI effective action we can obtain the standard IPI effective action. 
Using the definition, the IPI action can be thought as r[0, G]\k=o- In the language 
of the 2 PI action, the K — condition refers to a propagator, which satisfies 
the second stationary condition of (1.69) as a function of 0: 

G=G[^] 



SG 



0. 



(1.76) 



Therefore the IPI effective action is 



r U, G = G 



:i.77) 

Note that the equality is true only if the normalization is appropriately chosen (see 
above), otherwise a constant shift may appear. 

We can use the loop expansion of the IPl effective action determined in (1.47), 
since from (1.64) it is obvious that at leading order 

1 



F 



K\ 



S[<l>] + 




K{x, y)(l>{x)^{y) + - Tr log(l + {U (0) - K)i/^f), 



or equivalently 



r[0,G'] = 5[0] + -Trlog 



:i.78) 



(1.79) 



This is not the final form we are looking for, since the source K should be expressed 
with G and 0. This can be achieved using the second equation of (1.68). We find 
that at lowest order 



iG-^ ^i/^-p^ -U'\(t)) + K, 



(1.80) 
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which is consistent if we plug this relation to (1.79). With this we obtained the loop 
expansion of the 2PI effective action at leading order [6]: 



r[0,G] = 5[0] + ^I¥log ^ +lTr(D-iG-l), (1.81) 



2 



2 



where we introduced the short hand notation iD^^ = iAj} — U"{(f)), which is the 
classical (tree level) propagator at a given background in the absence of the bilocal 
source. If we want to go further in the expansion, we have to add 2PI diagrams con- 
taining at least 2 loops. To be very precise and taking into account the normalization, 
we have 

m G] = S[4>] + ^ Tr log + ^ Tr {D-'G - 1) + T^'^ - N, (1.82) 



2 



where iN — (every vacuum diagram with at least 2 loops built by perturbative 
propagators containing vertices of the original theory), and iF^^^ = (2PI vacuum 
diagrams with at least 2 loops built by propagators G with vertices of the shifted 
theory) . 

In the literature the and G independent parts are usually omitted, therefore 
the final form of the 2PI effective action is 

r2Pi[0, G] = S[^] + 1 Tr logG-i + ^ Tr (D-'G) + r^'\ (1.83) 

which is valid for every multicomponent scalar field theory. 

From stationary conditions (1.69) we get equivalent equations already obtained 
in the Dyson- Schwinger formalism. The equivalence is demonstrated as follows. The 
first observation is that the partial derivative of the 2PI action with respect to 4> is 
equal to the first derivative of the IPI action, if the propagator is set to G (i.e. the 
solution of the second stationary condition) : 



Sm _ Ar[4>,G[4>]] _ S^G] ^ ST[^,G] 



54) A(f) 5(t) 5G 



G=G 



where ^ refers to total differentiation. This shows that the condition of the physical 
field configuration being the extremal point of the IPI effective action is equivalent to 
the first stationary condition of (1.69) obtained in 2PI formalism, if the propagator is 
set to the solution of the second stationary condition of (1.69). The other important 
observation is that the latter equation in turn is found to be equivalent to the 
Dyson-equation. We have 
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which leads to 



iG {x,y) =T2{x,y) 



{d"^ + + U" {^))6{x - y) + 2 



(1.86) 



6G{x,y)' 



In (j)^ theory the diagrams of T^"^^ are 




+ ... 



(1.87) 



from which in the symmetric phase 




(1.88) is reminiscent of (1.61) and (1.36). At a first look it seems to be identical with 
(1.36), since a differentiation with respect to G removes a propagator line, therefore 
what remains is only 1-particle-irreducible, just as it can be seen on (1.36) (the fac- 
tor of 2 is also correct since two lines can be cut to obtain the very same diagram). 
However, at a second look we may think to have an apparent contradiction, since 
in r^^) we have diagrams which contain no self-energy insertions. This can lead to 
confusion, since in (1.36) every IPI diagrams are present, but in (1.88) not. The 
resolution is that we should keep in mind that in the right hand side of (1.88) the 
diagrams are made of full propagators (thick lines), and not of perturbative ones, as 
on the right hand side of (1.36). Indeed, we can organize the perturbative series on 
the right hand side of (1.36) in a way that the self energy insertions are resummed 
and full propagators appear. With this we arrive exactly at (1.88), which is called 
the self- consistent Dyson- equation. This equation appears in 2PI formalism in a very 
natural way. Since we have already argued that the Dyson-equation is identical with 
the Dyson-Schwinger equation of we obtained the same relation in three differ- 
ent forms. The difference between them is that while the ordinary Dyson-equation 
generates the 2-point function via perturbative propagators and vertices, the self- 
consistent Dyson equation and the Dyson-Schwinger formalism resum certain dia- 
grams. In Dyson-Schwinger equations diagrams are resummed into proper vertices, 
while in 2PI formalism (leading to the self-consistent Dyson-equation) perturbative 
vertices appear and only the self-energy insertions are resummed. 

Since T^"^^ contains infinitely many diagrams, analogously to the Dyson-Schwinger 
method, in practice we have to truncate the series. The truncation usually does not 
respect general identities of the n-point functions of the full theory. For instance we 
derive such an identity for the propagator (later we will mention a similar argument 
for the 4-point function as well). Since the 2-point proper vertex is T2 = ^g^^, using 
(1.77) and (1.86), we have the relation 




(1.89) 
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With the help of (1.84), the left hand side of (1.89) can be extracted: 



6^r 



2PI 



,G 



5'r2Pi[0,G] 



5r(2)[0,G] 



G=G 



5G 



;i.9o) 



G=G 



In the symmetric phase the second term of the left hand side of (1.90) vanishes, since 
in this case functions made of an odd number of field derivatives of r2Pi vanish (this 
is true in every theory with Z2 symmetry). It is convenient to define rint[0, G] = 
^2Pi[4>, G] — S[(l>] — I Tr logG"-*^, since when = (1.90) can be written as: 



which means that 



+ 



iD- 



<j)=0 



\<j>=0 



5G 



;i.91) 



0=0 



S^T,^t[^,G] 



(/>=0 



SG 



;i.92) 



<t>=o 



is an exact relation in the full theory. However, in general, the truncation of F^^^ does 
not respect it. It is possible then to define two propagators corresponding to the left- 
and right-hand side of (1.91). This plays an important role when one attempts to 
renormalize a specific 2PI approximation. (The left-hand side of (1.89) is sometimes 
referred as the curvature of the effective potential.) 

The truncation of F*^^) gives a partial resummation of the perturbative series. To 
get a taste how a specific resummation works, let us see the most simple example, 
the Hartree approximation. In this only the first diagram of (1.87) is kept. This leads 
to the following: 



whose iterative solution is 



id^ + m^)Six-y)-i _Q 



;i.93) 



F2(x, y) = —(9^ + m^)(5(x — y) — i 



+ ...).(1.94) 



(1.94) is the so-called super- daisy resummation. 



Chapter 2 

Effective scalar theories of QCD 



In this thesis we are investigating scalar field theories, particularly those which 
are effective models of the fundamental theory of the strong interaction, the quantiim- 
chromodynamics (QCD). The demand of building up such effective theories of QCD 
follows from the fact that it can not be solved at low energy scales due to a feature of 
the asymptotic freedom. This forces the running coupling of the theory to decrease 
monotonically with the energy, which is a nice effect at large enough scales, however 
around and below the scale of QCD (Aqcd) is rather inconvenient, since it breaks 
down perturbation theory. One way out is to give up this very successful method 
and apply lattice field theory, or we can build up effective models which are sensi- 
tive to specific features and based on the global symmetry properties of QCD. The 
spontaneously broken approximate chiral symmetry of QCD leads us to the idea of 
the linear sigma models, which we introduce in details in the following. 



At zero quark masses QCD with quark fiavors has a global Ul{N) x Ur{N) 
chiral symmetry, where L and R refers to the left and right projections of the Dirac 
quark field, respectively [1]. This means that a g = Qr + Ql quark field can be 
transformed without introducing changes into the Lagrangian as: 



where R and L refer to the left and right projections of the Dirac quark field, T" 
are the group generators and 9r and 6^ are transformation parameters. The basic 
properties of the U {N) group can be found in the Appendix. Alternatively we can 
write: 



2.1 The U(N)xU(N) model 



(2.1) 




(2.2) 
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where 9v,a '■— (^ii± ^l)/2 are the parameters of the vector and axialvector transfor- 
mations. 75 is the product of the Dirac matrices: 75 = i7o7i7273- This means that 
this global symmetry can be realized also in terms of vector and axialvector symme- 
try groups: Uv{N) x Ua{N), which can be decomposed into SUv{N) x SUa{N) x 
Uv{l) X Ua{1)- This is convenient, since the explicit, spontaneous and anomalous 
breaking of the symmetry is expressed through the vector and axialvector spectra. 
For instance, one can argue that due to the axial vector anomaly, in the quan- 
tum version of the theory the symmetry is reduced to SUv{N) x SUa{N) x f/y(l) 
[78]. Furthermore, nonzero, but equal quark masses reduce the symmetry further to 
SUv{N) X Uv{l), and with unequal masses only a Uv{^) symmetry remains. This 
refers to an arbitrary phase shift in the definition of the fields and is responsible for 
the baryon number conservation. 

In the U{N) x U{N) linear sigma model this symmetry is realized in terms of a 
new dynamical variable, an NxN matrix field M. A chiral transformation is defined 
as 

M — > RML\ (2.3) 

where R — exp (i^^T") and L = exp [iOIT""). For infinitesimal transformations we 
have 

SM = i{e''j{r''M - eiMT") = «(^^[t", m]+ + e^iT", m]_), (2.4) 

where [■,-]± refers to anticommutation and commutation, respectively. We have also 
introduced the vector and axial- vector transformation parameters as 9v,a — {Or ± 
Ol)/'^- From this, the pure vector and axial-vector transformations read as 

M — > VMV\ M — > A^MA\ (2.5) 

where V = exp (i6'^T'') and A = exp [id^T"). The U {N) x U{N) hnear sigma model 
is a chirally symmetric field theory of the M matrix. One can check that the most 
general renormalizable Lagrangian which preserves chiral symmetry can be written 
as 

£ = Tr {d^M^fd^M - rn^M^M) - gi[lh: {M^M)^ ^ - ^2 Tr [{M^Mf^ . (2.6) 

Since under pure vector and axial-vector transformations the determinant of M 
transforms as 



det (M) ^ det (yM0) = det (F Vm) = det (M), (2.7) 
det(M) det(A^MA^) = ( det (A^))^ det (M) = exp (2i^^ Tr T") det (M), 



2.1. The U(N)xU(N) model 
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the scalar combination det (M) + det {M"^) changes only with respect to the 

part of the transformation. This property is appropriate for describing the Ua{^) 

anomaly of QCD by introducing the following piece into the Lagrangian: 

£[^^(1) = c( det M + det M^) . (2.8) 

We note that this term gives birth to a non-renormalizable coupling for > 4. 
However, since above three quark flavors, the chiral symmetry of QCD is broken 
explicitly such that we do not expect that an effective theory based on this symmetry 
can be applicable at all. 

The simplest explicit symmetry breaking term is written as 

Ti{H{M + M^)), (2.9) 

which is built in a way to have the same transformation properties for the chiral 
transformation as the mass term of QCD. 

Introducing the scalar (s) and pseudoscalar (tt) fields as 

M = T''(s" + i7r"), (2.10) 
various pieces of the Lagrangian (2.6) can be written as follows [78]. 
jr = ^(aXa'^s" + aX5V-m'(s"s" + 7rV)) 

- \Faic<is''s\^s'' - ^F„;„,7r"7rV7r<^ - 2i/„b,cds"sV7r^ (2.11) 

where 

Fabcd -r i^ab^cd + ^ac^bd + ^ad^bc) ~l~ g {ddjjjyidjYiQd + dadmdmbc ~l~ ^acm'^m6<i))(2-12a) 

Hab,cd — -^^ab^cd ~l~ g {dabrndmcd facmfmbd fbcmfmad) ■ (2.12b) 

Alternatively, we can write 

(^s'^s^ + ttV)^- ^?7"^7", (2.13) 
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where 



= idaUs's'^ + TtV) - /„6esV. (2.14) 

The dabc totally symmetric and the fate totally antisymmetric tensors are defined in 
the Appendix. Furthermore, the term describing the anomaly is: 



^^(.(sV-TrV), A^ = 2 

G'^b,(s'^sV-37r«7rV), A^ = 3 (2.15) 

|G'afecd(6s"s''7rV - s^'s^s^s'^ - 7r«7rV7r'^), A^ = 4 
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where 



Gib = 2^SaoSbO - SalSbl - Sa23b2 - SasSbs), (2.16a) 

c / 3 9 \ 

^abc — g [dabc — g ('^oO^^Obc + <^60<^oOc + ^codabo) + 2 ^000<^aO<^60<^cO j , (2.16b) 

G abed -I 1^ ^ah^cd ^" ^ad^bc ^" ^ac^bc dabmdfncd dddfyidai^fjidfYibc dacmdmbc 

16 L 

+ 165ao5bQ5,i)5dO - 4:{6ao6boScd + ^aO^cO^td + SaoSdoSbc + SboScoSad 

+ ^bO^dO^ac + ^doScO^ab) + "/^i^aodbcd + ^bodcda + ^coddab + ^dodabc) .(2.16c) 

As it was already mentioned, the term containing the determinant only makes 
sense for N — 2,3,4 from the renormalizability point of view. Finally, with the 
parametrization H — hgT"', we have 

Ch = has\ (2.17) 

Although only low N values are physically relevant, we can try to attack the the- 
ory with the large- A?" technique and obtain quantitative results with continuing the 
results to these rather low values. In the end of this chapter, the reader can find a 
review of the results obtained in the past using large- A" technique. 



2.2 The 0(N) model 



Let us consider the U [N] x U [N] model with N — 2 and introduce the following 
multiplets: 

$ = (so,7ri, 772,773) = (s,7f), {'^0,^1,82,33) = {rj',ao). (2.18) 

With these notations, after some calculations, (2.11) extended with an explicit sym- 
metry breaking term can be written as 

L = ^(a^$^9''$-m|$^$)-^(<l>'^$)2 

+ l{d^.X^d''x-mlx'^x)-^{x'xf + Li{^:X) + haS% (2.19) 

where L/($, %) describes the interaction between the multiplets and 

m| = — c, = + c, X — 6gi + 3g2- (2.20) 

Since in the N = 2 case the strange quark mass is formally infinite, for consistency 
reasons 1]' and must also be treated as infinitely heavy particles, since for example 
from experiments we know that rj' is heavier than the strange kaon: ttik < rriri'. We 
note that it is the axial anomaly which pushes up the mass of the x multiplet. 



2.3. The large-N expansion 
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therefore for consistency reasons we should treat the combination + c as infinity, 
while — c remains finite. This means that x should be omitted from (2.19) and 
the SU{2) X SU{2) linear sigma matrix model transfigures into the 0(4) symmetric 
vector model: {ho = h) 

L = hd^^'^d^'^ - m|$^$) - :^($'-''$)^ + hs. (2.21) 

The fields of $ are usually called as sigma (s) and pion (tt). The 0{N) model is 
the generalization of (2.21), where instead of 3 pions, we have a number of (A^ — 1). 
This extension is motivated by the need of the application of the non-perturbative 
large- A?" expansion. In this case we factor out an N from the coupling: (m| = m^) 

L = _ _A_($r^^)2^^^^ ^2.22) 

since it can be argued that the coupling must scale with 1/A^ in order to have the 
free energy as an extensive quantity. In the large- A?" expansion we formally treat 
N ^ 1 and use 1/A^ as an expansion parameter. Surprisingly, the large- results 
evaluated at A^ = 4 provide phenomenologically quite interesting informations even 
at leading order. 

2.3 The large-N expansion 

Large- A?" expansion is a rather old classical non-perturbative tool of quantum field 
theory [16, 17, 18], in which instead of building up perturbative series in powers of 
the coupling, we treat 1/N as the small parameter. This leads to a certain partial 
resummation to all orders in the ordinary Feynman diagrams. 

Leading order solution of the large- A?^ expansion of the 0{N) model has been 
applied to interesting problems of finite temperature phenomena, in particular the 
study of the restoration of the SU{2) x SU{2) ~ 0(4) chiral symmetry of QCD 
[52, 53, 54]. One of its attractive points is that it preserves Goldstone's theorem at 
every order of the expansion, a feature not shared by all resummations of the original 
perturbative series. Being a resummation to all orders, it has some features which 
are absent at any given order of the perturbation theory, but which are believed to 
be true for the exact solution. On the other hand, related to the now well established 
triviality of scalar theories (sec e.g. [56] for a review), its spectra shows the presence 
of a tachyonic pole (Landau ghost) in the renormalized propagators. 

In early publications [17, 18, 19] the appearance of tachyons was considered 
as an inconsistency of the large-A^ approximation. In ncxt-to-lcading order (NLO) 
investigations of the expansion, started already in 1974 [19], the tachyonic problem 
seemed to be aggravated, because the renormalized effective potential appeared to be 
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complex for all values of the field. This led to the claim that the large- expansion 
breaks down, when one is going beyond leading-order. However, extensive studies of 
the 0{N) symmetric model [20, 57] established the consistency of the expansion 
for the effective potential and revealed its rich phase structure. Considered in a 
restricted sense, as a renormalized effective theory, the large-A^ expansion turned 
out to be a valuable nonperturbative tool, when applied to phenomena dominated 
by scales much lower than the cut-off. A strict cut-off version of the model was 
considered in [58] showing that with restrictions on the value of the background 
field the model has a phase with spontaneously broken symmetry, free of tachyons. 
A different solution to the tachyonic problem, called the tachyonic regularization, 
was proposed in [59, 60], where the tachyonic pole is minimally subtracted. 

Next-to-leading order results were also presented for the pressure of the rela- 
tivistic 0{N) model and applied to the physical (A^ = 4) sigma-pion gas [37, 38]. 
Some time ago the the method was initiated in the same model for the description 
of the chiral symmetry restoration [54, 89]. Further applications in modeling chiral 
symmetry breaking can be found even in technicolor theories [90, 91]. Applications 
to non-equilibrium phenomena, e.g. to the quantum dynamics of disoriented chiral 
condensates [105, 106, 107] also benefit from extending the number of flavors. Most 
of the studies were realized in a reformulation of the model in which the quartic self- 
coupling of the A^-component scalar field is replaced by an auxiliary field mediated 
interaction. 

An analogous development for the three-flavor linear sigma model model is hin- 
dered, because even the leading order solution of the large- A^ limit of the U (N) x 
U{N) symmetric model is unknown. To the author's best knowledge, no published 
attempts exist, which would go beyond the partial large- A?" treatments of the 0{2N'^)- 
symmetric part of the symmetry group (i.e. corresponding to the term proportional 
to s"s° -I- tt^tt"). The results obtained with such an approach are questionable in 
light of the rather different finite temperature renormalization group behavior of the 
0{2N^) and U{N) x U{N) symmetric models for A^ > 3 [77, 76]. In case of 0{2N^), 
one expects continuous (second order) transition, while in the U{N) x U{N) sym- 
metric model it is argued to be first order for A^ > 3 [77]. 



Chapter 3 



Renormalization of the 0(N) 
model in the large-N 
approximation 

In this chapter we investigate renormahzabihty of the 0{N) model in the large- 
hmit using 2PI formalism. Since most of the studies of the 0{N) model described 
in Chapter 2 were realized in a reformulation of the model in which the quartic self- 
coupling of the A^-component scalar field is replaced by an auxiliary field mediated 
interaction, we also start our investigation this way. 

We should keep in mind, that the renormahzation of the 0{N) model was per- 
formed at NLO in [19]. However, while pointing out all the divergences, this analysis 
did not attempt the explicit calculation of the counterterms, which were determined 
only at LO. Analyses of the renormalization of the self-coupling of the model showed 
that at NLO the /3-function is corrected by terms of order 1/N [61]. The success 
of the renormalization program in 2PI formalism [41, 43, 44] revived also the study 
of the renormalizability of the NLO corrections of the large- iV expansion [48, 62]. 
Despite all these efforts, the detailed and explicit knowledge of the counterterms 
was missing in the literature. Another, perhaps more serious motivation for the 
study is that recent publications even raised doubts on the renormalizability of the 
NLO approximation of the model formulated with an auxiliary (composite) field for 
arbitrary values of the field expectation values [37, 38]. 

In this chapter we provide an explicit construction of the renormalized effective 
action at NLO in the auxiliary field formulation of the 0{N) model. The leading 
order part of the effective action is proportional to the number of degrees of freedom 
A^, therefore our goal is to construct and renormalize it with 0{N^) accuracy. This 
goal will be achieved by separating divergences in the derivatives of the effective 
action and then integrating them to obtain the counterterms of the action and 
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eventually get a finite effective action. The identification of the divergences will take 
place by expanding the propagators in inverse powers of and omitting all terms 
which contribute beyond NLO. We note, that this step "kills" the self-consistent 
nature of the 2PI formalism, most of the propagators equal their tree level value at 
leading order, where by tree level we mean the expression which can be derived from 
the classical action including the auxiliary field. The construction of the counterterm 
functional will respect Goldstone's theorem, obviously obeyed in the NLO 1/N- 
expansion of the bare theory. 

The fact that the 2-point functions will be determined by following the hi- 
erarchy and not sclf-consistcntly has the definite backward effect of implying secular 
behavior in time dependent applications [64, 65]. It is however essential for ensuring 
Goldstone's theorem, since an 0{1/N) expansion on the level of the 2P1 generating 
functional was shown not to cure the violation of Goldstone's theorem by the self- 
consistent propagator [48]. Proposition for the preservation of Goldstone's theorem 
for the self-consistent (2PI) propagators exists at present only at LO (Hartrcc-level) 
[66] and the proposed solution follows efforts initiated in the framework of non- 
relativistic many-body theory [67] . 

By explicit construction of the counterterm functional in the auxiliary field for- 
mulation we demonstrate that the model is renormalizable at NLO in the large- A?" 
expansion for arbitrary vacuum expectation values of the fields. Substitution of the 
solution of the auxiliary field and the related propagators at the level of the 0{N^) 
accurate functional leads to the recovery of 0{N) and 0{N^) terms of the NLO 
2PI effective action of the 0{N) model in its original formulation [68], this time 
completed with all appropriate counterterms of the same accuracy. Alternatively, 
elimination of the NLO pion and LO sigma propagators produces the renormalized 
effective potential for the sigma and the auxiliary fields [37], again with correct 
counterterms. 

We remark that in the "complete" 2PI-1/N approximation, where one is looking 
for a self-consistent solution of the stationary conditions of the 2PI effective action 
with respect to the propagators (i.e. after building up the 2PI effective action, no 
further expansion takes place in 1/N) represents a different resummation strategy. 
For the sake of completeness this will also be presented in the last section of this 
chapter. 

3.1 Auxiliary field formulation 

At the level of the generating functional Z[J] (1.6), we introduce the auxil- 
iary field a into the Lagrangian (2.22) of the 0{N) model through the functional 
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Hubbard-Stratonovich transformation [70] : 
J Va exp Y J ^^"^ 



oc exp <i df'x 



X 



24:N 



{^\x)y 



(3.1) 



where $^ = + tt^. Then, the extended action without the exphcit symmetry 
breaking term and shifted in s by the background \/Nv reads: 



5'[s,7r„,a;u] 



d X 



+ 



m 



(^s^(x) + nl{x) + 2VNvs{x) + Nv'^^ 

a{x) (^s'^ix) + ttKx) + 2\fNvs{x) + Nv'^^ 



2 V 3iV 



(3.2) 



3.1.1 Approximate 2PI effective potential 



From the modified action (3.2), we construct an approximate 2PI effective action. 
Since the 4-point couphngs were changed to 3-point ones, the most simple truncation 
is to inchide only the 2-loop setting sun diagrams (check Fig. 3.1) containing 3 
propagators. As we have multicomponent fields and propagators, we have three of 
this type of diagram (sec (3.4) below). Because of translational invariance, the mean 
fields are spacetime independent and in this case we can always factor out a four- 
volume and work with the effective potential instead of the effective action: 



2PI 



r 



2PI- 



(3.3) 



Allowing the auxiliary field to have a saddle point value <a>= y using (1.83) 
Fourier representation of the propagators we have (the tt and the shifted s fields 
have no expectation values): 



V2Pi[a,v,GT,,g] 



( 



m — la 



+2A" -2 



{N-l){lnG-\k) 



+ D-\k)G^{k)) + Ti\Yig-\k) + Tr {V-\k)g{k)) 

+ '^11 Gaa{k)[{N-l)G,(p)G^(p + k) 

+ Gss{p)Gss{p + A;)) + 2GMGss{k)Gsaip + k) 

+ AV[a,v,G^,g]. 



(3.4) 
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Note that the momentum integrals are defined as = J In (3-4) m? and A 
are the renormahzed couphngs, and A\/[q;, f , G^^, Q] is the yet undetermined counter- 
term functional. Q and V are the variational and trcc-lcvcl symmetric 2x2 propaga- 
tor matrices with components Gss^ Gsa, Gaa and Dss, Dga, Daa, respectively. Their 
inverse matrices are denoted by and V^^. The elements of the inverse of the 
tree-level propagators appearing here are the following: 

t{D-Xa{k) = -1 , t{D-%,{k) = t^v, 

iiD-%,{k) = iD;\k) ^e- M^ (3.5) 

where we introduced the shorthand notation = — ia. 

It is important to mention that the 2PI approximation defined by (3.4) is identical 
to a specific truncation of the tower of the Dyson-Schwinger equations. We introduce 
the so-called Bare Vertex Approximation (BVA) as we close the Dyson-Schwinger 
hierarchy at the level of the 3-point function. In other words we reduce the tower 
to two class of equations (i.e. equations of the 1- and 2-point functions). The exact 
3-point vertices are denoted by TABc{x,y,z) = S^T / 6(f)A{x)6(j)B{y)S(j)c{z) (here (pi 
can refer to any of the fields TT,s,a) are approximated with their classical value: 



^' " Ssix)Liy)aiz) = 1 377^^" ' ^^^^^ " 

^7r„iTma{x,y, z) ~ T~r~r7~r~{ — ^nnJ^ ssa{x , y , z) . (3-6) 

dT:{x)d7:[y)a[z) 

The other 3-point vertices and all the higher ones are zero at classical level. With- 
out going into details, we just state that when we plug these expressions into the 
Dyson-Schwinger equations (actually they appear only in the equation of the prop- 
agator, check (1.57) for obtaining the derivatives of the effective action), we obtain 
the very same equations for the propagators and for the 1-point functions, as the 
ones which come from (3.4) after differentiating with respect to the propagators or 
the fields. This means that in this certain 2PI approximation it is possible to find 
the corresponding truncation of the tower of Dyson-Schwinger equations (or vice 
versa). Then, if we renormahzed the 2PI equations, we would immediately prove the 
renormalizability of the DS equations in the given truncation. 

We can easily figure out, that beyond the setting sun, no 2PI diagrams arc present 
which would contribute to the NLO effective action (see Fig. 3.1). In agreement, it 
reproduces the NLO 2PI effective action presented in (44) and (55) of [24]. Having in 
mind the 0{N^) (NLO) determination of the effective action we realize that only the 
setting sun involving pion propagators should be kept in (3.4) at NLO. Therefore we 
truncate further the functional in (3.4) and will work with the approximate effective 
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Figure 3.1: The first few 2PI diagrams contributing to V2P1. Tlie lines can represent 
various propagators allowed by the vertices. Only the first type (setting sun) can 
contribute at NLO of the large- iV expansion. 



{N-l){lnG^\k) 



potential 

1 37V^ % f 

V2Pi[a, V, G^, S] = - [m^ - ia) Nv^ + —a^ - - 

J k 

+ D-\k)GM) + TT\ng-\k) + Tr {V-\k)g{k)) 

+ ^^ JJ^GUk)G^{p)GAP + + AV[a, v, G^, G]. (3.7) 

This truncated form of the functional influences the propagator equations of the 
a — s sector (at LO) and of the pion fields (at NLO). We emphasize, that even if we 
preserved the omitted terms, the 0{1/N) solution of that approximation would not 
correspond to the complete NLO order solution in the a — s sector, since there are 
more 2PI diagrams proportional to 1/N. 

As already announced, renormalizability of this approximation will be investi- 
gated at the level of the derivatives of (3.7), then we attempt the construction of 
an appropriate counterterm functional AV[d, v, Gt^, Q]. This counterterm functional 
allows for a uniform treatment of the counterterms and also makes possible to keep 
track of the effect a counterterm determined from one equation has in the renormal- 
ization of another equation. Some new insights will be offered when the results are 
compared to other approaches, where the propagators are not considered as vari- 
ational variables [37, 38], or when the quantities related to the auxiliary field are 
eliminated [68]. We shall demonstrate the validity of Goldstone's theorem also for the 
renormalized NLO pion propagators. Before we start to determine the counterterm 
functional, we have to fix the renormalization scheme we are using. 



3.1.2 Renormalization scheme 

Following the spirit of Ref. [50] we will expand the propagators around appropri- 
ately chosen infrared safe auxiliary propagators (which has nothing to do with the 
auxiliary field!) and subtract the divergent pieces defined through these quantities 
(it reminds us to the minimal subtraction scheme). These subtractions define our 
renormalization scheme. 
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An important fact for finding appropriate auxiliary propagators is that we see 
from (3.4) that only the a — a component of the Q propagator receives a correction 
at leading order: 

i{G-Xo.{k) D^{k)D^{k +p)\^ = -l + ^I^. (3.8) 

Note that we took the finite part of the /^r integral immediately, since an appropriate 
counterterm of AV should kill the divergent piece. We will come to its construction 
later. The corresponding inverse is therefore (the index corresponds to the leading 
order expression): 

Gi°;(^) = (l - ^I^ik)^ G{k) , G<^lik) = -iD;\k)G{k), 

G^^Kk) = ^^vG{k), (3.9) 



where 



iG-\k) = {e -M')[l- ^I^ik)] - ^v'. (3.10) 



It will turn out that at NLO the asymptotics is determined basically by integrals of 
the LO propagator Dt^{p) and of the propagator G{j)) given in (3.10) which incorpo- 
rates the effect of the resummation of the chain of pion bubbles (Hint: expand the 
inverse of (3.10) in terms of A!), we need two auxiliary propagators. The first one is 

where Mq is an arbitrary mass scale. With this propagator we define the quadratically 
divergent integral 

T^f = / Go{pl (3.12) 
and the following one-loop bubble integral 

hip) = -^ ^ Go{k)G^{k + p) = + I^{p), (3.13) 
where the logarithmically divergent part of the integral above is defined as 

Tf^ --if Glik), (3.14) 



and Iq is the finite part of the bubble integral. The explicit cutoff dependence of 
the introduced divergent quantities are: 
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which were calculated in Euclidean space with po ip^ and a 4d cutoff A. The 

^0^ 



finite part of the bubble integral behaves asymptotically as /q (p) ~ In — 2 — i7r + 
O (p~^lnp), and defines together with Gq{p) the second auxiliary propagator: 



Integrals involving combinations of Go(p)) Ga{p) and /(f (p) will be fully included in 
the counterterms. 



3.1.3 Leading order (LO) construction of the counterterms 

We start our program with a short description of the leading order (pieces pro- 
portional to N) analysis. Since the pion propagator equals its tree level value at 
LO, first we will discuss the saddle point equation and obtain the corresponding 
LO counterterm AF"^'^. Then we turn to the equation of state and the propagator 
equations in the a — s sector. At this level of the approximation, corresponding to 
this sector only one piece of AV will be necessary, which is denoted by AT/"" and 
is proportional to the a — a propagator. 

Saddle point equation (SPE) 

Taking the derivative of the functional in (3.7) with respect to a we arrive at the 
expression: 



5a 



37V A/" / f \ 
[a, V, G,, g]^—a- i- U + D^{k)\ + (3.17) 



where we replaced Gtt, originally appearing in the integral above by introduced 
in (3.5). The last term is the contribution of the counterterm functional which has to 
be constructed to ensure the finiteness of this equation. As announced, the structure 
of the tadpole integral is given in the third section of this chapter. Since we have: 

/D.(fc) = Tf + (M^-Mo^)Tr+T;, 

Jk 

we can see that the finiteness of (3.17) (T^ is the finite part of the pion tadpole 
integral) for any value of a is ensured by the following 0{N) piece of the counterterm 
functional: 

/V" r / 1 \ 

(3.19) 



2 



2 
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Here, and are the quadratic and logarithmic divergences of the pion tadpole 
given by the propagator Go{p), see (3.12) and (3.14). Then, from (3.17) we obtain 
the finite saddle point equation: 

-Ya-i-{v' + T:)^0. (3.20) 



Equation of state and Goldstone's theorem 

The leading order term in the derivative of the functional V2P1 with respect to 
V is of order N 

SV2P1 , 



5v 



-[a,v,G^,g]^ NvM\ (3.21) 



The right hand side is finite in itself, it does not necessitate the introduction of any 
0{N) counterterm. 

Since the equation of state is obtained by equating the r.h.s. of (3.21) to zero 
and the leading order inverse pion propagator is of the form (3.5), we immediately 
see that Goldstone's theorem {D~^ {k = 0) = 0) is obeyed. 

Leading order propagator matrix in the a — s-sector 

As it was already stressed, the only entry of the 2x2 inverse propagator matrix 
which receives LO correction is {G~^)aa '■ 

A AAI/"" 
i{G-\a{k) = -1 + -h{k) - 2-—, (3.22) 

b dGaa 

where /jr is defined in (3.8). Its finite part can be obtained in terms of as fol- 
lows. Recall that the divergence of the one-loop bubble integral built by propagators 
Go is chosen to be T^°^ given in (3.14). Power counting shows that changing the mass 
of the propagators in this integral leads only to a finite difference. In other words, 
the divergence of a perturbative bubble integral does not depend on the masses of 
the propagators. Then we define as: 

Up) --if D^{k)D^{k+p) ^ Tf + ip), 



167r2 



(3.23) 
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Here, the momentum dependent function M) is 



2^1-^ arth (^^1 _ M|l j _ y 1 _ 4Ml^ p2 > 

- 1 arctan > > q, (3.24) 

2^1 - ^ arth (Jl - ^)"' p2<0. 



It can be shown that hm L(p^, M) — 2, which results that 

^(f-0)^^.„|. (3.25) 

Since /(f (p) has exactly the same form as I^{p), but with replaced by Mq, for 
large we have 

We see now that in order to make (3.22) finite we have to introduce another piece 
into the counterterm functional: 

/W'^'^^^Tf^ jcUk). (3.27) 

With this the LO renormalization is complete. Note, that in the broken symmetry 
phase all elements of the LO propagator matrix have common pole structure de- 
termined by G{k). This is the manifestation of the hybridization for v 7^ of the 
longitudinal field component s and of the composite field a ~ + tt^ [71] . This 
feature is relevant when studying dynamical aspects of the phase transition at finite 
temperature. 

Using the first entry of (3.9), we can derive the following relation between the 
LO sigma and pion propagators: 

This relation will be very useful for the divergence analysis at NLO. 



3.1.4 Next-to- leading order (NLO) construction of the coun- 
terterms 

The construction of the NLO counterterm starts with discussing the expan- 
sion of the pion propagator. Its asymptotics will be completely determined by certain 
integrals of the LO propagators. We shall see that the counterterm functional can be 
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determined by the asymptotics of the LO propagators and of the NLO self-energy 
of the pions. In order to exphcitly demonstrate the NLO renormahzabihty for all 
values of v and a we investigate the derivatives of the functional (3.7) with respect 
to these variables. 

First we analyze the divergence structure of the pion propagator equation, and 
determine the corresponding piece of the counterterm functional, denoted by AV^. 
Then we turn again to the equation of state (now at NLO). To have a finite version, 
we will determine the next piece: AV^. At this point we have the unrcnormalized 
NLO saddle point equation of a, and propagator equations in the a — s sector. It 
will turn out that both AV" and AV^ depend on a, therefore they give contribution 
also to the SPE. The mutual consistency of the countcrtcrms renormalizing the SPE, 
EoS and the pion propagator equation is fundamental for the consistent outcome 
of our analysis. During the renormalization of the SPE, this consistency is checked 
and also a new counterterm piece AV^'^ will be introduced, which completes the 
renormalization. 

We do not present here the renormalization of the partial NLO corrections which 
would occur in the a — s sector when we used the 2 PI effective action (3.4). All 
counterterms necessary for the renormahzation of these pieces are of 0{1/N), not 
contributing to the renormalized effective action at NLO (i.e. 0{N^)) in the 1/N 
expansion. 



Pion propagator 



The pion propagator at NLO in the expansion is given by 

iG-\k) ^ iD-\k) - I— jG^l{p)D^ip + k)-2—. (3.29) 

Since we need the pion self-energy to 0{1/N) accuracy, we replaced G^a with Gal 
and with D^^ in the above integral. 

In order to determine the counterterm contribution in (3.29), we have to study 
the divergence of the NLO self-energy. Using the first two entries of (3.9) together 
with (3.28) we can write 

= I ^ (s SO) 

""^^^ l-A/.^(p)/6 3 (l-A7,^(p)/6)^- ^ ' 

First we show that the divergence of (3.29) is momentum independent, i.e. there 
is no need for infinite wave function renormalization, in accordance with [59]. It is 
clear by dimensional reasons that we have to analyze that piece of the second term 
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of (3.29), which is proportional to 

Dnjp + k) 
,,l-XI^{p)/6- 



I 



(3.31) 



To study this one, we use the following expansion 



In order to prove that there is no infinite wave function renormalization, it is enough 
to look at the appearance of k'^ in the numerator, that is at terms n — 1, 2 in the sum. 
Keeping only terms up to and including 0{l/p^), but throwing away those which 
vanish upon symmetrical integration (note that I^ip) depends on p^) amounts to 
the following replacement at the level of the integrand in (3.31): 

1 . 1 , 4(p ■ kf - fcV 

{p + ky - M2 p2 _ 2 + (p2 _ ^2)3 ■ y ■ ) 

The second term gives vanishing contribution in (3.31) due to the following property 
which holds for any integrable function f{p^) upon the use of a Lorentz-invariant 
regulator: 

J d'pp,pj{p') = ^J d'pp'fi/). (3.34) 

It is clear that the other piece of the second term of (3.29), which is coming from 
the second term of (3.30) and proportional to does not have a divergent part with 
respect to the external momenta. In fact it does not contribute to the divergences at 
all, since after setting k — Q and iterating (3.28) we recognize that, by logarithmic(!) 
power counting the following integral is actually convergent: 

DI{P) 



I 



P (1 - XI^{p)/Q) 



2 • 



(3.35) 



Due to this fact we do not encounter any divergence proportional to v'^ in the NLO 
pion self-energy. 

It is instructive to point out here a peculiarity of the rcsummation procedure as 
compared to the order-by-order perturbative renormalization. Namely, when in the 
above integral the denominator is expanded in powers of A then at nth order of the 
expansion, we find a A"(logA)"+^ type divergence. Through formal resummation of 
this divergent scries a finite result is obtained. This argument explicitly shows that 
in resummcd perturbation theory the structure of the counterterms can be different 
from that seen at any given order of the perturbation theory. The same effect was 
noticed in [59] in connection with the wave function renormalization constants of 
pion and sigma fields which arise from imposing renormalization conditions on the 
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residua of their propagators. At NLO in the expansion they are finite whereas 
in an expansion to any given order in the couphng they appear divergent. 

Wc see now that the momentum independent divergence is determined by the 
first term of (3.30): 



div 



, 1 - XI^{p)/Q 



(3.36) 



div 



To find the expression of Tdiv(M^) using the auxihary propagators we start by taking 
into account that the one-loop bubble integral behaves logarithmically for large 
momentum. In the asymptotic momentum region the above expression allows us to 
write 



1 - XI^ip)/6 1 - A7o^(p)/6 ^ 6 (1 - A7o^(p)/6)2 

1 - A7o^(p)/3 



+ 



1 - A7o^(p)/6 p2 _ (1 - A7o^(p)/6)2 



Inp 
+ 



p2 2p 



,(3.37) 



where in the last line we used 7(f (p) ~ In — 2 — iTi + O {p ^ Inp) . The neglected 
terms give finite contribution in the integral of (3.36). Using there (3.37) and 



p2_M^ p2_M^ ip^-M^) 



(3.38) 



we obtain 



^div(M^) = Tj^^-{M'-Mly jGl{p) [l - ^I^{p) 
where the following divergent integrals were defined: 

T^f = / GM. TP ^-i r Gl{p)I^{p). 

J P dp 



-^{M^-Ml)Ti'\ (3.39) 



(3.40) 



Using in the remaining integral of (3.39) that JpGl{p) is finite, we obtain for T^iy 
the expression 



fdiv(M2) = - ^(M^ - M2)ry). 



(3.41) 



Since the necessary counterterm in (3.29) compensates Tdiv; we readily find the 
required counterterm functional piece upon functional integration of (3.41) with 
respect to Cj^ and multiplying it by A/6: 



A 



(3.42) 
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We conclude this part by giving the finite pion propagator at NLO in the 
expansion, including also the contribution of the counter-functional AV". It reads 
as 

iG-\k) ^k'-M'- ^^:{k), E^ik) = t J G^^l{p)D^{k+p) - f,^{M% 

(3.43) 

It is instructive to evaluate explicitly the cutoff dependence of the divergent integrals 
Ti^'' and Ta \ This can be also of some practical interest when we proceed to the 
numerical solution of the renormalized equations. As we did during the calculation 
of Tj"'' and T^^"*, we go to Euclidean space with po W% and use a 4d cutoff A, 
however we limit ourselves to an asymptotic analysis and expand the integrands for 
large momentum (check the asymptotics of the second entry of (3.23) for M — Mq). 
Exploiting the freedom to omit those contributions to T^^^ which are formally finite 
for A — > oo, we choose the scheme in which 

^ r d± ln(A:VMo^) 

{At^Y J A; (l + 2a-aln(A;VM2))2- > 

Here, we introduced a — A/(967r^). We notice that for k — Moexp(l -|-487r^/A) the 
denominator of the integral above vanishes. To avoid this non-integrable singularity 
to occur in the range of integration, that is for < A, we need A < Amax ~ 
Mocxp (l + 487r^/A). That means that for A 7^ the cutoff cannot be sent to infinity, 
there is a maximal value for it, which reflects the triviality of the theory. Performing 
the integral in (3.44), and obeying this restriction, we find 

With the same strategy we have 



Stt^A 



e^+-^AEifln^-2-^ 

A 

, , A A2 A 
+ 31n -;-r^ln— . + 1 + 



(3.46) 



967r2'"M2 ' ^ ' 487r2 
where Ei(a;) refers to the exponential integral function: 

Ei(x) = V j -dt. (3.47) 

An asymptotic series representation of this function looks as 

^ , , exp (x) n! / . x 

Ei(^) = ^^E^' (3-48) 

n=0 
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from which it is obvious that for large cutoffs Ti ~ A^/logA. In the hght of the 
previous two equations, it is important to mention that a similar consideration ap- 
plies to the divergence appearing in (3.36) to what was already stressed below (3.35): 
resummed perturbation theory can change the nature of the counterterms compared 
to order-by-order perturbation theory. Expanding the denominator of the middle 
expression of (3.36) into powers of A, we find the usual quadratic and logarithmic 
divergences characteristic for the perturbative contributions. Resummation however 
modifies these divergences as can be explicitly seen in the formulas (3.45) and (3.46). 

Equation of state and Goldstone's theorem revisited 

Next, we investigate the renormalization of the derivative of the 2PI effective poten- 
tial with respect to the background v. At NLO in the 1/A^ expansion this is given 



where for the second equality we used the last entry of (3.9). The counterterm 
functional AV^ determined above does not contribute (neither AV""^ nor AV^'-^ 
do), since its derivative with respect to v is zero. 

The equation of state is obtained by equating the r.h.s. of (3.49) to zero. Its 
unrenormalized expression obviously implies when confronted with (3.29) the valid- 
ity of Goldstone's theorem with 0{1/N) accuracy. Note that, as it is well known, 
Goldstone's theorem is not followed if we proceed in strict 2PI sense which requires 
the self-consistent determination of the full propagator without expansion in 

In a renormalization procedure which preserves Goldstone's theorem we need to 
construct a counterterm which does not depend on Gtt, therefore does not interfere 
with its already renormalized equation. Since the divergence in (3.49) coincides with 
the divergence of the NLO pion self-energy, the new contribution to the counterterm 
functional is obtained upon integrating with respect to v the expression given in 
(3.41) multiplied by Xv/3. We find 



by 




(3.49) 



AV 



A 



(3.50) 



6 



A remarkable feature of this 0{1/N) renormalized solution is that it satisfies Gold- 
stone's theorem for arbitrary values of a ! 
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Saddle point equation 

In writing down the derivative of the effective potential with respect to a, we must 
not forget about the contributions of the d-dependent counterterms AV^"'-'^, AV^, 
and AV constructed above (see (3.19), (3.42), and (3.50)): 



SVn 



2PI 



Sol 



+ I 



3N 
X 



N 



Tf + (M^ - M^,)T^ 



G^{k) - 



GUk) - G^{k) 



(0) 



+ 



2 

^Al^"'°(a) 
5a. 



A2 



(3.51) 



The contribution of the counterterms determined by the renormalization of the equa- 
tion of the inverse pion propagator (3.29) and of the equation of state is displayed 
in the last but one term. 

The last, yet undetermined part of the NLO counterterm functional, i.e. Ay"'°, 
provides the NLO completion to AF"'^. In order to determine it, first we have 
to evaluate the pion tadpole in the second term of the r.h.s. of (3.51) with NLO 
accuracy. Taking the inverse of G~^ given in (3.43) and expanding it to C(l/A'^), we 
obtain 



f Dl(k) jGl;>(p)D,{k +p) + ^if,UM^) ^ Dl(k) 



J{M') - ffJM^) / Dl(k) 



2^,2 



9N 



■J{M^), (3.52) 



where for the second equality we used (3.30) and introduced the following functions 



J(M2) 



Dlik) 



Dlik) 



f D^{p + k) 
D,{p + k)Gf}{p) 



(3.53a) 
(3.53b) 



„ J, (1 - \I^{p)l&Y 

Collecting all C(iV°) (NLO) divergent terms in (3.51) we see that 5AV'^'^{a)/5a is 
determined by 

5AT/"'0(a;) A 



5a 



+ 



Jdiv(M2) - fdiv(M2) / Dl{k) 

Jk 



i^^Vdiv(M2) 



j^{G^S{k)-D^{k)) 



div 



where Jdiv(M^) and Jdiv(M^) denote the divergences of the integrals defined in 
equations (3.53). Note that to the order of interest it was again allowed to replace 
in the last two terms the original C^r and Gss by -Dtt and Gf}, respectively. 
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The important question of consistency inquires whether the previously con- 
structed countertcrms cancel all sub divergences of J(M^) and the f^-dependent 
divergence of second and third terms in (3.54). Note that the double integral J(M^) 
has only overall divergence, since both k and p integrals are individually finite. 

We start with investigating the second double integral given in (3.53). Changing 
the order of integration we use (3.28) and the following relation which can be derived 
from (3.23): 



Dl{k)D^{p+k) = 



1 d 



to find 
Jdiv(M2) 



2dM2 



4M2 



167r2 



, M2 
In ^T7T + 



1 



Ml 



87r2 



, (3.55) 



;(l-A/^(p)/6)2p2_4M2 



div 



[Gl{p)Io^{p)^n^\ (3.56) 
Jp 



To obtain the second equality above we replaced with 1^ in view of (3.26) and 
used (3.38). 

We continue with the differences of the tadpoles involving LO sigma and pion 
propagators. Using (3.28) iteratively we find 



(G(°)(A;)-D.(fc)) 



A 



div 



Dm 



6 Al-A4^(fc)/6 



(3.57) 



div 



In view of (3.37), we can replace I^{k) with /(f (A;) in the denominator above. Then 
using (3.38) we obtain 



{GfXk)-Dm) 



div 



A 

— -i 
6 



G'm [ 1 - li!{k) 



div 



A2 
^36 



v^T^'\ (3.58) 



With this term the I'^-dependent divergences are canceled by the counterterm con- 
tribution proportional to appearing in the last term of (3.54). 

Let us now go through the detailed analysis of the divergence structure of J(M^), 
which is the first double integral given in (3.53). This one contains an overall di- 
vergence as well as subdivergences. By changing the order of integration and using 
(3.55), we have 



Jdiv(M2) 



4M2 



div 




4M2 1 - A/^(p)/6 



, (3.59) 



div 



where we have separated a divergence independent of I^{p), which by using a relation 
similar to (3.38) can be expressed as a hnear combination of t]^-* and T^^\ Since 



3.1. Auxiliary field formulation 



49 




Figure 3.2: The appearance of vertex type subdivergences in the first double integral 
of (3.53) as illustrated at leading A order in the expansion of the denominator of the 
integrand. The factor of three indicates the possible ways in which two lines of the 
setting-sun diagram form the bubble which corrects the vertex at 2-loop level. The 
cross denotes the associated lowest order counterterm. 

the form of the second integral in (3.59) is similar to the last one in (3.36), the 
calculation follows very closely the determination of T^iy. Using (3.41) and (3.25) 
the result is 



A ~ 
6 



Jdiv(M^) = irf + (4M2 - M2)irf 



+ ^ (^M^Tf ) - f,,(M^)) (l + ^Tf + ^ - ^4(p = 0)) .(3.60) 



The presence of subdivergences is reflected by the appearance of divergent terms 
proportional to ln(M^/MQ) (coming from /^(p = 0)). These should cancel if the 
approximation is renormalizable, since in this case only divergences proportional to 
powers of a (that is M^) are allowed. The result given in (3.60) for Jdiv(^^) shows 
the cancellation of the subdivergence Tdiv{M'^)lT,{p = 0), when it is combined with 
the second term of the square bracket of (3.54). The cancellation of this self-energy 
type subdivergence of the double integral J(M^) is expected in view of (3.43). 

The double integral Jdiv(M^) has also a vertex type subdivergence as illustrated 
in Fig. 3.2 at leading order of the expansion in A. This is cancelled by the last integral 
of (3.54). We can see this analytically by separating in the difference of the two terms 
of the square bracket of (3.54) a contribution proportional to T^^'T^ which on its 
turn cancels the contribution of the last tadpole integral in (3.54). 

With all the subdivergences and I'^-dependent divergences of (3.54) cancelled, 
SAV°''^{a)/6a is determined by the overall divergence of J{M'^) and that of the last 
tadpole integral. Its expression reads: 



5a 



12 167r2 

(3.61) 
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Since the above expression depends only on dt it will have no "back-reaction" neither 
on the propagator equations nor on the derivative of the effective potential with 
respect to the background. 

In conclusion, only divergences proportional to zeroth or first powers of 
remained which upon integration over a. determine the a-dependent counterterm 
functional Ay°'°(a). The counterterms induced by the renormalization of the NLO 
pion propagator and of bV^Pi I ^"^ played an essential role in the cancellation of sub- 
divergences and of f ^-dependent divergences present in the expression of 5V2Pi/5a. 
No limitations whatsoever showed up on the value of and/or f^, in contrast 
to the findings in [37, 38]. We shall return to the discussion of this difference after 
reducing the 2P1 effective potential to the ordinary IPl one depending only on the 
background fields. 



The explicit form of the counterterm functional 



Now we collect into a unique expression Ay all individual pieces determined in 
Eqs. (3.19), (3.27), (3.42), (3.50), (3.61), and express it in a conventional form, in 
which we associate them with the renormalization of different couplings appearing 
in the terms of V2Pi[6i, v, G^^, Q] (Eq. (3.7)). The counterterm functional reads: 

AV[a,v,GT^,Q] = - (^Sm^ — idgaj + iSKia + 6K,2Ci'^ 

+ \{?>rr? -ibga) [ G^{k) + lsKo I Gaa{k), (3.62) 
^ Jk ^ Jk 

where the counter-couplings are given by the following expressions: 



5k2 
5ki 



X 
3 



7V + 8 (0)^ A 

-d 



A. 



1 + ^rf + 



A 



647r^ 



6 » ' 
5Kn 



A. 



(0) 



Tf^ + (m^ - M^)rr + + (4m^ - M'^)T^ 



(2) 



.(0) 



-TP { 1 + ^Ti°) + 



487r2 



A 



.(0) 



+ ATi^) 



(2) 



+ ^ ( 2rr(m^ - M^) + rr + :^(3m^ - 



167r2' 



(3.63) 



It is interesting to note that the term proportional to 6k,i in (3.62) has no corre- 
spondent in (3.7). Moreover, the terms proportional to Skq and 6k,2 correspond in 
(3.7) to terms not proportional to any coupling of the original theory. However, the 
symmetry of the auxiliary field formulated model and the general concept of renor- 
malization allow the appearance of these terms. Skq and 5k2 correspond to different 
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possibilities of defining tfie two-point function of the auxiliary field (see (1.91)). They 
are calculated to different orders in since the corresponding terms contribute 
at different levels of the expansion, but after scaling back the fluctuating part of the 
field a. to a, the two counter couplings do agree to leading order in 1 /N . This feature 
is reminiscent of the renormalization of the operators corresponding to two different 
definitions of the n-point functions in 2Pl-approximations. The terms with 5g can be 
regarded as renormalizing the coupling g = y^XJJSN) through which the auxiliary 
field couples to fields of the 0{N) multiplet. We interpret the renormalization of all 
these operators not occurring in the original model as renormalizations of the two 
variants of the 2-point functions of the composite field 0^ and of the (fp' — s — s vertex 
(see ch. 30 of [2]). 

It is notable that 5k2 determines the renormalization of the coupling A following 
the formula 



6 6 4 6 iV + 8 (0)^ A A (0)^ A A 



This can be seen by looking at the terms proportional to 6? in the classical part of 
the functional. This relation is only intermediary, since it will change in the course 
of the ehmination of a. The part of the counterterm 5k2 proportional solely to Tf' 
appears already in the literature and follows the one-loop ^-function of the model 
[72]. The entire NLO part of 5k2 proportional to T^^^ is missing in the analysis of 
[37, 38] (see e.g. Eq. (23) of [38] for the expression of their counterterm hi). 
Introducing the following notations 

5n2^N5nf + 5n^2\ AB = A + 5Aa, a« = + ^a^, (3.65) 



we readily obtain 



= , A^) = A + 



6l + ATf/6' " l + ATf/6' 

5Ai^^ = -\{xff5K^\ (3.66) 

Comparing with Eq. (25) of [51] we observe that dXa'' is the counter-coupling of 
the 0{N) model formulated with its original variables and considered at LO in the 
large-A^ expansion. Likewise A^^ is the bare coupling of the model in the large-A^ 
limit. Wc shall sec in the next section, where the auxiliary field will be eliminated, 
that at NLO in the large- iV expansion the bare coupling of the 0{N) model differs 
from As, it turns out to be a combination of A^ and 5g (see (3.76) and (3.77)). 
With the counterterm functional explicitly determined by (3.62), we can 
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give now in a compact form the functional introduced in (3.7): 



[ [Tt In g-\k) + {N - 1) In G-\k)] 
Jk 

[k^ -ml + ica] [{N - l)G^{k) + Gss{k)] 
- iv^j^Gasik) + ^J^GUk) - ^ j^GUmik), (3.67) 



I 

2 
1 

2 



where Xb is defined in (3.64) and we introduced the following notations: 

1 1 
m| = + j:^^"^^-, c = 1 + ^'^5', c = 1 + 8kq, 

Ii{k) = -i I G^{p)G^{k+p). (3.68) 
Jp 

3.2 Elimination of the auxiliary field 

In this section the 0{N^) accurate rcnormalized functional will be established 
for the original formulation of the 0{N) model by eliminating the auxiliary field a 
and the propagator components related to it (i.e. Gas ^-nd Gaa)- In order to achieve 
this, we substitute the expressions of a, Gas, and Gal as found from their respective 
equations into (3.67). 

3.2.1 Determination of V2Pi[t', Gt^, Gs] 

For rewriting the terms depending on Gal and Gal we exploit their representa- 
tions which allow the expression of the result fully in terms of n(/c) and gI?. The 
latter will be replaced with Gs, the exact sigma propagator of the 0{N) model. In 
this way we find 

--jTr\ng-\k) = -^lndet^-^(fc) = -^ln [(^l-^/,^(fc))zG:^(^) 

= 4i^^^^^(^)-U^^^-^n(^))-UM^c) 

+ (3.69) 



In going from the first to the second hne above we used that to 0{1/N) accuracy 
-i I G^{p)G,{k+p) =-i [ D,{p)D,{k + p) ^Tf\ (3.70) 

Jp div Jp div 
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and therefore we have 

1 - ^I^ik) = c - ^L{k) = c - ^n(fc) + (^) , (3.71) 

where the neglected 0{1/N) contribution is finite. 

Using (3.30) for Gal we write the following chain of expressions for the last term 
of (3.67): 




In writing the second line above we use again (3.71). The first term on the r.h.s. 
above is cancelled against the last but one term of (3.67). Finally, for the last term 
in the second line of (3.67) we use (3.71) to write: 




(3.73) 



As a short digression from our current task we mention that we could proceed by 
further eliminating also the pion and sigma propagators using their respective NLO 
and LO equations, then we would obtain the renormalized version of the ordinary 
(IPI) effective potential as function of a and v studied in [37, 38]. We will come to 
that point later. 

Now we proceed instead with the elimination of a from (3.67) keeping the vari- 
ables of the original theory, i.e. v, G,r, and Gs- The simplest way to do this is to 
complete to a full square the functional depending quadratically on a and then use 
the saddle point equation 6V2Pi/6ct = 0. Replacing ci? by the exact propagator Gg 
of the 0{N) model to obtain from the a-dependent part 

2 

(3.74) 

Putting together all above pieces we also use that in view of (3.66) A/c = Ag^ and 



2m 



Ncv^ + c [ [{N- l)G^{k) + Gs{k)] - 25ki 

Jk 



54 Chapter 3. Renormalization of the 0(N) model in the large-N approximation 



obtain 



V2Pl[v, Gtt, Gs 



N 



2 XbcSki\ ^2 ^ jy^RC^ 



mo - 



3N 
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V 



24 



[(iV-l)lnG;^(^)+lnG7^(^)] 



si 



X(0) 



- ml + 



3N 



6 



-V 



Gs{k) 



6 Al-A^"^n(A;)/6 2 



[{N-l)G^{k) + Gsik)] 
U{k) 



\(0) 

In ( 1 - ^ 
6 



24 



G^k) 



12 



^.(A;) / Gs{p), (3.75) 



where we omitted terms of order 0{1/N) and a divergent constant ~ coming 
from (3.74) as well as the irrelevant last but one terms of (3.69) and (3.72). 

The bare couplings appearing in (3.75) have to be given with an accuracy cor- 
responding to that of the functional (3.7), therefore we write the counterterms as a 
sum of the LO and NLO contributions: 



Xbc^ = A + (5A, 



2 Xbc5ki 2 , e 2 

zrzz — = m + dm , 



3A^ 



(3.76) 



With help of (3.64), (3.65), (3.66), (3.68) and the separation 5ki = N5n>^^+5k\ \ 
the counter-couplings above can be given in terms of the counter-couplings of the 
model in the auxiliary field formalism as 



5m'^'^ = 



OKi , 



Sm 



2(1) 



Sm'^ 



1 



SX^^Sk^^ + Ag) {6k? + . (3.77) 



Using (3.76) in (3.75) we obtain at NLO in the 1/N expansion the renormalized 
functional of the original 0{N) model, that is without the auxiliary field, in the 
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following form: 



V2pi[v,Gt„Gs 



I 
2 
N 



{N-l){\nG-\k) + D-\k)G,{k)) 
\nG:\k) + D:\k)Gs{k) 



24 



24 



X(0) n x(0) 



GJk)+ . , 
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c^.(^)) --y^in(i-^n(.) 



+ ^^m\'' + iV^^;^ + ^ / [(iV - \)G^{k) + ^(A;)] 



+ 



2 

a 

12 



24 



[(iV-l)G,(A;) + 3G',(fc)] + 



G.{k) 



+ 



J^G,{k) J^Gsip)], 



(3.78) 



where we introduced the usual tree-level propagators for the sigma and pion fields 
as 

(3.79) 



Zi D 



If we forget about the counterterms, the expression in (3.78) coincides with the 
2PI effective potential obtained in [68]. We also could have figured out directly the 
diagrams present in the NLO large- iV 2PI effective potential of the 0{N) model, the 
auxiliary field however indirectly did this job. The terms above can be combined in 
a way which makes explicit that in the large- expansion there are only two bare 
couplings, namely -|- Sm^ and A + d\. When we would attempt to get a self- 
consistent solution of the equations arising from the variation of V2p/(f , Gtt, Gs), 
this would not be true (see the corresponding section later). We emphasize also 
that since the above functional is (9(A^°) accurate, in terms involving the sigma 
propagator Gg wc do not need the NLO part of the counterterms, i.e. 5m^*^' and 
5\'^^\ When differentiating with respect to G^, we should remember that also Il{k) 
is a functional of 6*,^! 

The interpretation of the terms in (3.78), which makes explicit the infinite series 
of diagrams summed up in the present treatment is as follows. The first three and the 
last three lines represent the 2P1 effective potential of the 0{N) model at Hartree 
level of the truncation and at NLO in its large- iV expansion. The remaining four 
terms represent the NLO contribution of the 2PI vacuum diagrams beyond Hartree 
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Figure 3.3: The two sets of vacuum diagrams which contribute beyond Hartree level 
and at NLO to the 2PI functional. Solid (dashed) line represents sigma (pion) prop- 
agator, while wiggly line represents the background v. The dots indicate any number 
of pion bubbles. 



level. The f ^-dependent part of these terms can be rewritten as 



A + 5A(°) 



6 



/ (k) 



(A + a(°))n(A;)/6 
1- (A + (5A(o))n(A;)/6 



A + aw 

6 



jGs{k)J2[^^^^m)) ■ (3.80) 



The v -independent part can also be written as a sum: 



(3.81) 



It is easy to show that disregarding counterterms, these terms correspond to the 
two sets of diagrams given in Fig. 3.3. Counterterm diagrams corresponding to the 
n — 2 term of the sum in (3.80) are displayed in Fig. 3.4. Similar diagrams with 
different number of pion bubbles can be drawn for the other terms appearing in 
the sums in (3.80) and (3.81). A direct way to obtain (3.78) consists of summing 
up all these diagrams with the associated combinatorial factors determined by the 
Feynman rules. 



3.2.2 Renormalizability checks on V2Pib, G^tt, Gg] 

The significance of (3.78) is that it displays all counterterms which guarantee 
the renormalizability of the resummation of the perturbative series, a resummation 
induced by the large- expansion. The finiteness of the equation of state and the 
self-energies obtained from its respective variations is ensured automatically, since 
this feature is "inherited" from the finiteness of the same quantities achieved in 
the formulation with the auxiliary field. Still, it is an instructive exercise to check 
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Figure 3.4: Diagrams corresponding to the n — 2 term in the sum on the r.h.s. of 
(3.80). A cross represents the counter-couphng 5\^^\ 



this feature directly. Exploiting the structure of our previous analysis done in the 
auxiliary formulation of the model we shall show the finiteness of the equations 
directly obtained from (3.78). 

The inverse pion propagator at NLO in the 1/N expansion is given then by 



where the nonlocal (E^) and local (M^) parts of the self-energy are: 



(3.82) 



P L 



= mo + 



1-XUf{p)/Q 3 (1 - Anp(p)/6)'. 



(^-^^0 ^'^'^^^ ~ ^ ^T'div(M2)(3.83) 



For the nonlocal part we used that — U — 6/X — Up/G and Tdiv is given in 

(3.41). 



Since the local part has LO and NLO contributions we write 



M^''^^ /N and expand the pion propagator to 0{1/N). With help of the integrals 
defined in (3.53) we obtain 
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\j^(Gs{k)-D^{k)) 



Dl{k) 



+ / D^{k) , (3.84) 



where in distinction to its definition given in (3.5) D^^ depends now on M^^"^. We 
used also that to leading order 11 (p) = lTr{p). 
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The equation for M^^*^^ in (3.84) is the usual gap equation at Hartree level of 
truncation of the effective action. This was analyzed in [51] and the counterterms 
which can be determined from this are SX^^^ and 6771"^^^^ given in (3.77). 



Observing that the left hand side of the equation for M^''^^ in (3.84) is finite, we 
obtain the following relation between counterterms and divergences: 
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Jdiv(M^^°)) - fdiv(M^(°) 
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Gs{k)-DM) 



(3.85) 



div 



Using the divergence analysis presented in the corresponding parts of the NLO 
renormalization, we have all divergences and integrals expressed in terms of M^^°^ 
for which we can substitute its finite gap equation M^*'"'* = + A(f^ + T^)/6. 
Requiring the vanishing of the coefficient of v"^ + determines 5\'^^\ while the 
remaining overall divergence determines 5m^^^\ Both arc in accordance with (3.77). 
The equation for the inverse sigma propagator obtained from (3.78) is 



iG-^{k) ^e-M 



2(0) 



A 



(0) 



(3.86) 



3 ^ l-Ag)4(W6' 

which is finite, since M^*^°^ is finite and 6/Ag^ — — 6/X — . 

Using the relation between the LO pion and sigma propagators we can show that 
the derivative of (3.78) with respect to v reads as 

SV2P1 , 



6v 



-[v,G^,Gs] = -NviG-\k = 0), 



(3.87) 



which is also finite, since we showed that G~^{k) is finite. It also displays the validity 
of Goldstone's theorem. 

We close this part by mentioning that the LO sigma propagator equation, the 
NLO pion propagator equation and the equation of state derived from (3.78) can be 
obtained also in the Dyson-Schwinger formalism, similarly as it was demonstrated 
in the auxiliary field formulation. Some details can be found in [73] (sec also [74] for 
the relation between the truncation of the Dyson-Schwinger equations and the 
expansion). This means that our investigation implies also the renormalizability of 
the Dyson-Schwinger equations at NLO in the large- iV expansion. 



3.3 Derivation of the IPI effective potential V[a^ v] 

Let us now go back to (3.73) and instead of eliminating the auxiliary field, let us 
focus on substituting the solutions of the propagators into V2Pi[a.,v,GT^,Q] (3.67). 
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In this case wc obtain the ordinary IPI effective potential (up to constants), which 
is defined as \^[q;,v] = — j^r[Q;, v]. 

Wc start from (3.67) and after using (3.72) and (3.73) we find that in view of 
(3.28) Ci? = Gs does not appear explicitly in the functional, which now becomes 

v[a,v,GT^\ — —[m^—tca)v +idKia-\ — —a 
2 2 As 



-kN-1) I {lnG-\k) + D-\k)G^{k)) 
--J In {k'-m^ + ta) (l - ^n(A;) j - 



Ad \ 2 



3 



(3.88) 



Here, the last term comes from (3.69). Next, we use (3.43), (3.41), and the definitions 
in (3.68) to write the inverse pion propagator as 

iG-\k) = k^-ml + tea - ^^n{k), (3.89) 

where S^(/c) is given by the integral of (3.43) calculated with the expression Gal 
taken from (3.30). Using this propagator in the first integral of (3.88), we easily see 
that when expanding it to 0{1/N), the contribution of the self-energy drops out 
and we are left with 

V[a, v] — —{rn^ — ica)v'^ + iSnia + ttt— 6^ — N- [ In (k^ — + ica) 
2 2Ab 2 

The radiative part of this functional has exactly the same form as the effective 
potential given in [19, 37, 38]. The difference in the classical part corresponds to 
slightly different ways of introducing the auxiliary field. More important is that the 
authors of [37, 38] restrict their counterterm functional only to terms proportional 
to pieces of the Lagrangian which are present already in the original formulation of 
the model. In their form of introducing the auxiliary field this restricts the countert- 
erms to those proportional to a and d^. By allowing all independent counterterms 
to appear which have dimension less than or equal to 4 in the auxiliary field for- 
mulation we will have enough fiexibility to ensure the renormalizability in arbitrary 
background. 

In order to demonstrate that (3.90) contains all the NLO counterterms, we sketch 
the renormalization of the SPE obtained by differentiating (3.90) with respect to a. 
Using that H, to be taken only at LO in 1/A^ expansion, depends on a through D^^ 
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defined in (3.5), we obtain 

- —a-iSK -i-(v^+ I ' 

\b ^ 2\ J — m\ -\- ica ^ 

-\ I (G2)(fc) - D^{k)) - ^ J(M^) + ^^.V(M^). (3.91) 

Here, we recognized the appearance of tlie expression of G^fJ wliicli can be read from 
(3.9). We used tlie relation (3.28) and for the last two terms also (3.71) and (3.53). 

All we have to do is to establish the connection between (3.91) and (3.51), the 
latter being already renormalized. The last three terms of (3.91) can be found in 
(3.51), if in that equation we take into account (3.52), so we have to work only on 
the first three terms of (3.91). Using the definition of the couplings, we expand them 
to 0{1/N) and obtains 



N f ^ [ i \ .{NX 
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2 V ' Jkk^-m% + icaJ " \2 ' 12' 

-^fa^{M^) J^Dlik), (3.92a) 



3A^ . , 3N ^ N I 

Xq X 2 



tP + (M^ - Mo^)ri°) 

+i5K^i^ + 25Ki^a. (3.92b) 



Since the last two terms of the second equality above coincide with 5A\^"'°(q;)/5q; 
of (3.51), the equivalence between (3.91) and (3.51) is demonstrated. 



3.4 Case of the complete 2PI-1/N resummation 

In this section we present the renormalization of the complete 2PI-1/N approxi- 
mation. By "complete" we mean that in the self-consistent equations derived from 
the approximate 2PI-1/N effective potential we do not expand further in 1/N for 
the propagators and the condensate but renormalizate within the self-consistent 
framework. 

Due to the complex self-consistent nature of the equations, we need a method in 
which the divergences resummed into self-consistent quantities can be extracted and 
eliminated. In the beginning of this section, we will review the method introduced 
in [44], in which these divergences are summed up by Bethe-Salpeter like ladder- 
equations in order to obtain counterterms. Although the method is elegant and well 
worked out, we will not have explicit, immediately calculatable counterterms as it 
was the case in the previous sections. However, we do not expect this at all, since 
we need to deal with more complex divergences as before. 
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Figure 3.5: Bethe-Salpeter equation of V . Ladder of A's are resummed. 



3.4.1 Theory 

Let us review briefly the results of [44]. Mostly we will not go into details, just 
state some useful parts of the referred paper. For detailed proofs, we refer the reader 
to [44]. 

As we have seen, in 2PI approximations it is possible to define the connected 2- 
point function in two independent ways. Considering the full theory, in the symmetric 
phase we had two identical representations of the self energy, check (1.92). If we 
denote the left hand side of (1.92) by — S, and the right hand side by — S, for 
approximate 2PI effective actions in general S 7^ S. 

Considering a theory with Z2 symmetry and staying in the symmetric phase, 
for the 4-point function we have three possible definitions, which are introduced as 
follows. Let us define the A function as: 



A(xi,a;2,a;3,X4) = 4 



5^T 



int 



5G{xi, X2)5G{x3, X4) 



(3.93) 



G=G 



where G is the solution of the propagator equation, as before. Let us denote by V the 
infinite series of ladder graphs made by A (see Fig. 3.5), which satisfies the following 
integral-equation (Bethe-Salpeter-like equation): 



V{xi,X2, X3, X4,) = A(xi, X2, X3, X4) + 



A has the following symmetry properties 



yzwu 



A{xi,X2,y,z)G{y,w) 

xG{u, z)V{w,u,X3,X4). (3.94) 



A{xi, X2, X3, X4) = A{x2, xi, X3, X4) = A{xi,X2, X4, X3) = A(x3, X4, xi, X2), (3.95) 



which is inherited by V. (3.94) can be written in a short-hand notation as: 



V = A + -AGV = a + -VG'^A, 



(3.96) 



where in the second equality we used the symmetry properties of the functions A 
and V. The first possible definition of the proper 4-point vertex is V. 
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Now we define the V function as: 

, (3.97) 



V(x X X X.)- ^^^^'^^'^ -2 ( 
[X,, X2, X3, X,) ^{x,)54>{x,) - A^{x,)A^{x,) [ SG{xs, 



int 
Xa) 



G=Gj 



where A/A0 refers to total differentiation, as before. It can be shown that V satisfies 
a very similar Bethe-Salpeter-like equation as y. In the previously introduced short- 
hand notation it reads as: 

V = K+l^VG'^K, (3.98) 

where 



A(X1, X2,X3,X4) = 2- 



(3.99) 



G=G 



b(\){xx)(\>{x2)?iG{x'i,Xi) 
It is possible to show that V can be expressed in terms of V: 

V = lv+\lv&V . (3.100) 



V is the second possible definition for the 4-point vertex. 

The last variant for the 4-point vertex is the most natural of the three possibih- 
ties: 



^(,) ^ 5^ ^ A^r[0,G[0]] _ ^^ ^^^^ 



Working out carefully the total derivatives, after a short calculation we arrive at 



r(4) ^ 



int 



G 



+ ^ (aG^V^ + perm. ) , (3.102) 



where V^{xi,X2,X3,X4} = V{x3,X4,Xi,X2), and the permutations must be taken 
in the last three arguments of F^^^ (there are two more possibilities). The most 
important fact is that it can be showed that in the symmetric phase of the exact 
theory r*^"*-* = V = V {in the broken phase there are additional terms in (3.98) and 
(3.101)), therefore V and V are as good definitions of the 4-point vertex as T^*K 
Similarly as it happened to the 2-point function, approximations of the 2PI effective 
action usually do not follow these relations. 

In [44] it is shown that with the renormalization of S, S and V, V, T^^\ the theory 
is completely finite, also in the broken phase. However, since in general S 7^ S and 
V ^ V ^ T^'^\ we have to introduce different mass and coupling counterterms 
corresponding to these quantities. Compared to ordinary perturbation theory, this 
is an unusual feature of 2PI approximations. 
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Figure 3.6: The counterterm diagrams corresponding to different definitions of tlie 
2- and 4-point functions. Neitlier wavefunction renormalization nor BPHZ countert- 
erms are indicated. A wiggly line represents f , while a plain one corresponds to the 
full propagator G. 



The conventional names of the counterterms which cancel overall divergences 
of these quantities are as follows. The Sm^ mass counterterm with the 5Zq wave 
function renormalization appear in the equation of S, therefore they can be found 
in the classical part of the effective action. 5m\ and 5Z2 are related to S, therefore 
they are included the tree level propagator of the 1-loop part ~ Tr Z^^^G. For the 
4-point functions, we introduce three different coupling counterterms: (5Ao, 8X2, 5A4. 
The first one, 5Xo is responsible for the renormalization of V , which means that it is 
related to the Hartree part ("eight" diagram). 8X2 renormalizes V , therefore it can 
be found in the ~ Ti D^^G piece, while 8X4^ makes F*^^) finite, it is included in the 
classical part. 

Since V and F'^'^^ can be expressed in terms of V ^ if the 2-point functions are 
renormalized, subdivergences can arrive only from A and A. In order to cancel these, 
it is shown in [44], that one has to apply the BPHZ method [8, 9]. This means that 
beyond the Hartree approximation we have to add corresponding counter-diagrams 
with unique counterterms related to A, and another one related to A. These are 
usually denoted by 5Xq^^^ and (5A^^^^, respectively. It is also shown that these 
counterterms are connected to 5Xo and 5A2, in general the latter ones include 5X^^^^ 
and 5Xf^^^ , therefore we can always introduce the following splitting: 

5Ao = 5Ao^™^ + AAo, 5X2 = 5X^^^^ + AA2. (3.103) 



In the next chapter, in a concrete example we will see the role of this splitting (see 
Fig. 4.1). 
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3.4.2 Effective potential revisited 

TBA 
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3.4.3 Propagator- and Bethe-Salpeter equations 

TBA 
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3.4.4 Renormalization 
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3.5 Concluding remarks 

We studied the renormalizability of the 0{N) model at next-to-leading order 
in the large-// expansion, at zero temperature. In the main part of the chapter, 
we constructed the 0{N^) counterterm functional of the model in auxiliary field 
formalism by studying the renormalization of the derivatives of the 2PI effective 
potential with respect to its variables. Expanding the propagator equation of the 
pion to 0{1/N) in the large-// expansion we showed that the renormahzation can 
be achieved for arbitrary values of the background and auxiliary fields, in a way 
that respects the internal symmetry of the model (e.g. Goldstone's theorem). This 
is expected on theoretical grounds, since divergences are determined only by the 
asymptotic behavior of the propagators and because any consistent resummation of 
the perturbation theory should resum also the counterterm diagrams associated to 
the perturbative series. Although, one can anticipate the consistency of the auxiliary 
field technique and of the large- expansion in dealing with perturbative series, 
the difficulty we face when trying to infer the renormalizability of the model in a 
given approximation from the fact that the model is perturbatively renormalizable 
is that one cannot easily keep track of what partial series of counterterm diagrams 
is actually resummed at a given order of the large- A" expansion. In consequence, the 
actual analytic check of the renormalization of a given approximation is unavoidable. 

We have also dealt with the elimination of the auxiliary field and the related prop- 
agators with keeping the dynamical sigma and pion propagators, which makes trans- 
parent the classes of diagrams containing also counterterms which are resummed in 
the original 0{N) theory at NLO of the large-A^ expansion. The explicit form of 
the counterterms is given for the first time for the theory using auxiliary field and 
also for the original formulation. In the original theory, the 2PI efi'ective potential 
contains only two counterterms, a coupling and a mass counterterm, both having 
LO and NLO parts. 

The two examples we worked out (i.e. V2Pi[f , Gtt, Gs] and \/[Q;,f]) demonstrate 
that the renormalizability of the broadest functional V2Pi[a,v,GT^,Q] implies the 
renormalizability of the functionals arising after the elimination of some subset of the 
variables. This result obtained at zero temperature can makes us confident that the 
renormalization goes through for T 7^ as well and that the counterterm functional 
determined here will prove helpful for phenomenological studies in the 0{N) model. 
Also, the method developed here can be used also to models with more complicated 
global symmetries. 

The final part of this chapter contained the complete 2PI-1/N renormalization 
of the 0{N) model at NLO, i.e. when we looked for a self-consistent solution of the 
propagator equations and the equation of state. In this case no further expansion 
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in 1/N was considered after deriving the equations for the 1- and 2-point functions 
from the NLO 2PI effective potential. The renormahzation was done in the sym- 
metric phase of the modeL In order to obtain well-defined and finite results also 
at the broken phase, beyond the 2-point functions we had to renormalize the 4- 
point functions as well (in the previous case, we renormalized directly the equation 
of state). As expected from [44], unlikely from our previous results, here we found 
that in order to renormalize the complete theory, we had to include various mass- 
and coupling counterterms regarding to different definitions of the 2- and 4- point 
functions, which is an unusual feature of (self-consistent) 2PI approximations. 



Chapter 4 



Numerical implementation of 2PI 
renormalization methods 

In this chapter our goal is to investigate the numerical applicability of the 2PI for- 
malism. This is motivated by the fact, that in an equilibrium setting and beyond the 
Hartree approximation (i.e. the "eight" diagram in theories with 4-point coupling) 
there are relatively few papers reporting on numerical solutions of renormalized 2PI 
equations, even in scalar models. 

The application of the 2PI method did not went much beyond the demonstration 
of its features in the simplest models (e.g. the A'"-component real scalar model with 
0{N) symmetry, in many cases restricted to N — 1). The explicitly finite self- 
consistent propagator equation of the 0^ model was solved including the setting- 
sun diagram at vanishing field expectation value and at finite temperature in [42]. 
The renormalized 0{N) model was solved at zero temperature and at vanishing 
field expectation value within the bare-vertex approximation of the auxiliary field 
formalism in [62]. The pressure of the one component 0^ model was calculated in 
Euclidean space and finite temperature in [29], and of QED in [30]. Solving the 
renormalized field and propagator equation of this model in Minkowski space the 
phase transition was studied in [49], where it was found that including the field- 
dependent two-loop skeleton diagram at the level of the 2PI effective potential results 
in a second order phase transition. It was proven analytically in [100] that in Hartree 
approximation the phase transition cannot be of second order. The study of thermal 
properties of the spectral function of the renormalized 0"^ theory in the symmetric 
phase was reported in [98]. In some phcnomenologically oriented studies [39, 99] the 
necessity of renormalization was neglected, e.g. in [39] the 0{N) model is solved in 
the Minkowski space with some approximations by taking into account momentum- 
dependent corrections only in the imaginary part of the self-energy and treating the 
tadpole with a UV cutoff. 
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Especially in the above listed finite temperature studies, little technical details 
are given on the numerics, in particular the rate of convergence of the iterative steps 
leading to self-consistent solutions. This fact does not allow us to easily infer the 
accuracy of the methods used in obtaining the results. In our opinion, the lack of 
standardized, well-tested algorithms might be the main reason for the rather scarce 
applications of 2P1 approximations in the study of the thermodynamics of relevant 
field theoretical models. 

Our aim in this chapter is to quantify the accuracy of the iterative solutions 
of 2P1 approximations and to compare different algorithms. We also would like to 
make connections between different renormalization methods. Our highly accurate 
solutions could serve as a benchmark for methods to be used for a finite temperature 
solution of these equations. 

In this chapter we stick to the (p"^ theory. We will derive the usual coupled set 
of self-consistent propagator equation and the equation of state. We will discuss 
in details their renormalization by adapting the renormalization method of [92, 47] 
to the broken symmetry phase (i.e. with nonvanishing field expectation value). We 
isolate the quantities which need to be renormalized, and imposing renormalization 
conditions, we recover the set of counterterms derived using the minimal subtraction 
method introduced in the previous chapter and which was also used in [50]. On the 
other hand, we derive explicitly finite equations as well. We will present different 
algorithms for solving the set of the finite equations and also the equations containing 
explicitly the counterterms. The solution will be given at zero temperature and 
in Euclidean space. We demonstrate that their numerical solutions coincide and 
compare the efficiency of the two methods. We show that in the case of solving the 
equations containing explicitly the counterterms, faster convergence of the iterative 
algorithm is obtained if the counterterms themselves are derived iteratively and 
develop parallel to the quantities to be determined (i.e. the propagator and the field 
expectation value). 



4.1 The 2PI effective potential at two- loop field- 
dependent level 



The 2PI effective potential of the one-component real 0^ model at the field {v) 
dependent two- loop truncation level is given in the following form [44, 94, 49, 50]: 
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4.1. The 2PI effective potential at two- loop field-dependent level 
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which can be obtained directly using the rules of 2PI formalism discussed in Chap- 
ter 1. The tree-level propagator is D{p) — i/{p^ — — \v^/2). The counterterm 
functional is given by 

where the different terms are represented graphically in the previous chapter on Fig. 
3.6. As already discussed there, the origin of two different mass counterterms and 
three coupling counterterms stems from the fact that within the 2 PI formalism it is 
possible to define two two-point functions and three four-point functions, which do 
not necessarily coincide within a given truncation of the 2PI effective potential. 

The stationarity conditions 5V[v^G\/ 5G = and 6V[v,G]/6v = give self- 
consistent equations for the full two-point function and the equation of state: 

iG-^(p)=/-E(p), (4.3) 

= v(m^ + dml + ^{X + 6X^y + ^ (-^ + S^2)T[G] + y 5(0, G)^ , (4.4) 

where the self-energy E(p) is given by 

E(p) = m^ + 5ml + hx + 6X2^ + ^(A + 6Xo)T[G] + ^AV/(p, G). (4.5) 

Zi 

Note that we call E(p) the self-energy, although its usual definition used in the 
2PI formalism does not contain the mass parameter and the corresponding mass 
counterterm. The tadpole integral TfG], the bubble integral /(p, G) and the setting- 
sun integral at vanishing external momentum 5'(0, G) which appear above are defined 
as follows: 

T{G\ = / G{k), (4.6a) 
I{p,G) = -i f G{k)G{k+p), (4.6b) 
S{Q,G) ^ -i [ I G{k)G{q)G{k + q). (4.6c) 

Jk Jq 

Just as discussed in Chapter 2., [50] also calculates the divergences of these integrals 
by expanding the full propagator around the auxiliary propagator 

with the mass Mq playing the role of the rcnormalization scale. The counterterms 
absorbing these divergences were obtained in [50] by requiring the separate cancel- 
lation of the divergent coefficients of i/K and of the environment (temperature 
and the v background) dependent finite function Tp[G], representing the finite part 
of the tadpole integral, both in the propagator and field equations. 
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4.2 Renormalization conditions and count erterms 

In this section we adapt the renormaUzation method used at finite temperature 
in [92] to the broken symmetry phase at zero temperature. We mention that the 
renormahzation of scalar models displaying spontaneously broken symmetries was 
discussed generally in the context of the 2PI approximations in [47]. Here, we sep- 
arate the field dependence from the divergent quantities in the same way as the 
temperature dependence was separated in [92] and give the renormalization condi- 
tions which, when imposed on the divergent quantities, lead on the one hand to 
exactly the counterterms determined in [50], and on the other hand to explicitly 
finite equations, with no reference to any counterterm at all. We note that the coun- 
terterms in question have exactly the same structure as those used in [49] and which 
were derived using the same minimal subtraction procedure introduced in Chapter 
3.^ In the next section we will check explicitly that the solution of the propagator 
equation and the equation of state obtained using explicit counterterms agrees with 
the solution of the explicitly finite equations to be derived below. 

4.2.1 Renormalization of the propagator equation 

We begin by splitting the propagator into "vacuum" and "matter" parts. 



where the vacuum part refers to a propagator defined in the symmetric phase (f = 0), 
while the matter propagator includes all explicit and implicit dependence on v. The 
truncation of the 2PI effective potential studied in this work leads to a momentum 
independent self-energy for f = which means that the vacuum propagator can be 
simply parametrized in terms of an effective positive mass Mq, as in (4.7), so that 
G^^\p) = Go{p). Note that we need this propagator with positive mass squared Mq 
because in the broken symmetry phase m^, the renormalized mass parameter of the 
Lagrangian is negative. 

As a consequence of the splitting (4.8), the self-energy (4.5) is decomposed into 
three parts: 



where the explicit expression of the different parts are given below. The vacuum 
part is independent of the momentum, depends only on the vacuum propagator and 
has divergence degree 2. The last two pieces of the self-energy appearing in (4.9) 

^Unfortunately the counterterms used in [49] appeared only on the corresponding poster, with 
no details on the calculations leading to their results. 



g'(p) = gW(p) + g(-)(p), 



(4.8) 



e(p) = e(^) + e(°)(p) + e«(p). 



(4.9) 
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are v dependent and are defined as follows. The part which does not result in any 
divergence when the full propagator is expanded around G^''^ is called regular part 
and is denoted by T.'^^\ while the remaining part, having divergence degree 0, is 
denoted by The splitting of the field-dependent part is somewhat arbitrary. 
Using (4.8) in (4.5), the different parts of (4.9) are identified as: 

= Ml + 6ml^ + ^^ jG^^\k), (4.10a) 

E(°)(p) = - Ml + 5ml^ + '^K^\q,p) + J^G^-\k), (4.10b) 

E«(p) = -iAV / G^'^\k)G^''\k+p)-i— [ G^'^\k)G^'^\k+p),{4:.10c) 

Jk 2 

where A^^^(0,p) is defined as 

aW(0,p) = A + 5A2-iA2 / G(")(A;)G'(")(A; + p). (4.11) 

Jk 

We will see that with the above choice for the difi'erent pieces of the self-energy, 
the expression of the counterterms given in [50] can be obtained through simple 
renormalization conditions. The mass counterterm 5m| was split into two parts: 
5m2 A is responsible for the renormalization of the vacuum part, while 5m2 Q has to 
remove the divergence generated by the E(°) dependence of G^™). In order to see 
how this latter divergence proportional to E^°) emerges, we expand the propagator 
around the vacuum part: 

^ i i (, I SW(p) + EC-)(p) \ 

p2 _ _ s(o) {p) - S W {p) p2 _ 5](v) + p2_ Y,iy) J 

= G^^\p) - z(G'(^)(p))'s(°)(p) + ^^(p). (4.12) 

The term proportional to E^"^ is 0(l/p^) (up to logs), therefore it gives divergent 
contribution after integrating over the momentum, while the regular part of the 
propagator goes with 0{l/p^) and leads to a finite contribution upon integration 
over the momentum. The sum of these two terms is the previously introduced matter 
part 

G^'^\p) = -i(G'(^^(p))'E(°)(p) + G«(p). (4.13) 

As already announced at the beginning of this subsection, the first renormalization 
condition is 

E(^) = Mj, (4.14) 
which determines ^ through the relation 

K. = -^rfi, (4.15) 
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where is defined in (3.12). 

Now wc turn to the rcnormahzation of E^^^. After plugging (4.13) into (4.10b), 
we obtain the foUowing integral equation: 



2 

We can search for the solution of Yj-'^^ in the following form: 

E(°)(p) = r(:) + ^r(^\0,p) + \v^^^ jG^\k). (4.17) 

Using (4.17) in (4.16) we find, that T^Z\ r^''^(0,p), and T^^) fulfill the following 
Bethe-Salpeter type equations: 

rL^^ = m^-Mo^ + KB-^^^rWy^(G(^)(^))', (4.18a) 
rr)(0,p) = Ar^(0,p)-i:^±^^(GW(A;))¥r)(0,fc), (4.18b) 

rw = A + ao-i^^rw /(GW(;t))', (4.i8c) 

2 Jk 

in which all integrals are divergent. To render these divergent equations finite, we 
have to impose some renormalization conditions, which will also determine the ap- 
propriate counterterms. Due to our introduction of the Mq scale and the splitting of 
the mass counterterm into to pieces, we need one more renormalization condition in 

^J]{v) 

addition to (4.14) to fix also ^m^^. By using (4.10a) we see that (m^ — M^] 
satisfies 



'0 



-I [m 



which is the same equation fulfilled by , up to counterterms. Since 

ml = <*-2») 

a simple renormalization condition on T^m is 

= - Ml (4.21) 
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This fixes the (5m| ^ counterterm: 

Sml, = (m^ - Ml)^-^^ = -^(m^ - M-)Tf \ (4.22) 

where Tf^ is the famihar logarithmic divergent quantity defined in (3.14). From 
(4.15) and (4.22) we obtain the complete 5m\ counterterm. As the sum of 5m|^ 
and 5m2_£, it leads to 

= + (r<=> + - Ml)Tr) . (4.23) 

Before moving to the renormahzation of (4.18b) and (4.18c) we note that Sm^ can 
be equivalently fixed to the same value by keeping the condition (4.21) but replacing 
(4.14) with the condition 

where "overall" refers to the environment independent part of the self-energy, which 
does not coincide with H^'"^ due to additional terms coming from E*^°\ which contain 
some overall divergences. We can identify this part of the self-energy as the one 
which remains after v and G^^^ are formally set to zero. Then, using the condition 
(4.21) in the right hand side of (4.18a), we obtain in (4.17) 

^^'\P)L. = - ^0 + Kb + (A + 5Ao)(m2 - Ml)Tf^ /2, (4.25) 
so that upon adding it to E*^"^) of (4.10a), the condition (4.24) becomes 

+ 8ml + + {m^ - Ml)Tf'^) = m\ (4.26) 

Next, we focus on the renormalization of and r*^^\ As a renormalization 
condition we impose 

r(^^ = A, (4.27) 

which using (4.18c) leads to 

= ao + ^A(A + ao)Tj°^ (4.28) 
Then, using (4.18b) and (4.11), we see that the difference 

^X'lPip) (4.29) 

is finite, where we defined l''p\p) as the finite part of the bubble integral in the zero- 
momentum subtraction scheme in which Ip^ (0) = 0, as in the previous chapter. This 



r(^\o,p)-n^^(o,o) = -iA 



(v). 
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shows that we only need to renormahze the momentum independent r2 (0, 0) which 
is achieved by imposing on it the renormahzation condition 

r^"\0,0) = A. (4.30) 

Then, the exphcitly finite r2^^(0,p) function is given by 

r^;\Q,p)^X + X'lPip). (4.31) 

Using (4.31) in (4.18b), we obtain the equation 

= ^2 + ^A(A + SXo) (rf + AT^)) + A^Tf , (4.32) 

where 

T^^-i f{G^^\p)yiP{p). (4.33) 

In conclusion, the imposed renormalization conditions give explicitly finite ex- 
pressions for the functions r^''Ho,p), and appearing in (4.17), and provide 
also the following finite equation for T,^^\p): 

E(°)(p) = - + ^(a + X'l^ip)) + \jG^'\k)- (4-34) 

Two interesting remarks are in order at this point. The first concerns the solution 
of (4.18b) which, as we can check iteratively, is easily expressed in terms of F^^^ as 

rr\0,p) = A(^)(0,p) - ^rW^A(^)(0,^)(G(^)(fc))^ (4.35) 

The second comment refers to obtaining a different expression for the 8X2 countert- 
erm. Using (4.27), (4.30), and (4.31) in (4.35) we obtain 

= <5A2 + ^A(3A + 5X,)Tf^ + ^A^ [{Tf^f + T^) . (4.36) 

Note that (4.36) does not coincide with (4.32), but the solutions for 8X2 do. Indeed, 
using the result of (4.28) for 5Xq in (4.32), the same expression as that coming from 
(4.36) is obtained. The equation (4.36) served in [50] as a consistency check of the 
renormalization procedure. 

Before we proceed to the renormalization of the equation of state, let us make 
a remark. So far in the thesis, renormalization was a more algebraic and formal, 
than a pictorial procedure, illustrated by diagrams. To help the reader getting closer 
to the more familiar BPHZ method, we also present a diagrammatic illustration 
of the cancellation of one of (sub) divergences. We choose the part proportional to 
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Figure 4.1: Diagrammatic renormalization of the I'^-dependent part of the self- 
energy. See the text for details. The symmetry factors are not indicated and we 
used the shorthand notation Aq = A + ^Aq. The boxes denote the divergence of the 
encircled part of the graph, the plain line denotes the propagator G^^\p) = Go{p), 
the wiggly one represents v, while 'F' denotes the finite part of the bubble integral. 
The index I in the sums goes over the number of bubbles (for the two classes of 
diagrams the lowest number of bubbles is zero or one, as indicated) . 



v"^ in Eq. (4.10b) for E*^^'', see Fig. 4.1. This example shows the role the coupling 
counterterms play in the removal of both sub divergences and overall divergences. 
The first line in Fig. 4.1 shows the u^-dependent diagrams generated upon iterating 
Eq. (4.10b) after the BPHZ subtraction procedure was implemented (see [44] for 
details.) The counterterm SX2 displayed in Fig. 3.6 is decomposed into the sum of 
two terms (JAf^^^ and AA2. ^Af^^^ is used to absorb the divergence of the bubble 
integral present in A2^''(0,p). The remaining finite part of the bubble integral is 
denoted by T'. In order to obtain the second line of the figure we localized all 
the subdivergences. Then, we use the renormalization condition (4.27). To obtain 
the third line of the figure, a cancellation between the third and fourth term of 
the second line was exploited and the consequence of the renormalization condition 
(4.27) (which is A - Aq = (AoA/2) j^G^"'\k), where Aq = A + 5\o) was rewritten 
in a diagrammatic form. The remaining overall divergences of the diagrams in the 
third line are removed by AA2. Since the divergent integrals associated with these 
diagrams are Tf^ and t'^^\ keeping in mind that we also have to add (^Aa^^^, from 
the condition of vanishing of the square bracket in the last line of Fig. 4.1 we obtain 
the expression (4.32) of 8X2 determined previously. 

At this point, we obtained the counterterms 5m|, SXq and 8X2. These are deter- 
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mined by (4.23), (4.28) and (4.32), respectively. Explicitly finite equations were also 
derived for E^^^ and E^*") (and therefore for the complete self energy), these are (4.34) 
and (4.10c). 

4.2.2 Renormalization of the equation of state 

We proceed with applying a similar procedure on the equation of state (4.4). We 
use again the splitting of the propagator introduced in the previously in (4.8) and 
obtain 

l2 



Jp ^ Jp Jk 

-i^ JJ^isG (p) + G (p)) G {p + k)G (k) 



A + SXa n 



+ 1 I A^^\0,k)G^'''\k). (4.37) 



K 



2 

Using (4.13) and (4.17) in the second hne, we obtain 

= m' + Sml+^^^^ fG^^\p)-i^ [ f G^''\p)G^^\p + k)G^^\k) 

^ Jp ^ Jp Jk 

-^rW / Ar)(0,^)(G(^)(A;))V^rf) + ^ / V^\Q,k)G^^k) 

ipL ^ ^ ^ {P + k)G ("^^ (A;) , (4.38) 
where we used the expression of T^^ given in (4.35) and defined 

r5") = A + 5A4-^^ / A2{0,k){G^''\k))Y;\0,k). (4.39) 

^ Jk 

The integral in the expression above is divergent, and we impose the renormalization 
condition 



^ = A. (4.40) 



This leads to the equation 

= SX, + lx{X + 6X2 + A^Ti°)) (Tf) + AT^)) + ^A^ (t]^) + ATj^'^^) , (4.41) 
where 

Tf^) = -z/(G(^)(p))^(4^)(p))l (4.42) 
Jp 

Since the last but one term of (4.38) is already finite in view of (4.30) and (4.31), 
and the last term is a finite integral, we are left with the divergences of the first line 
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and the first term of tlie second line. In order to complete the renormalization we 
only have to find the quantity on which the renormalization condition determining 
Sml can be imposed. It is straightforward to check that 

A-'V2Pi[v,G,] 



AvAv 



-i^JJ^G,=o{p + k)G,=o{p)Gv=o{k): (4.43) 

where G^ on the left hand side is the f-dependent solution of the stationarity condi- 
tion SV2Pi[v,G]/SG\a=Gv — the subscript indicates explicitly that it depends 
on the vacuum expectation value v, so that the chain rule has to be applied when 
taking the derivatives with respect to v. With the splitting of the propagator G„=o 
into "vacuum" , "matter" , and "regular" parts, we can proceed exactly as at the be- 
ginning of this subsection, but we have to omit everywhere the terms proportional 
to v'^. Then, we have 

I 2 Jp 

-i^ [ [ G^^\p)G^^\p + k)G^^\k) 
6 JpJk 

-i^^Z^ I (0, k) {G w + (0, k)Gi%{k) 

JJ^ (3GtUp) + g2(p))g'2(p + (4.44) 



AvAv 



A2 



The last term above is finite and since r2^''(0, k) was already made finite by imposing 
the condition (4.30), the last but one term above is also finite. Therefore, in principle 
we could choose our last renormalization condition as ^.J^i^h^l = ^2 then 

AvAv lv=0 ' 

(5mQ determined through this condition would differ by a finite term from 6m'^ given 
in [50]. In order to obtain the desired expression for Sml we choose the following 
condition 

A-'V2Pi[v,G, 



AvAv 



= m\ (4.45) 

v=0, overall 



where the concrete expression on the left hand side is given by the right hand side 
of (4.44) without the last two terms, and the "overall" notation indicates that the 
divergence of this term is the usual overall divergence, that is the divergence which 
remains after all subdivergences are removed. In other words, the quantity on the 
left hand side of (4.45) is obtained formally by setting v and G^^^ to zero after (4.8), 
(4.13), and (4.16), (4.17) are used in the explicit expression of ^J^p^!:^ Then, it 
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is easy to check that (4.45) gives 



= 5mg + -(A + 5A,)[Tf + (m2-Mo2)Tf 



(4.46) 



where 



G^^\p)G^''\k)G'^^\p + k). 



(4.47) 



J p J k 



The two renormahzation conditions (4.39) and (4.46) together with the expression 
of r2^^(0,A;) given in (4.31) lead to the exphcitly finite version of the equation of 
state 



With this the renormahzation is complete. We obtained equations for the re- 
maining counterterms Srri^ and 5A4, these are (4.46) and (4.41). The finite version 
(4.48) of the equation of state is also known. 

The imposed renormalization conditions give exactly the same expressions for 
these counterterms (see (4.23), (4.28), (4.32), (4.41) and (4.46)) introduced in the 
functional (4.2) as those appearing in [50]. This demonstrates the equivalence be- 
tween the two renormalization procedures. We note that even though we formulated 
our set of renormalization conditions on the quantities F'-^-', r2^'*(0,p), T^^ \ corre- 
sponding to the three four-point functions V ^ V and F'^^-' introduced in the last 
section of the previous chapter, and also on the curvature ^ ^'»a1''^"^ ^"^^ ^^^^ 
energy E, they do not look as natural as those imposed in [44, 49] at T = T* and 
in the symmetric phase of the model (f = 0). (This method was the one adopted to 
the renormalization of the complete 2P1-1/N approximation of the 0(N) model in 
the previous chapter.) In our case G*^^^ is not the full propagator at w = and in 
order to obtain given expressions of the counterterms (which contains G in place 
of the full propagator), we needed a specific way to formulate the renormalization 
conditions in (4.24) and (4.45). We imposed these conditions in order to facilitate the 
comparison between the solution of the finite equations and those which use coun- 
terterms (sec next section). On the other hand, as discussed in Ref. [100], finding a 
set of natural renormalization conditions, prescribing e.g. the mass defined from the 
variational propagator and the curvature at the minimum of the effective potential 








p J k 



(4.48) 
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(as we know, these two quantities do not necessarily coincide in a given truncation 
of the 2PI effective potential) together with the value of the independent four-point 
functions is problematic at T = 0, that is in the broken symmetry phase, even at 
the Hartree approximation of the 2PI effective potential. 



Before discussing the algorithms and presenting the results we make some general 
statements about the numerical method we use. 

Since we solve iteratively the set of coupled integral equations consisting of the 
self-consistent propagator equation and the field equation, we have to store and up- 
grade in each iteration step the value of some functions on a grid and with the use of 
interpolation, approximate the value of these functions at the points required by the 
integration routine. For this purpose we use an equidistant grid for the modulus of the 
Euclidean four-momentum, the one-dimensional Akima spline interpolation method 
and the numerical integration routines of the GNU Scientific Library (GSL) [101]. 
When we solve the set of finite equations, we store the "regular" and the "matter" 
part of the propagator (see (4.13)), while when solving the set of equations with 
the countertcrms it is the full self-energy which is stored on the grid (see (4.9) and 
(4.10)). In the first case the integrals are convergent in the UV, but nevertheless, for 
the functions which are stored on the grid we have to take into account the fact that 
there is a maximal value of the modulus of the four-momentum. The convolution 
of these functions should be calculated accordingly, using ^-functions to restrict the 
momenta of the propagators: 



where 9 is the step function, for the moduli we introduced p = \pe\, Q = \kE+PE\- L 
is the maximal value of the modulus stored on the grid. The g-integral is a remnant 
of the angular integration in a 4d spherical coordinate system and is obtained with 
the change of variable cos^^i = (g^ — — p'^)/{2kp). 

For functions of the momentum with known analytical expressions (e.g. the vac- 
uum propagator), there is no need to store them on the grid, and in consequence 
a simplification is encountered in the convolution involving such functions. In this 
case there is no need for the second 9 function, the entire angular integration can 
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(4.49) 
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be done analytically, and with the change of variable t — tan(^i/2) we have 

^ 

+ kl + 2\pe\ \kE\ cos^i + M2 



TT 



4« l^^l +Pe + K - \/ {kl + pI + M',Y - iklplj . (4.50) 

In contradistinction to the previous type of functions, when counterterms are 
used, all momentum-dependent functions arc cut at the maximal value of the mod- 
ulus of the four-momenta, that is at the physical cutoff A, irrespective of the fact 
that they are stored or not stored on the grid (in the latter case their expression 
is known analytically). We mention here that cutting the sum of momenta in the 
case of the bubble integral and as a consequence in all the double integrals encoun- 
tered seems to be the cleanest way to proceed. In fact, in the case of solving the 
set of finite equations we could afford to simplify the calculation, by cutting only 
the loop-momenta in the integrals involving the propagator G^^\P) = Go{P) which 
is not stored on the grid, because the integrals are all convergent in the UV, and 
in consequence the corrections are suppressed. We will see indeed that the solution 
of the finite equations nicely agrees with that obtained by solving the equations 
with counterterms. Moreover, cutting only the loop-momenta, when solving the set 
of equation containing the counterterms, produces some small differences compared 
to the case when the sum of momenta is also cut, at least up to the largest cutoff 
we investigated, and in consequence is not a viable method for obtaining accurate 
numerical results. 

4.3.1 Algorithm for solving the finite propagator equation 
and the equation of state 

We will solve in Euclidean space the exphcitly finite equation of state (4.48) and 
the propagator equation pieces (4.34) and (4.10c) obtained before. A Wick rotation is 
performed in every integral appearing in our equations: with the continuation ko 
ik^ for example the tadpole integral reads G{ko, k) — )■ i J^,^ G{ik'^, k) = J^^ A{kE) 
where the Euclidean four-momenta is A';^ = (k, k^), such that k'^ — li^ + k\ and the 
Euchdean propagator is defined as A{kE) — i/{k% + I^{kE))- 

Using the expansion (4.12) of the propagator around the vacuum propagator 
A^^^ks) = 1/(^1; + Mq) and the definition (4.13), the matter and regular parts of 
the Euclidean propagator reads 

A^^^HM = .2 , ..2 , ^ , w.w. . -yAyJI^ (4-51) 



kl + M^ + m{kE) + ^('KkE) k% + M^ 



A<'.(*.) ^ A<»-.(*.) + ^1^. (4.52) 
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From (4.34) and (4.10c), the two pieces of the self-energy which appear above are 
given by 

E(°)(M = m'-M^ + x(l + XlP{kE)y-^ + ^j^ AW(M, (4.53) 

\ 2 2 

E«(A;£;) = X''v^I^'^^\kE) + ^I^""^\kE), (4.54) 
where the "matter-vacuum" and "matter-matter" bubble integrals are defined as 

7(-v)(^^) ^ _ /■ a(-)(ps)AW(pb + M, (4.55a) 
l(^^\kE) = - / A(^\pE)A(^\pE + kE). (4.55b) 

The finite part of the vacuum bubble integral, Ip\ appearing in (4.53), is defined 
in (4.29). After analytical continuation to Euchdean space, using (3.23) and (3.24) 
it has the following expression (recall that Go = G^^^): 



iPikE) = 



+ 2 



(4.56) 



Note that a self-consistent equation for the regular part A*^"") can be obtained in 
principle if wc substitute the expression of A("^) from (4.51) into (4.52) because 
E^°^ contains the integral of A*^'') and A*^™) appearing in the integrals of S*^''-' can be 
expressed in terms of A^'') and Nevertheless, we store on the grid two quantities, 
A*^™) and A^^\ which have to be solved simultaneously with the equation of state 
(4.48), which in Euclidean space reads 



o ^ 



[ (l + Xli:\kE))A^'\kE) + ^(SS^"""^"^ + )), (4.57) 



where the corresponding pieces of the setting-sun integral are defined as 

-S^"^"^") = - / / ^^'^\vE)A^'^\kE)A^''\vE^kE), (4.58a) 

Jve J^e 

^(mmm) _ _ /" /" A^"^) (pfi) A^'") (A;^;) A^"^) (p^; + A;^;) . (4.58b) 
J vw. Jkp. 



Using (4.50) in the first term of (4.54) and for the setting-sun contribution to the 
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equation of state (4.57) which contains one vacuum propagator we obtain 

X (^kl + p2 + m2 - ^(A;|+p2 + M2)'-4A;|p2^ , (4.59) 

^(mmv) _ __1_J^ dkkA^'^\k)J^ dppA^'^\p) 

X ^A;2 + p2 + - Y^(A;2+p2 + M2)^-4/tV^ • (4-60) 

For the second term of (4.54) and the setting-sun contribution to (4.57) which con- 
tains only matter propagators we use (4.49) and the results of the angular integra- 
tions are 

1 fL ;>min(p+|A;£;|,L) 

7(mm)(^^) = -— ^ / dppAM(p) / dqqA^^\q) 

87r-^k% Jq J\p-\kE\\ 

x^Ap^q-^- {p'^ + q^-kl)\ (4.61) 

^(mmm) _ dk k A^"'\k) dp p A^'^\p) 

pirLm(p+k,L) 

X / dqq^Ap-^q^ - (p2 + g2 _ A;2)2A('^)(g). (4.62) 



Our iterative algorithm for solving the coupled set of equations for A*^'"\ A^™\ 
and f goes as follows. Wc start at zeroth order by neglecting all the integrals in the 
self-energy and the field equation, including the finite part of the bubble I^p\ By 
doing so, from (4.57), (4.51), and (4.52) we obtain: 



A' ^0 (M- (A;2 +Mo2)(A;| + m2 + |^2)' 

^^^^ = ,.2 \.2,2?.2^ "Ia 2V 4.63 

(/c|; + M^Y{k% + m? + f t;2) 

Then, starting from the first order, the iteration is done in two steps. At a generic 
order n > 1 wc first upgrade Ar™"* and A^'^'' on the grid using (4.51) and (4.52), 
respectively, where the self-energy, which includes now J^^\ is calculated with the 
lower order quantities: fn-i, A\J^, AJ_-^. As a second step of the nth iteration, we 
upgrade the vacuum expectation value by calculating f„ from (4.57) with the in- 
tegrals done with the already upgraded An™"* and An''. This procedure is repeated 
until the iteration converges. 
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4.3.2 Algorithm for solving the propagator equation and 
equation of state using counterterms 

If we decide to solve the propagator equation and the equation of state using 
counterterms, whose role in this equations is to cancel the divergent part of the 
integrals, then we need a method to determine them. Such a method was developed 
in [50], and we review it here because it will be slightly modified for numerical 
reasons. In terms of the remaining finite part of the tadpole Tp [G] and bubble integral 
/f(P, G) the propagator equation (4.3) reads 

iG-\p) ^p'-M'- ^X'v'hip, G), M'^m' + + ^Tp[G]. (4.64) 
Using the expansion of G around Go 

G{p) = Go(p) - iGlip) (^M' - Ml + \\\^h{p, G)^ + . . . , (4.65) 

the tadpole, bubble, and setting-sun integrals can be explicitly decomposed into 
divergent and finite parts (see Ref. [50] for details): 

I{p,G) = Tf^+h{p,G), (4.66a) 
T[G] = rf^ + (M2-M2)ri°^ + ^AVry^+rF[G], (4.66b) 

>S(0, G) = 5o(0) + 3Tf Tf[G] + 3(M2 - Ml) [i^ff + 

^-\\\- (rfTf + rf^)) + S,iO, G), (4.66c) 

where the divergent integrals Tj^\ T^^-* , \ Tj^ '^^ and 5'o(0) were already defined 
in (3.12), (3.14), (4.33), (4.42) and (4.47), respectively. Plugging the decomposed 
expressions of the integrals (4.66) into the propagator equation and the equation of 
state, after subtracting the corresponding explicitly finite equations from them, writ- 
ten in terms of the renormalized parameters m^, A and the finite part of the integrals, 
we obtain relations between the counterterms and the divergent integrals. Requir- 
ing that the divergent coefficient of the v^, v^, and Tp[G] vanish independently, we 
obtain the counterterms determined in [50], which were also obtained from appro- 
priate renormalization conditions previously (see (4.23), (4.28), (4.32), (4.41) and 
(4.46)). Actually, these requirements lead to six equations, instead of five. The extra 
condition arises when the outlined procedure is applied to the coefficient of f ^ in the 
equation of state. We obtain (4.36), but as stated below this equation, its solution 
for 6X2 coincides with the one obtained from (4.32), as it should. 

In Euclidean space these counterterms are functions of the 4d rotational invariant 
cutoff A appearing through the integrals can be calculated explicitly. They read as 
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follows: 
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(4.67a) 
(4.67b) 
(4.67c) 

(4.67d) 
(4.67e) 



where for Ca[AW, A(^)](A;) we have to use (4.49) with /\^^\k) = l/{k'^ + M^) and 
with L replaced by the physical cutoff A. 

The counterterms Sm^, 5m|, 5Ao, 8X2-, 8X4^ can be used to solve iteratively the set 
of Euclidean equations 

E(p) = + + ^(A + 8X2^ + \{X + 5Xq)T{Y) + ^AV7(p, E), (4.68a) 

Zi Zi z 

= + + i(A + ^4)^;' + ^(A + 5A2)r(E) + ^5(0, E), (4.68b) 

in which they have to cancel the divergences of tadpole, bubble and setting-sun 
integrals, which in Euclidean space read as 

r[E] l\kk'A{k), I{p,E) = -Ca[A, A](p), 



5(0, E) = 



87r2 



dA;A;^A(A;)CA[A,A](A;). 



(4.69) 



However, since the counterterms are designed to cancel divergences of integrals 
produced by the full (i.e. iteratively converged) propagator, they cancel the complete 
cutoff dependence of T[S], I{p, S) and S{0, S) only when the solution of this set of 
equations has converged. This is not a problem in itself, since, as it will be demon- 
strated, the iterative procedure converges, but the convergence of the solution of 
the set of equations is rather slow. It is possible to improve the iterative procedure 
if, instead of using the obtained counterterms, we rcdcrive them at each order of 
the iteration using the procedure outlined above. These counterterms, which come 
from the requirement to have finite equations at each order of the iteration, will 
evolve during the process of iteration toward the value of the counterterms already 
obtained. Since, as we will see, the intermediate counterterms are now environment 
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dependent, that is they depend on v and TVfS], this convergence can only occur, if 
V and S also converge to their appropriate values, which has to be checked numer- 
ically. The final countertcrms then become environment independent and coincide 
with the ones determined by (4.23), (4.28), (4.32), (4.41) and (4.46). 

In what follows we present the improved iterative procedure which uses evolving 
counterterms. At each order the iteration starts by the upgrade of the self-energy, 
followed by the upgrade of the vacuum expectation value determined from the equa- 
tion of state. 

At zeroth order of the iteration the field and the self-energy are equal to their 
tree-level expressions: 

^0^ = -^, Mp)=ni' + ^vl (4.70) 

At this order, since there are no quantum fluctuations at all, there are no divergences 
to cancel and in consequence all the countertcrms arc zero, that is SuIq q = g = 
^^0,0 = SX2fi = 6X4^0 = 0, where the counterterms also carry the index of the 
iteration number because all of them will change during the iteration process. 

The general formulas for the nth order counterterms to be given below encompass 
all orders with the following prescription: at n = we start from the tree-level 
values, and formally T(S„) = /(S„) = for n G {—1,-2} and = vq. Suppose 
that the (n-l)th order of the iteration is done, i.e. v'^_^ and E„_i(p) together with 
the corresponding counterterms are already known, then at nth order we start by 
upgrading the self-energy using (c.f. (4.68a)): 

+ ^-^^n-2 + ^-^^(Sn-l) + Y^n-JyP, ^n-l) ■ (4.71) 

The finite part of the tadpole and bubble are defined through 



+^^ry^+ri^(Sn-i), (4.72) 
7(p,E„_i) = + 7^^(p,E„_i). (4.73) 



(0) 
d 



The renormalization procedure requires the independent vanishing of the overall 
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divergence, field-dependent subdivergence and Tp-dependent subdivergence: 








6ml „ + 



(Tf + (m^-Mo^)Tr), 
2 



SX2,n + (A + 5Xo,n) + ^TP(X + 5Xo,n)(l - Snl) 



+XX 



2rri0) ^n-l 



V, 



n-2 



(A + 5Ao,„)^Tf T^(E,_2) + 5Ao,nTF(S„_i). 



(4.74a) 

(4.74b) 
(4.74c) 



The counterterms Sml,^, 5X2,n and 5Xo,n are obtained from (4.74), so that S„(p) 
is finite and can be calculated. Note that for n = 1 the second term in the equation 
for 5A2,n vanishes as it should because in this case there is no divergence due to the 
bubble, that is there is no T^^^ in T'(Eo), whose expression is 



(4.75) 



Note also that for n = 1 we have (5Ao,i = 0, since Ti?(E_i) = by convention. 

We continue with the determination of from the equation of state, which at 
nth order of the iteration reads 







g T(E„) + ^5(0,E„), (4.76) 



where ^(E^) is given by (4.72) with n — 1 — > n and the setting-sun is 

5(0, E„) = So{0) + 3rf rF[E„] + 3(m' + ^vl_, + ^T^(E„_i) - M^) 



X 



T!.y + ri") + iA^„L. (ri"rj»> + t^) + s,{o, s„). (4.77) 



Applying the same requirements as for the self-energy, we obtain: 



+y(rn^-Mo^)((Tr)^ + Tf) + ^5o(0), 



= 



3A 



A2 
6 

(0) 



a4,„ + y (A + 6X2,n + A^rf )) (ri°) + ATi^)) 



+ 



3A3 



AT. 



(/,2) 



Tf(E^ 



X 



A 



2(A+a,„)Tr+^((Tr)^+Ty)) 



(4.78a) 



(4.78b) 



(4.78c) 



Note that a new counterterm 6X2 has been introduced instead of 6X2, since in the 
nth iteration 5X2,n does not eliminate properly the appropriate subdivergences of 
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the equation of state. Extracting Sm^j^, SX^^n and SX2,n from equations (4.78) we 
can obtain from (4.76). The counterterm 5A2,n will be equal to 5A2,n deter- 
mined from the propagator equation only when the iteration converged, so that 
TfC^ti-i) /TFC^n) — ^ 1- In this limit and when Vn/vn-i 1 the counterterms will 
converge to the values obtained from the renormalization conditions. This happens 
only asymptotically and in practice the iteration stops when the relative change 
in the values of the field and the self-energy are smaller than some given value dic- 
tated by our requirement of accuracy. The same iterative procedure is used when the 
counterterms arc not evolved during the iteration, but used as determined from the 
renormalization conditions: we first upgrade the self-energy and then the equation 
of state and check for their convergence. 

4.3.3 Numerical results 

In what follows we present the results on the numerical solution of the propa- 
gator equation and the equation of state in Euclidean space. As it was stressed, we 
wanted to build up a demanding numerical framework, in which very highly accurate 
solutions can be obtained. This concerns in particular the equations containing the 
counterterms explicitly, since in these the increasing cutoff requires very large can- 
cellations, which is hard to be handled numerically. Previous attempts on obtaining 
the cutoff independence of renormalized quantities for large cut offs (e.g. [49]) could 
be only moderately achieved, therefore such a numerical framework in question also 
carries theoretical importance as well. 

When solving the explicitly finite set of equations the maximal value of the 
modulus of the four-momentum stored on the grid is L/ ^/2\m?\ = 500, while in the 
case of solving equations with counterterms the highest value of the modulus stored, 
that is the 4d physical cutoff, is A/ ^y2\m'^\ — 200. In the numerics, we choose |2m^| 
as the unit mass, while the value of the renormalization scale is Mq/ ^/2\m^\ = 2.1. 
Actually, to have an idea on the physical scale, we can require to have v/T.^{p ~ 
0) ~ 0.66, which is the value of the ratio fn/triT^ of the pion decay constant to the 
pion mass for /^r = 93 MeV and = 140 MeV, and T,^{p = 0) = m^r. Then, from 
the left panel of Fig. 4.2 we see that the first requirement can be meet by choosing 
A ~ 8, for which value the second one selects y^2\m'^\ ~ 158 MeV. In units for which 

= —0.5, the step size on the grid in general was 0.01. We have checked that by 
halving the step size the change in the results is comparable with or smaller than 
the precision required by our convergence criterion. The iteration was stopped when 
the change in the vacuum expectation value was smaller than 10~^. The precision 
of the numerical integration routines used was higher than that, since we required 
relative error bounds between 10~^ and 10~^. 
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Figure 4.2: The solution of the exphcitly finite equations. Left paneh v/M and 
S5(p = 0) as function of the couphng obtained for L/^^/2\m?\ = 500. Right 



M 



paneh The converged value of v as function of the maximal value of the modulus 
of the momentum stored on the grid for different values of the coupling constant A. 
An appropriate constant was subtracted from the values of v obtained for A = 3, 10, 
and 20 such that the resulting values coincide at Lj 
V obtained for A = 30. 



500 with the value of 



First, let us discuss the solution of the exphcitly finite equations. The field ex- 
pectation value and the self energy at zero momentum can be seen in Fig. 4.2. In 
the right panel we see the converged solution of the equation of state obtained for 
various values of the coupling as function of L, the maximal value of the modulus 
stored on the grid. As explained in the figure caption, v was shifted appropriately 
in order to make all four curves meet at L/y^2\m^\ — 500. In all cases a plateau 
can be observed for increasing values of L. However, with increasing value of the 
coupling A the plateau starts at a higher L. With the same convergence criterion 
the number of iterations increases with increasing A and for A = 30 we need twice as 
many iterations until convergence occurs compared to the A = 3 case. Interestingly, 
at a given A, the number of iterations does not depend on L. 

Next, let us focus on the numerical solution of the equations containing explicitly 
the counterterms. We show, that we can get the very same results for the field 
expectation value and the self energy, as which obtained from the finite equations. 
Wc also quantify the extent of the achievable improvement given by the use of the 
iterative ly evolved coupling and mass counterterms, as compared to the solution of 
the equations which use the counterterms originally derived in [50]. In order to see 
this, we compare in Fig. 4.3 the number of iterations needed for the convergence of 
V at four different coupling constants. Wc use a reasonably large cutoff A/ ^y2\m'^\ = 
200, around which the convergent results are expected to show cutoff independence 
within some desired accuracy. We observe the advantage of evolving the counterterms 
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Figure 4.3: A comparison between the convergence of the iterative algorithms using 



fixed and evolving counterterms for a fixed value of the physical cutoff: A/ W2\m 



200. The value of f / ^'l\m^\ is shown as a function of the number of iterations for 
different values of the coupling. 



even for lower values of the renormalized coupling. For larger values (A > 10) the 
improvement is significant, as we can see from the large iteration number needed 
by the algorithm using fixed counterterms. Therefore, it is expected in general that 
solving self-consistent equations iteratively in a fully non-perturbative regime with 
fixed counterterms is not efiicient from a numerical point of view. For large couplings, 
A > 10, the number of iterations needed in case of using fixed counterterms are at 
least a factor of two larger than the number of iterations needed for the explicitly 
finite equations to converge. 

In Fig. 4.4 we compare for A = 5 the ways how a cutoff independent solution 
is reached using fixed and evolving counterterms. The left panel shows that with 
increasing values of A, a plateau in the w(A) functions is approached rather slowly, 
when the fixed counterterms of [50] are used, despite the fact that the coupling is 
not large. Looking at the plot with values of A spanning a large range, after several 
iterations a plateau is finally reached. However, the structure in f (A), better seen in 
the inset which uses a smaller A-range, shows that the convergence is not uniform: 
more iterations are needed for larger cutoffs to get close to the true solution, that is 
the solution of the explicitly finite equations (horizontal line). This non-uniformity of 
the convergence in A is more pronounced for larger values of the coupling constant. 
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Figure 4.4: A comparison between the convergence rate of the iterative algorithms 
using fixed (left panel) and evolving counterterms (right panel) as reflected in the 
change of the vacuum expectation value at different values of the physical cutoff A. 
The horizontal line represents the solution of the explicitly flnite system of equations 
obtained with L/ ^J2\iv?\ = 500, where L is the maximal value of the modulus of the 
momentum stored on the grid. The insets show the convergence at large iteration 
numbers. 




Figure 4.5: The convergence of the full self-energy obtained for A = 5 with evolving 
counterterms for A/ ■\j2\m?\ — 200. The inset shows the behaviour of the self-energy 
at large momenta. 



Contrary to this, the right panel of Fig. 4.4 shows a uniform convergence in A when 
the iterative algorithm uses evolving counterterms. Since by the very construction of 
the evolving counterterms a plateau is observed already for moderate values of the 
cutoff at every order of the iteration, we see again that with this improved approach 
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the final, convergent result is obtained in a more efficient way. 

The convergence of the self-energy can be seen in Fig. 4.5. The function converges 
smoothly, i.e. for all momentum values the self-energy reaches its converged value 
almost simultaneously. Note however that, the convergence for ~ is a slightly 
slower as compared to larger values of the momenta. 

4.4 Concluding remarks 

In this chapter we presented an approach which can be used to obtain a very 
accurate numerical solution of the self-consistent propagator equation coupled to 
the equation of state, both derived from the two-loop level approximation of the 2PI 
effective action in the 0^ model. We showed that the renormahzation method devel- 
oped in [50], which removes the divergences using the minimal subtraction scheme 
introduced in Chapter 3 is equivalent to imposing nontrivial renormahzation con- 
ditions on the self-energy, the curvature of the effective potential and on kernels 
of Bethe-Salpeter equations related to these quantities. The use of renormahzation 
conditions also allowed us to construct explicitly finite propagator and field equa- 
tions. 

We compared the convergence of the iterative method apphed to solve, on the 
one hand, the explicitly finite equations and, on the other hand, the equations which 
contains the exphcitly calculated counterterms of [50]. The very same results were 
obtained numerically, which proves the correctness of our analytic study on renor- 
mahzation. It turned out that, especially for larger coupling constants (A > 10), 
the convergence rate of the algorithm which uses the counterterms is quite poor. 
This should not come as a surprise given that the counterterms were determined 
based on the asymptotic behavior of the propagator. Within an iterative procedure 
the asymptotics develops progressively and is reached only when the solution con- 
verged. Only then the divergences present in the equations match the counterterms 
constructed to cancel them, at intermediate steps of the iteration process the use of 
fixed counterterms results in an oversubtraction. In order to cure this, we managed 
to develop a different algorithm which rederives the counterterms at every step of 
the iterative procedure from the condition that the counterterms cancel the diver- 
gent part of the integrals at every order of the iteration. This method, which uses 
evolving counterterms turned out to be very effective in obtaining the converged 
solution within a few iterative steps also for larger couplings. 

The presented approach and some of the methods applied in our present work 
may be directly extended, first to a zero temperature numerical solution of the 
model in Minkowski space, then to a finite temperature study. It is expected that the 
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solution which can be obtained with the method described here could serve as a good 
benchmark for other, even more effective methods used in finite temperature studies 
(see e.g. [29, 49]), and will represent a good basis for obtaining highly accurate and/or 
numerically controlled solutions of the 2PI approximations in more complicated 
theories as well. 



Chapter 5 



An approximate solution of the 
U(N)xU(N) model in large-N 

Universal features of finite temperature and finite density variation of the QCD 
ground state realizing the U{3) x U (3) approximate chiral symmetry were investi- 
gated with the help of the corresponding meson model [77], introduced already in 
Chapter 2, for arbitrary flavor numbers. There were numerous attempts to treat 
quantitatively the flnite temperature restoration of chiral symmetry of the strongly 
interacting matter in the framework of this model. A common goal of these investi- 
gations is to find out the nature of the phase transition, when the baryonic density is 
varied. Another central issue is the description of the quark mass dependence of the 
finite temperature symmetry restoration [78, 79, 80]. For its investigation one couples 
constituent quarks carrying baryon number to the meson model [81, 82, 83, 84]. 

Optimistically one could say, that the location of the characteristic points of the 
QCD phase diagram (i.e. the T — ^ plane) determined with different variants of the 
model and in different approximations do agree with each other and with the results 
of the lattice field theoretical simulations within a factor of 2. In particular, improved 
agreement with lattice determinations of the QCD phase diagram were reported, 
when the Polyakov loop degree of freedom is coupled to the quark- meson model [85, 
86, 87]. In realistic approaches the effective models are all strongly coupled, therefore 
one usually experiences large variations in their predictions, when the simplest mean- 
field treatments are improved by taking into account quantum fiuctuations of the 
mesons and the constituent quarks [55, 88]. 

As discussed in Chapter 2, more recent investigations for the three-fiavor meson 
model in the large-A^ approximation are missing from the literature. Although the 
N = oo solution of the 0{N) model is the textbook example of the application 
of the 1/N expansion [2], the leading-order solution of the large- iV approximation 
of the U (N) X U{N) model is still unknown. There were various attempts to han- 
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die the model in the large- iV limit, but none of these went beyond the related 
0(2A'"^)-symmetric nonlinearities (i.e. corresponding to the ~ (s")^ + (vr")^ part of 
the Lagrangian). A first step in this direction was made with the matrix general- 
ization of the Hubbard-Stratonovich transformation involving the complete quartic 
potential of the U{N) x U{N) symmetric model [108]. Based on this reformulation 
of the model, a saddle-point solution was constructed using an auxiliary matrix field 
[109]. It was conjectured that this solution describes the leading order behavior of the 
ground state for large flavor number. The approximate solution was applied directly 
to the thermodynamics of the three-flavor effective meson model (sec also [83]). In 
the actual solution of the matrix saddle point equation the authors restricted their 
Ansatz to a condensate proportional to the unit matrix, which reduced once again 
the quartic potential to its 0{2N'^) invariant part. As a result, the part of the vacuum 
condensate pointing in the algebra space along the direction of the longest diago- 
nal generator remained proportional to the corresponding component of the explicit 
symmetry breaking external source. To the author's actual knowledge, no published 
further step was attempted towards the exact large- analysis of the U{N) x U {N) 
symmetric matrix model. 

The numerical details of the pattern of chiral symmetry breaking in theories of 
three-flavored strong interactions are essentially dictated by the observed mesonic 
spectra [102, 78, 103, 104]. Potentials of the effective models are constructed with 
the requirement that in the ground state a condensate proportional to the 3x3 unit 
matrix should arise spontaneously, and split the degeneracy of the parity partner 
states. This fundamental structure is slightly modifled by an external source also 
proportional to the unit matrix, which generates mass for the pseudo-Goldstone 
pseudoscalar mesons. Guided by these aspects, in the following section flrst we will 
use the symmetry breaking pattern U{N) x U{N) U{N), corresponding to a 
condensate of the dynamical field which is proportional to the unit matrix. However, 
a possible parametrization of the linear sigma model has to refiect also some realistic 
mass information concerning the pseudoscalar sector. This leads us in the second 
half of this chapter to a more general condensate when U{N) x U{N) U{N — 1) 
breaking can be realized. This is induced by the presence of the part of the external 
field proportional to the longest diagonal generator and results the kaon-pion mass 
splitting via the formation of a supplementary condensate along the same direction of 
the Lie algebra. It is a rather nontrivial question if this component of the condensate 
arises fully from the external field, or can evolve spontaneously. 

The goal of this chapter is to describe an approximate leading order (A^ = oo) 
solution of the linear sigma model with U (N) x U{N) symmetry, and to obtain the 
ground state at zero temperature and density. Although the solution takes into ac- 
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count only two-loop contributions of the 2PI effective action, it definitely goes beyond 
the restricted 0(2iV^) symmetric solution. For the restriction of the contributing di- 
agrams it exploits in addition to the large- expansion an assumption in the mass 
spectra of the model, concerning the relation of the scalar and pseudoscalar sector. 
The region in the coupling space, where this assumption is valid can be estimated 
by determining the spectra from the approximate solution self-consistently. In the 
first part of the chapter we also present the construction of the renormalized version 
of this solution in some details, and provide an illustrative investigation of its range 
of validity. 

The idea to impose an extra assumption on the mass spectra stems from the way 
one deals with the 4 scalars and 4 pseudoscalars defining the U{2) xU (2) symmetric 
meson model. As we have seen in Chapter 2, there one simply omits half of the fields, 
i.e. the 3 components of the scalar-isovector triplet and the pseudoscalar-isoscalar 
singlet with reference to their higher mass. We shall assume an analogous feature to 
occur in the spectra of the leading order solution of the U (N) x U{N) meson model 
at large-A^. Then, introducing a more general symmetry breaking, there will be 
splittings in the pseudoscalar masses as well. Based on the "heaviness" assumption 
concerning the scalar mesons, in a first approximation to the = oo solution we 
retain only the quantum fluctuations of the light pion flelds. We do not attempt the 
anomalous realization of the Ua{1) symmetry, therefore in the first limiting case all 
pseudoscalars will have the same mass. 

After finding the propagators of the scalar fields, we can calculate corrections 
to the pion propagators and condensates arising from the heavy scalar fiuctuations. 
It is possible then to upgrade also the scalar propagators and in principle, we can 
iterate this procedure until the solution of the large- 2-loop 2P1 approximation is 
reached. It will be demonstrated, that the pion fields still obey Goldstone's theorem 
after the scalar corrections are included. We shall explore the divergence structure 
of the equation of state, the pion self-energy and the saddle-point equations (SPEs), 
and determine the counterterm pieces in the effective action necessary for their 
renormalization. 

For the second part of the chapter we address the question of the determining 
the ground state of the system. We will describe a general symmetry breaking, then 
turn to the case corresponding to the U{N) x U{N) — > U{N — 1) pattern. Then 
we construct the 1-loop effective potential and show that in a large part of the 
parameter space there are two local minima of the effective potential, in which non- 
trivial n-scalings are found corresponding to the field expectation values. This will 
result that instead of the usual power series in n^, we will also found some parts of 
the effective potential which are of odd powers of n. 
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5.1 Renormalized solution under the assumption 
of heavy scalars 

Let us recall the U{ri) x U{ri) symmetric Lagrangian from (2.13) and extend it 
with the explicit symmetry breaking term in the zeroth scalar direction. Taking into 
account the n scalings of the couplings, we have 

C = l(^5^s«5'^s" + aX'9''7r''-m2(s"s'' + 7rV)) 

- ^ fs"s" + 7r"7r"V - ^?7"C/" + /ioS°. (5.1) 
4n"' V / 2n 

Before proceeding, we introduce two auxiliary (composite) fields representing the 

combinations of C/" = ^dabds'^s'^ + tt^'tt'^) — fabcS^T^"^ and + tt^tt". This can be 

achieved by adding the following constraints to the Lagrangian: 



2 , , I X 2 



A£ = -l(x-y|L(,V + .V«)) -l(y«-,Y|ir.J , (5.2) 

and treat X and as variables. In the sum £ = L + AL, the U{n) xU {n) U {n) 
symmetry breaking pattern corresponds to the shifts 

s" ^ + V2^v5''^, U" 2^A^vs" + nV2n5''%'^. (5.3) 

After the introduction of the auxiliary field variables and the shifts, the full La- 
grangian has the following form: 

L = ^ [{d^s'^f + {dXf - (^n'^v'' + 2\/2^^;s° + (s«)' + (tt")') ] + 
- ^-X^ - ^(y«)2 + iy^x(2nV + 2^2^^^;^° + (s")^ + (Tr")^) 

+ iW— r"(C/" + 2vrivs" + ^2^5'^%^) + V2^/i(s° + ^2^^;), (5.4) 
V n 

where \^2n?h = Hq. Next, wc shortly describe the assumed structure of the solution 
and introduce some notations. The classical constraint equations show that the back- 
ground V induces nonzero X, y° values and introduces mixing of the pair s^,X and 
also of s",y" for every value of the index a. Wc construct correspondingly a quan- 
tum solution, where the saddle-point values of X, y° are nonzero, = 0, a 7^ 0, 
and only the following 2-point functions do not vanish: 

[GgOX, GgOyO, GxyO, GgOgO, GyOyO, Gxx], 
[Gs^s^,u=vj^Oj Gg^yvy^^yjLQ^ Gyuyvy^^^jLg]^ G^"7r'',w=i;- i^-^) 

For arbitrary fields A{x) and B{y), the 2-point function Gab refers to {0\T{A{x)B{y))\0). 
As in the previous chapters, we represent these propagators in Fourier space. The 
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Figure 5.1: Vacuum diagrams contributing to the leading order equation of (cf. 
the "setting sun" integrals of (5.11) containing and proportional to n^). One pion 
line is cut corresponding to the variation of the 2PI effective potential with respect 
to the pion propagator. The factor of 2 in front of the first diagram indicate the 2 
ways of cutting the vr lines. 



sets in square brackets form mixing sets of fields: there is a 3-dimensional mix- 
ing sector, and there are — 1 identical copies of 2-dimensional mixing two-point 
functions. The vr sector is diagonal. The equations below show degeneracy of the 
2-point functions for m = f 7^ 0, therefore it is convenient to introduce the following 
short-hand notations: 

gOgO = UgO, LrxX = (^X, CryOyO = CryO, Lr-^O^o = G^o, 
G-j^u^u = GjT) GgUgu = Gg^ G guyu = G gY ^ Gy^Y^ — Cry^ U 0. (5.6) 



5.1.1 Approximate 2PI effective potential 

In this subsection a 2PI effective potential is given, from which equations of the 1- 
and 2-point functions can be obtained by functional differentiation. The functional 
depends on the background v, the auxiliary (composite) fields (X, y°) and the 2- 
point functions. We use the same approximation as in Chapter 3: only the setting-sun 
type of diagrams are kept. Since we have more fields and vertices, there are a lot more 
of these diagrams. Since we perform a leading order large-n analysis, only those parts 
of the expressions are needed, which contribute to the leading order of the equations. 
For example on Fig. 5.1 those diagrams are indicated which are needed for the pion 
propagator equation. 

We introduce the following notations: 



Q{s,Y) 




G{X,Y'^,sO) 



^ Gx GxYO GgOx'^ 

GxY'^ GyO GgOyO 
\GgOx GgOyO GgO j 



(5.7) 
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The tree level versions of these quantities are: 



2iy/g^v -1 



(x,yo,sO) 



/ -1 






-1 



\ 2iy/g{v 2iy/g^v iDq ^ / 

(5.8) 



where we used the notation 



77/ 

which turns out to be the tree-level pion propagator: 



(5.9) 



(5.10) 



— i- 



— i- 



The approximate 2PI effective potential without counterterms reads as: 
^2Pi = n^{m^v'^ - 2hv) + ^ (^^ + {Y^f^ - in{^2iiX + ^2^1"°)^^ 

lnG;^(A;) + D-\k)G,{k)) (lnG'-^(A;) + D-J{k)GAk)) 

Tr[ In e;^- V) ik) + {k)g^s,Y) {k)] 
- i\ I Tr[lne?(3^^^^o,,o)(A;) + I^(l^_yo,,o)(A;)^?(x,yo,,o)(A;)] 

J k 

-(72 I [{G,{k)G^{k + p) + Gs{k)Gs{k + p))Gy(p) 



n 



+ iy^2y j GsY{k)GsY{k+p)Gs{p) 



^(?2 I (G'y(A;)G',(A; +p) - G,y(A;)G',y(A; + p))g^(p) 

i / (G.(A;)a(A; +p) + G^{k)G^{k + p))Gx{p) 
iy/^2jJ^(G,{k)Gs{k+p) + G^{k)G^{k + p)')GxYo{p) 
^-92 ^ (G.(^)a(fc + p) + G,(^)G.(fc + p)) Gyo (p) 
+ 2t,/^ I I GsY{k)Gs{k+p)G,ox{p) + 2tg2 1 I G,Y{k)Gs{k + p)G,OYo{p) 

J k J p J k J p 

i^g2 ^ J (Gs{k)Gy{k+p) + GsY{k)GsY{k+p))GAp) 

(5.11) 



+ i 
+ i 
+ i 



+ i 



+ 192 GY{k)G^{k+p)G^o{p). 

Jk Jp 
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Note, that since the tt and s — Y sectors are of multiphcity (n^ — 1), some 0{1) 
terms in (5.11) are not relevant for their leading order equations. Also, there are more 
of 0{1) terms of the complete effective potential, which contribute to the tt and s — Y 
sector only at NLO, these are omitted from (5.11). However, for the remaining fields 
it is necessary to take into account the proper 0{1) terms (otherwise we would have 
trivial equations for the corresponding propagators), which are therefore included 
in (5.11). This means that (5.11) is not a large-n expanded 2P1 effective potential 
truncated at 2-loop level. It is a potential from which leading order equations can be 
derived in each sector of the theory. Note, that the structure appearing in the seventh 
line corresponds to a setting-sun diagram with antisymmetrized vertex functions (i.e. 

~ fabc)- 

Wc stressed in Chapter 3 (during the renormalization of the 0{N) model), that 
the setting-sun approximation of the (unrenormalizcd) 2P1 effective potential corre- 
sponds to a specific truncation of the tower of the Dyson-Schwinger equations. Here 
an analogous statement can be made. If we close the tower at the level of the 3-point 
functions and substitute their classical value to the equations of the 1- and 2-point 
functions, we obtain exactly the same equations which can be derived from (5.11). 
These nonzero classical 3-point couplings are: 

^x-,:-^^, y, z) = Fxs-sbix, y, z) = i—^5abS{x - y)5{x - z), 

TYasbgc{x, y, z) = VYa^b^c{x, y, z) = —dabJ{x - y)S{x - z), 

V n 

Tyasb^olx, y, z) = -iJ—fabcS{x - y)S{x - z). (5.12) 

The derivatives of the unrenormalized effective potential with respect to the 1- 
point functions are: 

— X — m - ' - I 



r° - in^2v'' + i^n (^j^ Gs{k) + ^ G^{k) 



sx 

^V2PI /7r~..2 , 192 

5Y° 

^^2PI _^2(..2^. U\ r-2 



6v 



2n^{M\ -h)- 2i^n^ I Gsvik). (5.13) 

Jk 



When we look for the physical values of the 1-point functions (mean fields), coun- 
terterms must be included and the derivatives must be set to zero. In this case, 
comparing the solutions of the equations for the two auxiliary fields, we find the 
relation: 

= (5.14) 
which should also be valid after renormalization. 
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Now we turn to the propagator equations. Zeros of the derivatives of the effec- 
tive potential with respect to the 2-point functions give the appropriate equations. 
In the pseudoscalar sector the following equations are found at leading order (the 
abbreviated notation J = Ji^G[,]{k)G[,]{p + k) is used): 

iGj = iD^^ - 2ig2 j GyG^, 

iG-' = tD^'-tg2 J GY{Gs + G^)+tg2 J GsyGsy. (5.15) 

We note the potential violation of Goldstone's theorem when comparing G^{k — 0) 
or GT^o{k — 0) with the equation of state. Nevertheless, our forthcoming approxi- 
mate solution obeys this theorem due to a rather nontrivial relation between the 
relevant tadpole and bubble contributions. The couphng between s" — Y"- can be 
seen explicitly in the structure of the corresponding sector, where (if a 7^ 0) we find 
— 1 identical mixing 2x2 equations: 

iG:^ - iD^^ - ig2 j GsyGsy - 192 j Gy{Gs + G,), 

iG;^ = 2i^2V - ig2 j Gsy{Gs - G^). (5.16) 
In the 3x3 mixing sector of {X, y°, s°), the following 6 equations are obtained: 
iGx' ^-1-igi J {GsGs + G^G^), iG'^yo = -^^/gIg^ J (G.G, + GA), 
^^so^x = '^^Vdiv - 2i^gig2 j GgGgY, iGyl = -1 - ig2 J (GgGg + G^G^), 
iG^YO = 2i^v - 2ig2 J GsGsY: iGj = iD^^ - 2ig2 J {G^Gy + G.yG.y). (5.17) 

5.1.2 Counterterm functional 

In view of the equations for the 1-point functions and the propagators, we col- 
lect all the counterterms, which should ensure the finiteness of these equations. As 
before, we introduce the most general counterterm functional, which is allowed by 
2PI formalism. 

In this case it means that we can introduce counterterms to all independent 
pieces appearing in the mean-field part of the effective potential and independently 
of them also to the terms showing up in the inverse tree-level propagators. Note, 
that we might have to introduce counterterms also to pieces which are chosen to 
have an apparently fixed numerical (zero or unity) renormalized coefficient. In this 
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sense the most general form which we will need and is allowed by the structure of 
V2P1 is the following: 

1 
2 



Ct 



0\2 



+SxyXYQ 



— I 



n 



(5yyGy(A;) + 5xxGx{k) + 5yOyoGyo(A;) + 25xyoGxyo(^) 



n 



(^Ztt/c^ — Sm^ + i- 



n 



n 



n 
n 



(5.18) 



Note, that the self-energy corrections to G~y turn out to be finite, therefore no 
counterterm is needed there. In the second line, countercouplings proportional to 
XY^,Y^,X are also needed despite the fact that the corresponding renormahzed 
values are chosen to be zero. Similarly, countercouplings in the third line belong to 
the contribution of the (pure) auxiliary propagators occurring in terms of the type 
TrT)~^Q, in which all their renormalized coefficients are fixed to unity except the 
coefficient of Gxy", which is zero. In the fourth line, countercouplings are introduced 
into the same type of terms, they appear in Dj^ and D~^. As we know from Chapter 
2 and 3, in the exact solution of the theory only unique quadratic and quartic 
counterterms should occur, but in any finite order of a 2PI approximation we have the 
freedom to choose the countercouplings in the counterterm functional independently 
[44]. For the determination of (5.18) we should analyze the divergence structure of 
the integrals appearing in the propagator equations. 

Prom now on, we make an assumption on the mass hierarchy of the model. 
We assume that the scalar sector is considerably heavier than the pionic, which 
simplifies the solution of the limiting form of the coupled 1- and 2-point equations 
valid at n — )■ cxD. In addition, we also assume that the scalar masses are considerably 
larger than the amplitude of the symmetry breaking vacuum condensate. Since each 
component of the propagators in the mixed sectors will have a common denominator 
displaying the corresponding heavy mass, the only consequent way to neglect the 
heavy sector in a first approximation is to erase the bubble contributions containing 
at least one component of G(s,y) or ^(x,yt',s")- Then only bubble diagrams exclusively 
built with arc included. This means that in (5.15) all bubbles arc suppressed, 
while in (5.16) and (5.17) only pure pion bubbles remain. In this case all pion 
propagators arc equal to their tree-level value, and therefore have the same mass 
(check (5.15)). 

The equations of the 1-point functions are not that simple, they also need renor- 
malization. In the forthcoming procedure we use the scheme introduced in Chapter 
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3, i.e. we expand the propagators around an auxiliary perturbative propagator Go 
with mass Mq and identify and subtract the divergences in terms of the integrals of 
Go- The explicit form of the saddle-point equations can be written with the help of 
splitting the pion tadpole into finite (Zjf ) and divergent parts as (see (3.18)) 

/ G,{k) ^ [ Do{k) = T: + Tf + (M^ - Mo2)Tf . (5.19) 

Jk Jk 

With the use of the counterterms, it becomes the following: 

- = i (2v' + + (M^ - Mo^)rf + T^) - (6,x + + 2i5,o: (5.20a) 

91 ^ ^ ^ y/9l92 ^ 

^ = ii^v" + Tf + (M^ - Ml)Tf + - (5,/ + + 2i5,o. (5.20b) 

92 ^ ^ ^ \/9i92 ^ 

Here we introduced the notations 

/2^^ ^,0 V2^,.o 



ix, |/° = 1^Y\ (5.21) 



X — - 

n " n 

and we note that the counterterms ?)xv,^xi:i^xsi^yvi^yK ^-nd 8ys can be chosen zero, 
therefore the sake of transparency we did not write down them explicitly in (5.20). 
Taking into account the definition of M^, the following counterterms renormalize 
both equations: 

8x = ^y = T^, Sxy = ^/(ng~2Tf\ 

-^^--\{Tf' + {m'-Ml)Tf^). (5.22) 



Because of the omission of the last ("heavy") term in the third equation of (5.13), 
the equation of state does not require any extra counterterm. The finite equations 
for the 1-point functions read as follows {5mQ does not receive any contribution at 
this level): 

-=i(2^;2 + rj), ^=i(2^;2 + rj), M\^h. (5.23) 
9i 92 

In the similar approximation, the 2x2 propagator matrix of the (s, Y) sector 
simplifies to 



(5.24) 



where 7,^ is the usual bubble integral: lTr{k) = —i Do{p)Do{p + k) and we have also 
written down the contribution obtained from V^^^ . By choosing 

SvY^i^-^Tf, (5.25) 
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we ensure the finiteness of the matrix elements. The squared scalar mass is deter- 
mined by the zero of the determinant: 

Mj,,, = M^ + 4,..^^-^^^j^-^. (5.26) 

where I^{k) — I.,^{k) — Tf \ check (3.23). The mass matrix of the {X, Y^, s°) sector 
also becomes more transparent: 



/ -1 + 2iSxx + gih{k) 2i5xY» + y/9i92lAk) 2iy/giv \ 

2i5xYO + y/gml^k) -1 + 2i5YOYO + g^hik) 2i^v 
\ '^iy/giv '^iy/giv iDo^{k) J 



(5.27) 



which after introducing the obvious counterterms 



Sxx = yiT^'\ SxYO = ^V^rf , 5yoyo = ^^2Ti°) (5.28) 
leads to a determinant equation completely analogous to the previous one: 

o,,o) = + 4{g, + g^y- . ^ -^^ (5.29) 



l-(5i + 52)/^(A; = M(x,yo,,o))' 



We see, that the spectra resulting from the assumption we made for the mass hier- 
archy might be consistent with the outcome in the sense that the scalar sector (and 
the auxiliary fields hybridized with it) can be indeed heavier than the pionic one. 

In the next subsection we investigate if there is a region of the parameter space, 
where the scalar masses become much heavier than the pseudoscalars ensuring the 
self-consistency of our approximate solution. 



5.1.3 Validity of the heavy mass assumption 

First one has to note, that in the case when no explicit symmetry breaking term 
is added to the Lagrangian {h = 0) the mass assumption is true, since our approxi- 
mation preserves Goldstone's theorem and therefore makes the pions massless, which 
means that they are "infinitely" lighter than the scalars. 

The interesting case is when h ^ 0. Let us define rg y := M(s,y)/M and rx,yo .,0 := 
M(x,yo,sO)/M. In order to obtain a proper region of the parameter space, we introduce 
a heaviness criterion: the scalars are heavy enough if relations r^^y > tq, Tx.yo.s" > tq 
hold simultaneously, where tq is a given number. It is somewhat arbitrary what value 
to choose for Tq. In an exploratory study, we work with the convenient choice Tq = 2, 
because the quantity I^{k,M) appearing in both gap equations (5.26) and (5.29) 
develops an imaginary part just for > 4M^ (two-pion threshold). Above the 
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threshold the masses are defined as real parts of the complex solutions. Expressing 
V from the equation of state in (5.23), the two relevant equations are the following: 



2 _ 2 ^-f^I^iM^sY)) 



^M^ + An — l-gum^{M^x,YO,sO)) 



M^ = mV,:,(^ + ^log^-^^), (5.32) 



where gi2 = gi+g2 and 3?, 9" correspond to the real and imaginary parts, respectively. 
The finite part of the bubble integral was given explicitly in (3.23) and (3.24). With 
the use of (5.23) and (3.18), we get 

where Mq is the renormalization scale. It is convenient to express all masses in 
proportion to the absolute value of the renormalized mass {^/[^\). For a fixed Mq, 
(5.32) determines M'^ /\m?\ as a function of /i/|m|^, (71, (72- Plugging it into the scalar 
gap equations, they can be solved for M(s^y) / -^/jm^ and M(^x,yo,sO)/ \/\Tn?\- Then, 
we can trace out the region, where the heaviness criterion is fulfilled. This region is 
the part of the positive gi — g2 octant (the stability region of the n — )> 00 theory) 
below the surface displayed in Fig. 5.2. As expected, the projection of the allowed 
region onto the gi — g2 plane shrinks for increasing value of h/\m\^'. We have varied 
Mo/^/\m^\ in the interval (1, 10) and a mild displacement of the allowed region was 
observed. This change can be balanced by an appropriate renormalization group 
transformation of the quartic couplings, i.e. using g^ = gfj(Mo). 



5.1.4 Heavy scalar corrections 

Using the propagators of the coupled (s, y)-sector, we can start to systematically 
take into account the effect of the heavy degrees of freedom on the pion propagator, 
the EoS and the SPEs. In particular we now invoke the tadpole and bubble con- 
tributions to these equations evaluated with (5.24). For this we write explicitly the 
components of the heavy scalar propagator matrix: 

Gsrik) = -^fV^2V, 

d{k) = (1 - |/,^(A;)) {k' - M') - Ag2v\ (5.33) 

With the help of these expressions, we readily write the scalar tadpole contributions 
to the SPEs. The correction of the EoS is of imminent interest, since it has an 
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Figure 5.2: Region of the parameter space, where the mass assumption holds (using 



renormahzation scale Mq/ y jm^l = 1.0) 

important role in the discussion of the validity of Goldstone's theorem. Using the 
expression of G^y, we find the unrenormalized equation of state: 

M^ + 2gJ -^ = -. (5.34) 
Jp d[p) V 

Now we proceed with the unrenormalized pion propagators: 



J p 

iG-\k) = e-M'-g2 f 

J n 



i p"^ - M'^ 
d{j)) {k - pY - M'^ 



2 2 f ■ f {p'-M')il-g,I^{k-p)/2) ^^^^^ 

ldip)dik-p)-''^l diMk^) -^'-''^ 

When we set k = 0, both propagator equations turn out to be the same. Compar- 
ing to the equation of state we find, that the approximation fulfills the Goldstone 
theorem characterizing the U{n) x U{n) — t- U{n) symmetry breaking. It can be 
shown, that both equations receive the same counterterm contributions, therefore 
the theorem remains valid after renormahzation. 

Since the tadpoles of Gs and Gsy play a role in the EoS and the SPEs, it is 
worthwhile to investigate their corrections, although these contributions should be 
considered only as NNLO "heavy" corrections. Substituting the leading order prop- 
agators into the second term on the right hand side of the equation of Gjy (5.16), 
we find for the S^y self-energy: 

(5.36) 
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which is actually finite. The heavy correction to ^ is very similar to that of ^ 
in (5.35): 

.2 ,.2 f i P'-M' 



J V 



2„.2 



f i [ ip'-M-){l-g,I^{k-p)/2) 
J^d{p)d{k-p) '^Vp d{p)d{k-p) 



The scalar corrections are taken into account also in the SPEs. They appear 
partly directly via the scalar tadpole, and also by the scalar correction of the pion 
tadpole. The structure of the two SPE's in (5.13) is identical in view of (5.14), 
therefore it is sufficient to investigate the equation, which determines x. The un- 
renormalized form of the saddle-point equation reads as follows: 

= - - 2iv^ - I [ {G,{k) + G^{k)). (5.38) 

Qi Jk 

The expression of the scalar tadpole is readily written down. The pion tadpole is 
expanded to linear order in the self-energy contribution: 

-i ^ GM ^-^1 Do{k) + ^ ^k2^MY ^-' ^^-^^^ 



where 



^n{k) =92 

Jv 



p^ - Ag^v' ^ (p' - M^){1 - g^I^ik - p)/2) 



pd{p) [{k-py-M^ d{k-p) d{k-p) 

(5.40) 

One can achieve the explicit (resummed) renormalization of the equations (5.34), 
(5.35), (5.37) and (5.38). which will not be discussed here. In this procedure the 
remaining counterterms ^ttIq, ^^tji ^aiTD ^rrsj ^j/dj ^j/td ^i/s and dZ^^^dZg are determined. 
In a careful analysis, we find the following relations between them: 

^XV ^XTT 5yy ^yji 5yg. (5.4l) 

The reader can find the corresponding calculations in details in one of the papers 
(Phys. Rev. D82: 045011, 2010), whereupon this thesis is based on. 



5.2 Discussion of a more general breaking of the 
symmetry 

In this section we develop the previously introduced approximate large-n solution 
with a more general condensate. In this way we can investigate the stability of the 
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U (n) X U {n) — )■ U{n) symmetry breaking pattern. Also wc arc motivated by possible 
phenomenological applications announced in the beginning of the chapter. 

It turns out that for a more general description of the symmetry breaking, it is 
convenient to introduce not two, but three auxiliary fields. In the previous section the 
auxiliary field Y"" did not have a definite parity, however in the generalized analysis 
below it is worthwhile to work with quantities with a definite refiection property. 
Using the Jacobi identity, we easily recognize the decomposition 

where C/f = \dabc{s^ s"^ -\- t^^t^"^) and U2 — fabcS^'^'^- Exploiting this in the Lagrangian 
(5.1), every piece of it can be expressed with fields with definite parity. Then, [/" and 
U2 are perfect combinations in addition to + tt^tt" to be replaced by auxiliary 
fields. It will turn out that one of them mixes with the scalar (as in the previous 
section), and the other with the pseudoscalar elementary fields. The composite pseu- 
doscalar (i.e. corresponding to C/f ) is assumed to be similarly light as the elementary 
tt", therefore its dynamics will be included into our analysis, while the auxiliary fields 
corresponding to C/f and + tt^tt" are classified to be "heavy" , therefore their 
fiuctuations will be omitted from the contributions to the equation of state and 
the saddle point equations. In the present refined analysis, the mixing in the pseu- 
doscalar sector results in a spectrum, where also heavier excitations ( "kaons" ) might 
appear naturally. The self-consistency of this solution requires that even the heavier 
pseudoscalars should be considerably lighter than the lightest scalar excitation. 

5.2.1 Introducing three auxiliary fields 

Using (5.42), the Lagrangian (5.1) with a more general explicit symmetry break- 
ing term reads as: 

L = Jp^s'^)' + (aX)'-^'((s")' + (7r")')] + V2;^V 

- + + (5.43) 

where "q" indices refer to the diagonal generators of U (n). Then, three auxihary com- 
posite fields are introduced by adding the following constraints to the Lagrangian: 
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In the sum C = L + AL a general symmetry breaking pattern corresponding to 
diagonal directions in the algebra space is realized through the shifts 

+ n\/2dabqS^Vq + dapqVpVq. (5.45) 

The indices "q" correspond again to the diagonal generators. The index g = refers 
to the unit matrix and g = 1, n — 1 refer to the elements of the Cartan subalgebra 
of SU{n). 

After the introduction of the auxiliary field variables and the shifts, the full 
Lagrangian has the following form: 

1 



y Th 

'y^U!^ + nV2v,fa,cY^7^'] . (5.46) 




The background Vq introduces mixing of the pair s", Y"" for every value of the index a. 
In addition it induces nonzero values for Y^ and also a (2n + l)-dimensional mixing 
sector is being formed for the propagators of these following degrees of freedom: 
s'^,X,Y^. On the other hand pions will mix (for some index values at least) with 
the composite field Y^, which has negative parity. Since the ground state of strong 
interactions respects parity, we expect to have a quantum solution of the model, 
where Y2 = 0. Note that if Y2 had nonzero value, scalar and pseudoscalar fields 
would mix and the analysis of the mass spectrum would be a lot more complicated. 
Instead of going through the most general analysis, we assume from the start that 
7r° and fields do not mix, and below we show that this is consistent with Y2 = 0. 

We construct correspondingly a quantum solution in the same approximation as 
before (i.e. using a large-n approximated 2PI effective action and neglecting heavy 
excitations), where the saddle point values of X.Y^ are non-zero, Y2 = 0, and the 
2-point functions corresponding to the above mentioned coupled (mixing) sectors are 
assumed not to vanish. Our final goal is to construct the ground state of the system 
(which is determined by the equations of the 1-point functions), and since according 
to the heaviness assumption, in a first approximation only the quantum fiuctuations 
of the light (tt — Y2) fields enter to the equations of state of the condensate and the 
saddle point equations, therefore in this section we fully avoid the analysis of the 
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effect of the more complex symmetry breaking pattern on the scalar sector. We note 
however, that due to the fact that now three auxiliary fields have been introduced, 
the same approximation of the 2P1 effective action introdTiccd in the previous sec- 
tion leads to slightly different equations in the scalar sector than the ones found 
in (5.16) and (5.17). We will not reanalyze these objects as it was done there, but 
simply expect that the heaviness criterion remains valid in the essential part of the 
previously explored parameter space. For any specific phenomenological applications 
it would be recommended to investigate the correctness of this expectation, however 
in this study we limit ourselves to determine the features of the ground state of the 
system. 

The TT-sector breaks up into two separate groups of fields. The first group contains 
pions which mix with each other. From now on, we say that an index is diagonal, 
if it corresponds to a diagonal generator. Our first observation is that only pions 
with diagonal indices can mix with each other. This can be seen using some useful 
identities of the Appendix, or can be understood as follows. 

Since a potential mixing could arise from the term Y^dabc (i-e- Y^U^ in (5.46)), 
assuming that Y-^ 7^ 0, if a refers to a diagonal generator (this will be checked 
later), then the first index of dahc must be also diagonal. If the other two indices {b 
and c) are not diagonal, then they are necessary equal, since (check the Appendix) 
dabc ~ 9? Tr (A^A^A'^) (here A's are the generalized Gell-Mann matrices), therefore 
A^A'^ must be a real and diagonal matrix. From the generalization of the Pauli algebra 
(see the Appendix) this requirement is fulfilled only when b = c. This indeed shows 
that only diagonal pions can mix with each other. This mixing is determined by the 
specific pattern of the symmetry breaking and the corresponding structure constants 
can be found in the Appendix. 

The second group of pions are those which can mix with Y2. Here we find that 
only pions with non-diagonal indices can enter the game and each of these can mix 
with exactly one Kg"- This means that these arc disjoint from the previous pions and 
form 2x2 blocks. Again, using the identities for the structure constants given in 
the Appendix, one can immediately see this, however for the reader's convenience, 
we also give a short argument here. 

The explanation is very similar as the previous one, but now the mixing comes 
from Vafuav (sec the last term of (5.46)), where a 7^ is a diagonal index. It is easy 
to sec that independently of the value of a, we always get nonzero fuav for the same 
appropriately chosen pairs of u — v indices. In other words, if we have chosen a "m" 
index, for a nonzero fuav, we must point to the same "f" index independently of 
a. This is due to the fact (see the Appendix), that fuav ~ "^Tr ^[A", A^]A"j, which 
says that if A" is diagonal, then the commutator [A", A''] must be also diagonal 
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and purely imaginary (it is possible however that even in this case the matrix in 
the argument of the trace gives zero). This immediately fixes both u and v, since 
they must correspond to the non-diagonal generators, which are the generalized 
Pauli matrices (i.e. generators with indices {x,jk), {y,jk), for the notation see the 
Appendix) . 

We introduce notations for the propagators of fields described above. For the 
diagonal pions we have a huge n x n matrix, denoted by = G^^aT^b. The 

same quantity for the non-mixing ¥2^8 are (^yj'^^)^^ = Gyayb- The inverses of these 
at tree level multiplied by i are 

ziVt'-yJip) = (/ -m' + tJ%X)8a, + t\[^da,cY^ (5.47a) 

The 2x2 blocks of non- diagonal pions and the corresponding Y2 fields are denoted 
by: 

where u ~ {x,jk) and v = {y,jk) or vice versa. Since there are n(n — l)/2 possible 
choices for u and v, taking also account the possible interchanging between them, 
we have a number of n(n — 1), in general different sectors of this type. The inverse 
matrix at tree level multiplied by i is: 



\ W^g^nfvQuVQ -1 



(5.49) 



We see that sectors connected via the interchanging of the indices are almost the 
same, only the off-diagonal elements differ due to a minus sign (this comes from the 
total antisymmetry of the structure constants /^q^ and the fact that u and v are 
not independent, therefore dauu = daw)- 

We needed the expressions of the inverse tree level propagators, since they ex- 
plicitly appear in the 1-loop part of the 2PI effective potential. As it was already 
stressed, now we would like to focus only on the pseudoscalar sector, therefore pieces 
corresponding exclusively to the scalar sector are not present in the expression below. 
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The approximate 2PI effective potential witfiout counterterms looks as: 



- ^Ei^log(^i^))-^-l^IVlog(^-)-^ 

(u,v) 

- ^ / IV log fc^-)" -iY.h^ ip^r'G^ 

(p + k) 

J p J k 

Jp Jk 



n 



^ ■ I '^9x92 



. .92 

+ I — 
n 



Gn„^^ip)G,.Yb{k)GYa^?ip + k) 
i^J^^daa/^S^s J J^G^^^^{p)G„g^^{k)GxY-{p + k) 
^^afiS^s J J^Gxx{p)G^^^g{k)G^^^^{p + k) 
+ i^daapdbys J J GYaYb{p)G^^^^{k)G^^^^{p + k) 



+ i 



. ■ 91 
+ I 



' p J k 
Ip Jk 



+ i^faiafbjp I I Gs^s,{p)GY-nf,{k)GY-nAP + k)- (5.50) 



5.2.2 Equations in the general framework 

Prom the effective potential (5.50) we can derive equations for the 1-point func- 
tions and for the propagators of the 7r — Y2 sector. The last seven hues correspond to 
the setting sun diagrams which give contributions to this sector. In these expressions 
we did not simplify the index notations as we did in (5.11), since we will not make 
use of these expressions in the actual equations. This is due to our heavy scalar 
assumption: checking the setting sun diagrams in (5.50), we see that they contain at 
least one scalar (heavy) propagator, therefore in a first approximation the quantum 
corrections of the pseudoscalar propagators can be neglected. It means that the 2- 
loop truncated effective potential together with the heavy scalar assumption leads 
to tree level propagators in the pseudoscalar sector: 

^diag ^ pdiag ^ ^diag ^ ^diag ^ gM ^ _ ^^^^^ 
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Now let us work out the equations for the 1-point functions. The saddle point 
equations of the auxiliary variables are: 

0^^-^^X-i./2^,{v''f-z^l^J^G^a,a, (5.52a) 
= = ^1 - in^/g^dabcVbVc - \\j^^<^bc J G^b^c, (5.52b) 

Note that for Y2 we obtained a trivial relation, since only scalar-pseudoscalar mixing 
propagator G^a^b could have given contribution, which was previously excluded from 
the solution on the basis of the parity conservation of strong interactions. This shows 
that the Y2 — assumption is consistent. Analyzing (5.52b) we obtain that for non- 
diagonal "a" indices, the right hand side is identically equal to Y"". This means that 
Yi inherits the structure of v, i.e. only diagonal indexed field expectation values are 
non-zero. For the equations of state we obtain 

= ^ = 2n'{m\a - K) - 2i./%Xn\« 

-2in^/ng^dabcVbyi - iV^^fvau / Gy-n^ip)- 

Jp 

It is convenient to treat the zeroth index independently. For the saddle point equa- 
tions we get 

X = i^/2^l{vl + vD + J {Gno^roip) + G'^q^q(p)), (5.53a) 

y° = iV2^2{vl + vl) + {G^oM + G.^^^ip)), (5.53b) 
Y^ = i^^^(AnV2nvoVQ + 2n^dQBcVBVc 

+dQuu J <^7r„7r„(p) +2^^ j G^^^^ + dgBC J G^B7rc(p))- (5.53c) 

Here capital letters refer to nonzero diagonal indices, while "w" corresponds to non- 
diagonal ones. With the use of the previously introduced definition of (check 
(5.9)), the equations of state simplify to 

= V2n{ho - A'Pvo) + 2t^Y^^VQ, (5.54a) 
= V2n{hQ - M^vq) + 2i,/^Y^^Vo + i^/2^dQBcyfvc 

+i^l^fvQu j Gy^'Kuip). (5.54b) 
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In order to solve these equations and to determine the ground state of the model, 
we need to obtain the propagators Gy^n^- We saw that our approximation led to 
tree-level propagators in the pscudoscalar sector, therefore we only need to invert 
the 2x2 matrices given in (5.49). We obtain 



- ^Lv) XV^fubvVi, (p' -m^ + 2i^X)/i + ^nd. 




where 



M^u,v) - - 2^y - ^dauuY^^ - 2g2n{UuVbf, (5.56) 

which is symmetric under the change u v (recall that dauu — daw)- 

Before we choose a specific symmetry breaking pattern, we have to analyze the 
scaling behavior of the quantities, since we would like to work in the large-n approx- 
imation. Let us therefore introduce the following rescaled variables: 

XQ^V^VQ, yf^Y,^/V^, x^X/n, yl^Y^/n. (5.57) 

Note that the definitions of x and y1 somewhat differ from the ones introduced in the 
previous section, check (5.21). Using these notations, we approximate our equations 
at LO of the large-n limit. For the equations of state we have 

= ho- M^VQ + i^2y'^XQ/n, (5.58a) 

= /iQ - M^xq + i\[^2Vx'"Q + i-JlhdQBcVxXc 

-g,Y^h9^^^T{Ml^^^). (5.58b) 

Note that we have also rescaled the diagonal indexed external fields: hq hq / ^Jn. 
With the new variables can be written as = — i{^J2g[x -\- ^y2g^y^) and 
T{fjp) — Jpi/{p^ — /U^) is the usual tadpole integral. 
The saddle point equations read as: 

= iV2i2(vl + XAXA/n + ^ ^ ^(^(l^))) ' (5-59a) 

{u,v) 

y? = i\^(^'2V2voXQ + dQBcXBXc 

+^dQBC J G^B^cip) + ^ J] dQuuT{Ml^^))j , (5.59b) 

P {u,v) 

with ^x = y/gly^. 

From now on, we choose a symmetry breaking pattern, where vq = VsSqs- The 
index "8" refers to the longest diagonal gcncrator(!) of the U{n) group, reminding 
to the n = 3 case. We stress that this is just a notation, our analysis is a large-n 
study and has nothing to do with the U{3) case. 
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5.3 Ground state of the system with a breaking 
in the 0-8 sector 

First we discuss the mixings when the only spontaneously broken direction be- 
yond the zeroth one is the "8", corresponding to the longest diagonal generator. 
Recalling (5.47a), we assume that Yi inherits further the structure of v (it will be 
checked later) and therefore — Y^6ao + Y^Sas- Then we obtain 

K'^-)ab(P) = - ^^-b + tV9~2dsabyl (5.60) 



ab 

Using the Appendix, when A, S 7^ 0, 8 are diagonal indices, then 



dsAB = \l —, T\^AB- (5.61) 

" n[n — 1) 



Also, 




(5.62) 



This shows that the pion mixing is reduced to the — 8 sector. 

Now we turn to the non-diagonal pions, which can mix with Y2. Recalling (5.49), 
we need /„8t) 7^ (since /„ot, = 0). Using the Appendix, the only possible choice is 
to have u — {x,jn) with v = {y,jn) or vice versa. 

We have (n — 1) choices for the value of j, and since the interchange oi u -ir^ v 
doubles the number of nonzero structure constants, therefore we have 2{n — l) pions 
which mix in the large-n limit with the corresponding Y2 fields. These pions will be 
called the kaons. Using the Appendix, independently of the u, v values in question 
(i.e. (xjn), (yjn)), we obtain 



/.8.--^y^. (5.63) 

The dsuv values for the non-diagonal indices are as follows. For indices which do not 
count in the mixing with Y2 (i.e. u,v — {x,jk), {y,jk) where k ^ n): 



dsuv — \ 1/ tt'^uv, (5.64) 

" nin — 1 



and for the ones which mix (i.e. u,v — {x,jn), {y,jn)): 



(n — 2) 

dsuv — / , =Suv (5.65) 

^y2n{n — 1) 
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This shows that apart from the kaons, the "88" and the "00" pions, every other pions 
are identicaL With the help of (5.61) and (5.64), we find (a, b ^ {x,jn), {y,jn), 8, 0) 

KS-Va^'iP) = iP' + ^—yl)Sab ^ ip' -M'- -AM^)Sab, (5.66) 

where we have written down the NLO correction to the mass, which is needed for 
the NLO construction of the effective potential (see later). RecaUing (5.55), for the 
kaon-la sector we obtain: 



(5.67) 



where 



Ml^M' + t^^yl + g2xl (5.68) 

Note that we have changed tt to X in the notation of the propagator matrix, and 
erased the indices of it, these should be thought as appropriate {u, v) pairs regarding 
to the previous analysis. The ± sign in the off-diagonal elements signals the doubhng 
which corresponds to interchanging oi u v indices. This however has no effect in 
the mass M^, as it was already stated in the general analysis. 

Now let us see how the equations of state and the saddle point equations simplify. 
Using the rescaled variables introduced in (5.57), sticking to the large-n hmit, from 
(5.58) and (5.59) we have 






= ho 


- M\o, 


(5.69) 





= hg 


- M\s + iV^2y\{vQ - Xs) - g2XsT{Ml), 


(5.70) 


X 


— i\ 


f2iM + \T{M'')), 


(5.71) 




— i\ 




(5.72) 


y? 


— "^V 


%{2V2voXQ + dqssxl) 





+ ^ ( E dQuuT{Ml.)) + J2 dQAAT{M')) . (5.73) 

{u,v) A^8 

Again, capital letters {Q, A ^ 0) correspond to diagonal indices. First we should 
prove that inherits the structure of v, i.e. only Q = 8 gives nontrivial equation. 
The first two terms of (5.73) are obviously zero for Q 7^ 8, therefore we have to 
analyze the last two. In the third term, almost every pion has mass, only the 
kaons are exceptions. In the last term, if we include A — 8 into the sum, the error 
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we commit counts only at NLO. In this case we have 

{u,v) A 

= dQ^^{T{Ml)-T(M'))+J2dQaaT(M'),(5-7^) 

u = kaon a 

where in the second sum of the right hand side the summation goes over every 
possible index. Since Ylia ^Qaa — for every Q, this last term cancels and we must 
analyze the sum X^^j = ka,on dquu- Let us assume that Q ^ 8. Using the Appendix, we 
can easily obtain that 



dquu — 0. 



(5.75) 



u = kaon 

This proves that yf is nontrivial only for Q = 8. In the latter case the sum is 



dsuu ^ -{n - l){n - 2) 

= kaon 



n(n — 1] 



therefore the equation of yf is 



yf = i^/2i,{2voXs -Xl- \nMl) + \t{M^)). 



(5.76) 



(5.77) 



Wc can solve the equations (5.70), (5.77) for the "8" components with fo, y^, x 

taken from their respective equations (5.69), (5.71), (5.72) for a set of h^, gi,g2- This 
latter part of the calculation coincides with the solution presented in the first part 
of the chapter where only fQ-condensate was considered. 

Equations (5.70), (5.71), (5.72) and (5.77) are divergent, therefore need rcnor- 
malization. We would like to have the very same equations, but with the tadpole 
integrals substituted with their finite part in the renormalization scheme we used 
before. In this one T^(/i2) = T {^i^) -Tf^ - {^i'^ - M^)T'f\ check (3.15) and (3.18). We 
can apply the method described in details in Chapter 3: identify first the divergent 
pieces in the equations, then integrate back to obtain the appropriate counterterm 
functional of the 2P1 effective potential. Without going into the details, wc just state 
that every divergent equation can be made finite in the previously mentioned scheme 
with the following counterterm functional: 

Vct^i- / ln(V + Mo') 



n 
— 11 



1 

2' 



^ 2 



(M2-m^)(rf -Mo^rf 
^AM^ + Ml - M^^ {rf^ + (M^ - Ml)T^ 



(Oh 



(5.78) 
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Although (5.78) is written in a rather compact form, it is a function of vq, xs, Hi, yi,x 
and is completely allowed by the formal requirements of renormalizability [44] . Note 
also the important simplification + Mj^ — Af^ = g2X8- The first term of 

(5.78) is field independent, therefore gives no contribution to the equations for the 
1-point functions. It serves only for eliminating the overall divergence of the effective 
potential generated by the zero point fiuctuations of the pions. 

Plugging (5.51) to (5.50), we obtain the IPl effective potential. The obtained 
functional depends also on the auxiliary variables and as it was already stressed, we 
include only the quantum fluctuations of the light particles (i.e. pions and kaons). 
Since every setting sun in (5.50) contains at least one heavy scalar line, we neglect 
these diagrams. The corresponding formal expression is the following with 0{n) 
accuracy: 



V 

Vcl 



= Vcl + Va 



n 



+n 



quant i 



V, 



quant 



{n'-2n) / In (-p' + M;) + 2n / In (-p' + M^) 
Jp Jp 



(5.79) 



where Vquant is the standard expression of the one-loop part of the effective potential. 
It is important to note that for the accurate evaluation of the quantum contribution 
of the pions, we have to make use of the 0{n) accurate value of the pion mass, that 
is Ml — + ^AM^, see (5.66). The divergences of Vquant are ehminated by Vd- 

The renormalized effective potential {Vr = V+Vd) in the renormalization scheme 
defined through specified in (5.78) above looks hke 



V, 



R 



n 



+ n 



M'vi + -{x' + {y',)')-2hoVo 



+ 



+ 



n 



327r2 

n 
167r2 
1 



M'xi + ^iVi) - 2/i8X8 - W'2g2ylXs{'2vo - Xs) 
'-{M^ - M') + Ml{M' - Ml) + -M' In 



1 



(AM^-M^m^ + lM^m^ 



^2X8^0 



-AMJm2 + -(M^ 
2 V 



M'k) 



(5.80) 



The fact that the "8" condensate contributes only to 0{n) means that the fg contri- 
bution will be subleading relative to that of Vq. Still it belongs to the leading large 
n piece, since the NLO contributions which were omitted from the very start would 
contribute 0{vP). 
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5.3.1 Numerical results 

We have explored the gi, §2 plane and in a large region we found three solutions 
for xs, y\ in a fixed background of Vo,x,y values. It turned out that one of them is 
the trivial solution (when hs = 0), and the other two represent a local minimum 
and maximum of the effective potential. In order to select the true ground state of 
the system between the trivial and non-trivial local minima, we must look for the 
energetically more favorable one. It is clear that we must investigate only the Xs, uf 
dependent parts of the effective potential. We can immediately see that only the 
0{ri) part of Vr is what we have to study. This part coincides with the set of the 
(XsjUi) dependent terms. It is important to stress that due to this fortunate effect, 
the selection of the ground state does not depend on the explicit value of n, since 
it appears as an overall multiplicative factor. In the forthcoming figures we refer to 
y^"^ as the 0{ri) part of the full effective potential divided by this overall n factor. 

In the first part of this chapter we mapped out the region of the three-dimensional 
parameter space {gi,g2,hQ), where the scalar masses turned out to be at least twice 
as large as the mass of the pseudo-Goldstone fields (pions) . In this part of the inves- 
tigations, relative to the case of the pure Vq background, we have chosen rather low 
values of hg (all dimensional quantities are measured in proportion to the appropri- 
ate powers of the absolute value of the renormahzed mass). The results below are 
obtained in those regions of the parameter space, where we found the heavy scalar 
assumption valid (i.e. where the scalar masses are at least twice as large as the pion 
mass) . 

Let us start the discussion with the mass splitting in the vr-sector, in particular 
the kaon/pion mass ratio. The most detailed study was achieved for /?.8 = 0, ad- 
dressing the question of the vs (i.e. X8)-dependence of the effective potential after 
eliminating the yf dependence of Vr with the help of (5.77). We scanned through 
the allowed region at different values of Hq. As announced, in addition to the trivial 
solution Vg, = yf = 0, also a non-trivial (meta) stable solution (local minimum of the 
effective potential) of (5.70) and (5.77) with positive squared masses of all 7r-fields 
was found in a large part of the region allowed by the heavy scalar mass assumption. 
In Fig. 5.3 the ratio M'^/M^ is displayed over the (gi, g2)-p\a.ne for 2 different values 
of hf). We see, that in a quite large region of the allowed parameter space the mass 
ratio varies mildly, but the kaons arc getting heavier when ho is lowered. Also it 
can be observed, that for small gi values increasing g2 induces very large kaon mass 
values, while for larger gi the same does not lead to significant changes of the mass 
ratio. 

The 0{n) part of the full renormahzed effective potential Vr could be plotted in 
an extended interval around the non-trivial minimum, where all mass-squares are 
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Figure 5.3: Kaon/pion squared mass ratio as a function of the couplings {gi,g2) for 
two different ho values. Kaons become lighter as ho increases, at the same time the 
region where non-trivial solutions are found narrows. 



g,=5.5, 92=9.0, M„=|in| 




Figure 5.4: The 0{n) part of the quantum effective potential divided by n as a 
function of the xs condensate. The real part of the potential is plotted with dotted 
lines for xs values where some negative mass squares lead to complex values of the 
potential. The curves from top to bottom follow the order of the ho values indicated 
in the upper- left-hand corner. 
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positive. An example is shown in Fig. 5.4, where we see that the energy density of the 
non-trivial solution grows monotonically with Hq increasing. (As already mentioned, 
the 0{n) part of the full potential divided by n is denoted as \/("), and the same 
holds for the classical part as V^^K) There is a critical value /io(5'i, 5'2), where the 
non-trivial local minimum disappears. For example we can verify on Fig. 5.3 that for 
the gi,g2 values of Fig. 5.4, there does exist a non-trivial solution of the equations 
for ho = 0.005 but not for = 0.05. This circumstance explains why we used much 
smaller values of ho compared to the first part of this chapter. 

The dotted pieces of the curves represent the real part of the effective potential 
in regions where it becomes complex due to the fact that one of the mass-squares 
becomes negative. This phenomenon is a well-known feature of the loop expansion 
[4, 110, 111]. Although the effective potential should be a real quantity by construc- 
tion, in regions where the classical potential is chosen to be non-convex, the loop 
expansion breaks down and produces an imaginary part. Mathematically, a non- 
convex classical potential can be understood as an analytic continuation of a convex 
classical potential through the mass term: m im. The loop-expansion creates a 
divergent asymptotic series, but there is no such theorem which would state that 
the analytically continued series corresponds term-by-term to the series of the ana- 
lytically continued effective potential. In other words, the analytic continuation and 
the expansion in h are not interchangeable operations. From a mathematical point 
of view, this is the origin of the appearance of the imaginary part. Nevertheless, 
the naive calculation leading to complex effective potentials has a physical meaning. 
First, it turns out that if we want to understand the effective potential at a 0c point 
as the expectation value of the minimum energy density of the class of quantum 
states in which the field expectation value is (pc (which property can be derived from 
the usual Legendre-transform definition of the effective potential) , then we have to 
take the real part of the loop-expansion and apply Maxwell's construction. Further- 
more, Weinberg and Wu showed [112] that even the naively calculated real part in 
itself contains physics: it is the minimum energy density of a class of quantum states 
with a (pc field expectation value together with the more restrictive constraint that 
their wave functionals are concentrated on the uniform configuration (pc- The state 
chosen this way may differ from the state which minimizes the energy density of 
states with a (pc field expectation value. In this case the state in question is unstable 
and the imaginary part of the potential can be interpreted as half the decay rate per 
unit volume. 

We have also plotted the curve of the classical part of the 0{n) piece of the poten- 
tial in a typical point of the coupling space for two values of ho, which demonstrates 
that with the present renormalization scale Mq the quantum fluctuations only mod- 
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h(,=0.005|m|3, g,=5.5, g2=9.0, M(,=|m| 



ho=0.01|m|3, g,=5.5, g2=9.0, Mg=\n\\ 
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0.008 
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Figure 5.5: The effect of the quantum ffuctuations is only moderate, as indicated by 
comparing the V^^^ (dashed hues) with V^"'\ The dotted part displays the real part 
of V^"'\ where it becomes complex. 



erately modify the value of the potential, see Fig. 5.5. In the whole allowed region 
of the coupling space we found for the potential in the non-trivial minimum positive 
values, therefore we conclude that with hg = 0, there is no /io,fo,M^ combination 
which would induce spontaneously nonzero vs condensate into the true ground state 
of the system. The effect of the existence of such a metastable state might still show 
up in out-of equilibrium situations. 

In a number of energetically "promising" coupling points we have studied the 
effect of the application of a non-zero on the effective potential. In these cases 
also the trivial solution shifts in proportion of the external field, and the effective 
potential can be displayed in a finite, but small region around the shifted trivial 
minimum too (i.e. no negative mass squares appear in its neighborhood). If Hq is 
small, the value of the effective potential in the non-trivial minimum varies faster 
with /ig than near the shifted trivial minimum, therefore we find a critical external 
field hl{gi,g2), which drives the minimum characterized by a large vq value to be 
the true ground state of the system. In Fig. 5.6 we show the typical variation of the 
potential as a function of hg for the value of {gi,g2, ho) chosen as in the left panel of 
Fig. 5.5. In Fig. 5.7 we display the variation of the ratio Mj^/M^ as a function of hg 
and g2 for given gi and values as taken in the non-trivial minimum. We see that 
apart from a jump to unity happening at lower g2 values, the ratio is not varying 
too much. This ratio equals unity in the region where the non-trivial minimum is 
missing. Choosing higher values for involves smaller kaon masses (as it can also be 
seen in Fig. 5.3), which points to the direction in the parameter space with a chance 
of phenomenological applications. We also note, that choosing other gi values, the 
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ho=0.005|m| , g,=5.5, 92=9.0, Mo=|mI 



v'"'/|m|'' 0.01 



h8=0 
h8=0.004 |m|3 
hg.O.OOS Imp 
h„=0.012 Imp 




Figure 5.6: Smooth variation of V^"'^ with increasing hs displaying the symmetry 
breaking pattern Uv{n) — )■ Uv{n — 1). The curves from top to bottom follow the 
order of the applied hg values indicated in the upper-left-hand corner. 




Figure 5.7: The variation of M'^/M^ with ho and g2- The relative variation of the 
mass ratio is very small with respect to hg. Note that the abrupt jump downwards 
is occurring when the non-trivial minimum is missing. 
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shape of the figure changes only mildly. The figure also shows that the mass ratio 
is increasing linearly with hg at fixed g2, however this change is almost negligible 
compared to its full magnitude, i.e. the relative change is very small. This means 
that the kaon mass is induced dominantly by the action of the vq condensate, and not 
by the explicit symmetry breaking proportional to hg. This circumstance is rather 
different from what comes in the n = 3 flavor models from various perturbative 
discussions. 

5.4 Concluding remarks 

In this chapter, a leading order approximate large-n solution of the U{n) xU (n) 
model was presented in auxiliary field formalism. In the first part of the chapter a 
certain symmetry breaking pattern was assumed, which reduces the U{n) x U{n) 
symmetry of the Lagrangian to U{n). During the construction of the solution, a 
light pseudoscalar/heavy scalar hierarchy of the spectra was assumed. The renor- 
malization of the saddle point equations, the equation of state and the propagator 
equations was performed with explicit counterterm construction, which was followed 
by an investigation of the consistency of the additional heavy scalar assumption. It 
turned out that in a large region of the parameter space the existence of heavy scalars 
occurs. It was also shown, that the proposed solution explicitly fulfills Goldstone's 
theorem. 

In the second part of the chapter we examined the possibility of a more general 
symmetry breaking, in particular when the remaining U{n) symmetry is reduced 
further to U{n — 1). In this case new solutions were found. After renormalization, it 
was shown that in the direction which corresponds to "8" in the 5*6'^ (3) classification, 
in addition to the trivial minimum, a mctastable local minimum of the renormalized 
effective potential appears, even without introducing explicit symmetry breaking in 
the corresponding direction. This means that a condensate can arise spontaneously 
in the direction in question, which breaks the remaining U (n) symmetry to U{n — 
1). It turned out that the application of a moderate external field transforms the 
nontrivial minimum into the true ground state (i.e. the global minimum of the 
effective potential) of the system. 

The different ground states arc characterized by different pscudoscalar excitation 
spectra, also explicitly constructed. The mass ratio of the heavier modes ("kaons") 
and the remaining pseudo-Goldstone modes was found sensitive mainly only to the 
external field in the zeroth direction (which is proportional to the unit matrix), 
i.e. to the symmetry breaking pattern of the starting solution. If it turns out that 
phenomenology can be built on the nontrivial vacuum condensate, this observa- 
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tion would imply also rather different finite temperature befiavior compared to tlie 

perturbativc descriptions. 

Altliougli tlie transition of tlic condensate in tlic zerotli direction in itself would be 
mucfi similar to that of the 0{2n^) symmetric model, the insensitivity of the location 
of the new nontrivial minimum with respect to the external field in the "8" direction 
(see Fig. 5.7) hints a first-order transition. We observed that it is not the external 
field but rather the condensate in the zeroth direction which determines the strength 
of the symmetry breaking to U{n — 1). The tightly connected variation of the two 
condensates probably leads then to a phase transition with both fields going through 
a discontinuous change of similar amplitude. This expectation complies well with 
renormalization group considerations which classify the restoration of the broken 
symmetry in this model quite differently than in models displaying symmetries of the 
orthogonal group [76, 77]. An obvious next step is the study of the finite temperature 
symmetry restoration from this nontrivial ground state, but this is not included in 
this thesis. 

Further consolidation of the discovered structure should require the investigation 
of the fluctuation effect of the heavy scalar fields. On a general basis, one expects 
the modification of the different coefficients in different EoS, SPEs, and self-energies 
by terms of 0{M^/Mg) (where Mg is one of the scalar masses), which would just 
slightly deform the coupling region where the nontrivial vacuum structure is self- 
consistently present. 

A final goal could be to make use of this solution at n = 3, when also the 
contributions from the effective term reflecting the axial anomaly (the 't Hooft de- 
terminant) would be included perturbatively using the large-n propagators of the 
different fields. The comparison with the treatments of the n — 3 three flavored me- 
son model in different perturbative calculations [81, 82, 84, 113, 114] and functional 
approaches, which go beyond perturbation theory [118, 119], partly based on the 
exact renormalization group techniques [115, 116, 117], could give some insight into 
the importance of strong coupling effects. 



Conclusions 



In this thesis we apphed 2PI formahsm to various scalar quantum field theories, 
at zero temperature. The main scope of the work was to investigate the renor- 
malizability of the partially resummed perturbative series induced by different 2PI 
approximations. Beyond the theoretical realization, the study of the numerical fea- 
tures of 2PI renormalization was also a major motivation of this work. In a specific 
theory [U (N) x U (N) model) we also managed to present new ground state solutions 
of the 2PI equations, which could lead to phenomenological applications. 

First we dealt with the 0{N) symmetric model in auxihary field formulation. We 
constructed the 2 PI effective potential at next-to-leading order of the 1/N expansion, 
and showed its renormalizability with explicit counterterm construction. We found 
that the renormalization can be done without any restriction on the background 
fields, stated earlier in [37, 38]. We also eliminated the auxiliary field to get the 
next-to-leading order renormalized 2PI effective potential in terms of the original 
variables of the model. During this step we could see in a transparent way how the 
partial series of counterterm diagrams are resummed in the original formulation of 
the model. The renormalization of the 2PI effective potential formulated with an 
auxiliary field led to renormalized functionals after the elimination of any subset of 
the variables. 

Wc also discussed algorithms of the numerical implementation of 2PI renormal- 
ization. We employed the one component 0'^ theory at 2-loop level of the 2PI effective 
potential. Our main goal was to implement very accurate numerical realizations of 
the renormalization procedure, and to quantify the rate of convergence of the numer- 
ical solutions. In order to develop different algorithms, we used different theoretical 
approaches to obtain renormalized equations. We pointed out that counterterms of 
[50] are equivalent with appropriate renormalization conditions imposed on the 2- 
and 4-point functions. We developed accurate numerics regarding both variants of 
the renormalization programme and found the very same result numerically, which 
confirmed our theoretical considerations. The method using the exact counterterms 
turned out to be inconvenient at large coupling constants: solving the equations 
required large computational resources. We managed to overcome this with a new 
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algorithm which did not work with fixed coiintcrtcrms, but iteratively developed 
them together with the equations. This led to a significant improvement in the con- 
vergence rate, therefore the method is expected to be a valuable tool for solving bare 
(i.e. counterterm included) equations of strongly coupled theories. 

In the last chapter, we worked on the solution of the U (N) x U (N) meson model 
in the large-A^ limit. Wc built up the 2P1 effective potential at 2-loop level using 
auxiliary fields. We presented a solution of the coupled set of equations of the 1- and 
2-point functions in the broken phase of a one component condensate realizing the 
U{N) X U{N) — > U{N) symmetry breaking. Our solution involved an additional 
assumption on the scalar/pseudoscalar masses: the scalars are considered to be much 
heavier than the pseudoscalars. Due to this, we neglected terms (diagrams) of the 
equations, which contained at least one scalar propagator. This approximate solution 
was the first large- iV study of the model, which went beyond the Ansatze related to 
the 0(2iV^) symmetric piece of the Lagrangian of the model. The renormalizability 
of the solution was presented with explicit counterterm construction and we showed 
that it fulfills Goldstone's theorem. We also solved the scalar propagator equations 
numerically. With this we demonstrated that in a large part of the coupling space, the 
heavy scalar assumption is valid, confirming the self-consistency of our approximate 
solution. 

We also discussed the stabihty of the symmetry breaking with respect to a more 
general condensate. After deriving equations in the same approximation for a general 
breaking of the symmetry, the study of the renormalized effective potential showed 
that the field equations have further solutions in certain part of the parameter space. 
One of them turned out to be a new local minimum of the effective potential, which 
is a metastable vacuum even without the existence of an explicit symmetry breaking 
term. However, moderate application of a conjugate source, the metastable vacuum 
transforms into the global minimum, i.e. the true ground state of the system. This 
nontrivial structure of the effective potential is a signal of a more completed treat- 
ment of the model in the sense that, our solution dealt with the full U{N) x U{N) 
structure, not just the partial 0{2N'^) nonlincaritics, which had been a typical ap- 
proximation before. Our analysis hinted a first order finite temperature transition in 
the cases where the nontrivial minimum is the true ground state. This is expected 
on general grounds of the renormalization group analysis in the large- limit [77], 
however until now no publications are known which could constructively reproduce 
this conjecture. 



Appendix 



structure of the U(N) algebra 

The U{ri) algebra is dimensional, the generators are denoted by T". They are 
traceless and normalized in a way that Tr(T"T''') = Sab/'^- The generalized Gell- 
Mann matrices A" are defined as 



rpa 



(Al) 



therefore Tr (A" A'') = 26ab- The generators can be classified into two sets: there are 
diagonal, and non-diagonal ones. There are n diagonal Gell-Mann matrices: 



/I 



A° 



/ 1 



A^ 



1/ 



n{n — 1) 



/ 1 



v 



-(n-1) J 



(A2) 



We note that in the Chapter 5, we use the "8" index to A" ^. The number of the 
non-diagonal ones is n{n — 1) and they read as 



(A3) 



where we introduced the compact index notations {x,jk) and {y,jk) {j < k). This 
arises from the fact that the matrices defined in (A3) are related to the ax and ay 
Pauli matrices: they are mainly the same but the elements of them appear in the 
j — k row- column. 

The structure constants of the algebra are defined through the relation 



A" a'' — {dabc + ifabc)^'^. 



(A4) 
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where dabc is totally symmetric and fabc is totally antisymmetric in their indices. 
From this we have 



{X%X'} = 2dabcX'' 



Tr[[A^A^]A^] =AtU„ (A5a) 
Tr[{A«,A''}A<=] =446c- (A5b) 



Here [., .] and {., •} refer to commutation and anticommutation, respectively. Alter- 
natively, we can write 



fabc " 
dabc 



$5Tr(A"A^A"), 



sftTr(A"A''A"). 



(A6a) 
(A6b) 



We can immediately see that 



^Oab 



and in can easily be proved that 



Kb, foab — 0, 



(A7) 



(A8) 



Now we list some useful identities which can be derived from the properties (A5) 
or (A6). We assume that Q 7^ corresponds to a diagonal generator. Then, for the 
antisymmetric structure constants the following can be derived: 



Q 



2(Q+i)' 
1 



f{x,jk),Q,{y,jk) — * 



y^2Q(Q+l)' 
2Q ' 



[0, 



j < Q + 1 < k, 
j <Q + l^k, 
j <k<Q + l. 



(A9) 



A special case when Q — n — 1, that is "8", then only the third and fourth lines 
make sense. The nonzero values are: 



f{x,jn),8,{y,jn) — f{x,jn),n-l,{y,jn) 



We also note that when j ^ j' and k ^ k\ 



n 



2(n- 1)' 



(AlO) 



f{x,jk),Q,{x,j',k') - f{y,jk),Q,{y,j',k') - 0, 
f{x,jk),Q,(y,j',k') = f{y,jk),Q,{x,j',k') = 0. 



(Alia) 
(Allb) 
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For indices corresponding fully to diagonal generators, the symmetric tensor looks 

like (g,p,ry^o): 



dqpT — 



T(T+1) 



Q(<9+i) 



Q<P<T, 
Q <P = T. 
Q^P<T, 



Finally, we have (Q 7^ diagonal) 



Q,{3:,jk),{x,jk) — dQ,{y,jk),{y,jk) 



Q 

2(Q+i)- 



^2Q(Q+1) ' 
y/2Q{Q+l)'' 



and 



(A12) 



Q + 1 < j < k, 

Q + l=j <k, 

j<Q + l<k, (A13) 

j <Q + l = k, 

j < k < Q + 1, 



dQ,(x,jk),{x,j',k') - dQ,(y,jk),{y,j',k') - dQ,(yjk),{x,j',k') - 0. 



(A14) 
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